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PREFACE 



This volume presents the description and program documentation of 
the three-dimensional, free-surface mathematical model for thermal pollu- 
tion analysis and prediction for shallow water bodies, for example, lakes 
and coastal waters. The program was developed by the Thermal Pollution 
Croup at the University of Miami, and was successfully verified through 
application to several sites. This success was made possible by funding 
and technical assistance provided by the National Aeronautics and Space 
Administration (NASA) and the Environmental Protection Agency (EPA). 

The model is time dependent, and the leap-frog and DuFort-Frankel 
schemes are adopted for solving the predictive equations based on the 
conservation principles of mass, momentum and energy. The model has 
been developed with minimal physical and site restrictive assumptions, and 
its algorithm has sufficient generality to allow for different boundary con- 
ditions specified at open boundaries. The program shows both the tem- 
poral and spatial variations of the surface water height. It computes 
three-dimensional velocity and temperature fields. The model can serve 
as an effective means for hydrothermal analysis and prediction. Plotting 
programs employed for representing the numerous results are also in- 
cluded* 

The volume is intended as a user's manual and, as such, presents 
specific instructions regarding data preparation for program execution. 
To illustrate further, an example case is included here with its input 
data, hard copy printout and plots. The complete listing of the program 
and its accessories is also included. 



ABSTRACT 



A mathematical model that can be used for the analysis of thermal 
discharge from power plants into tidal estuaries and coastal waters is 
described. This transient, free-surface, three-dimensional model can 
be applied to predict the water temperature as a function of time and 
position in a specified region. 

In situations of practical relevance, the specific coastal or off- 
shore region will be a water body of irregular bottom topography with 
possible islands or keys. The user specified the boundary and boundary 
conditions, as well as the water depth distribution. Semi-diurnal tide is 
considered in the model. Hourly weather data is needed for wind stress 
calculation and heat exchange between water and the atmosphere. The 
ambient temperature is assumed of a sinusoidal form of 24 hour period. 
The ambient turbulence is included by an eddy viscosity and ^'^fusivity 
formulation. The appropriate values are to be calibrated against mea- 
sured currents. 
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SECTION 1 
INTRODUCTION 



The analysis of thermal discharges is important in order to minimize 
the environmental impact and to manage efficiently and safely the waste 
heat problems* The study of technological solutions to the problems of 
heated water disposal involves complicated relationships, such as the 
location, geometry end types of the discharge outlet, the flow condition 
and temperature of the receiving water body, the meteorological conditions 
of the site, the waste heat output of the power plant, etc. A thorough 
understanding of the thermal effects and physical processes of heated 
water dispersion to the environment is an essential part of serving the 
rapidly growing demand for electrical energy, while reducing the possible 
impact on the receiving ecosystem. 

The thermal effluent from a power plant will have variable conse- 
quences on the aquatic life of a receiving water body and the adjacent 
environment depending on the temperature rise. Therefore, the prime 
objective of the heated water discharge system is to bring the discharged 
water into thermal equilibrium with the surrounding water by bringing 
thermal outfall to the mainstream of the water body, whereby the mixing 
and convective processes will increase the surface heat transfer to the 
atmosphere. Thus, the temperature rise within the tolerance of natura! 
environmental conditions is very important on the disposal system de- 
sign and the standards for regulating thermal effluent. 

Under limited circumstances, in-situ measurements can serve for 
diagnostic and monitoring purposes for meeting the need of analyzing 
thermal impact on receiving water body. However, to provide a priori 
information about the nature and extent of thermal impact for site selec- 
tion and discharge system design, numerical modeling for simulating hy- 
drothermal behavior of the water body is imperative. 

During the past years, many numerical models have been developed 
for hydrothermal studies. Dunn et al. (1 975) gave a presentation of 
different models developed up to then. They applied those models to 
the Point Beach Nuclear Power Plant and compared the performance of 
various models in predicting a standard data base. A general conclusion 
that can be made from their analysis is that though some models may 
perform well under certain conditions, a generalized model which accounts 
for wind, current, tide, bottom topography and diverse meteorological 
conditions is yet to be developed. 

1 



Since 1974 the Thermal Pollution Croup at the University of Miami 
has endeavored to develop a mathematical package for hydrothermal studies . 
The primary motivation behind the effort was to develop a series of models 
which make minimal site restrictive assumptions enabling application to 
diverse basin and discharge configurations. Two separate formulations 
were made, one with the rigid-lid approximation and the other with the 
free-surface included. The details of the package and formulation are 
presented in a number of reports by Lee et al. (1978). 

The present report concerns the UM's free-surface model and its 
application to the Anclote Anchorage in Florida for waste heat discharge 
from a power plant. The features of the model are: a) three-dimensional, 
b) nonlinear, c) time-dependent, d) irregular topography, e* driving 
forces including wind, tide, heat and mass flux, f) graphical represen- 
tation of results of velocity and temperature fields, g) prediction of 
temporal and spatial variation of water surface. 

The descriptions of the main program and its subroutines, main 
algorithm and flow chart, program symbols, input data and logic para- 
meters, as well as the description of associated plotting programs are 
contained herein for the ready access of the computer program package 
to the user. A preliminary review on existing three-dimensional free- 
surface models, basic concepts of the present model, assumptions, ap- 
proximations, governing equations, initial and boundary conditions, 
finite difference implementation and numerical solution methods is pre- 
sented in Lee et at. (1979) and Carter (1977). 

The report also contains a plotting program which is used to analyze 
the results of the main calculation. A subroutine to compare the calcu- 
lated temperature field with that obtained by IR scanning is presented in 
the program. Note that the IR temperature is interpreted by hand from 
the mosaic film and then read in for isotherm plotting and comparison. 

The model has been tested for its adaptability. That is, the model 
allows for program modifications so that different initial and boundary 
conditions could be considered. The main program has several flag 
statements which make different usage of same program possible. Any 
program modification for the purpose of model transition should be made 
with care, and ih^t new program should be validated by sample runs to 
assure that the t feet of the modification is as desired* The same is 
also applicable to the plotting program. 

The program therefore contains two parts. Part 1 is primary and 
performs calculations; Part 2 is secondary and is for analyzing and plot- 
ting the results of Part 1. ANCMN is the driving program of Part 1. 
The input contains parameters, geometrical and initial data, or tape which 
stores intermediate results; output contains printout of results at preset 
time intervals, in both hard copy and tape form. The hard copy printout 
provides the base for analysis, upon which decisions and choice of the 
plots needed for further analysis and detailed comparison with measured 



results can be made* PLOTMN is the driving program of Part 2. The 
tape output from Part 1 is the main input. In addition, control cards 
assigning choice of plot, plot size, simulation hour and measured result 
for comparison are also required. Output is in printout and plot tape. 
The latter is used for plotting by CALCOMP plotter. 



SECTION 2 
RECOMMENDATIONS 



The model can be enlarged to handle any passive constituent, dis- 
solved or suspended, possessing arbitrary decay characteristics. The 
formulation of the constituent transport is based on the convection- 
diffusion equation, which is analogous to the thermal transport equation 
in the present model. The enlarged model would then be an ideal tool 
to study the ecological response of aquatic biota to the thermal effluents 
of the power plant. 

Th*. f >del can be modified to include bouyancy effects caused by 
fresh water /salt water sensing by including a salinity-dispersion equa- 
tion. This equation will be of similar form to the energy equation. 

The code is written for a constant grid size. Modifications can be 
made to incorporate a coarse grid for the complete field in comparison 
with a fine grid near the discharge location. This will allow a more 
accurate prediction of plume behavior in the near field A penalty in 
computational cost will be incurred. 



SECTION 3 
PROGRAM DESCRIPTION AND FLOW CHART FOR MAIN PROGRAM ( ANCMN ) 

DESCRIPTION OF PROGRAM ALGORITHM 

The governing partial differential equations are given in Table 1; 
the symbols and definitions are referenced in Nomenclature. The problem 
is set up as an initial-boundary value problem, so values of dependent 
variables are assumed known initially and prescribed on boundaries- 
Values at successive time steps are obtained by using a true explicit 
scheme. The leap-frog finite difference format is used to calculate sur- 
face elevation n and two horizontal velocity components, u and v, at 
time n+1, where n is the present time step. The variables at limes n-1 
and n are all known* The sequence in which calculations are performed 
is as follows: 



1. Integrate the surface elevation equation using central-time central- 
space (CTCS) differencing. That rs, n°" , n° and u n , v n and h. 
Note that h is independent of time while the present water depth H 
= h+n n js needed in tfjls calculation. In the subroutine BETA, not 
only n n but also n (i, j, k) , the modified vertical velocity in 
transform (a) plane, is accomplished immediately after the n compu- 
tation. 

2. The next task in the sequence is to calculate the nonlinear inertia 
terms that appear in the horizontal momentum equations. Here, two 
subroutines are involved: BNRT1A is 

TABLE 1. Governing Equations 

Continuity Equation*: 



3Hu , 3Hv , u 3 n , 3 n _ A 



u Momentum Equation: 



|Hu + |Huu + 3Hvu +H 3QU + 3u3 r, _ fyH 

3t 3x 3 y 3a 3 a 3t 



p3x 3 8x 3x' H 



Av 9 : u 



v Momentum Equation: 



3 t 3 x 3y 3c 3a 3t 



- H 3P nil f.-, 3H 1 3n 1 t AV 3? Y- 



p 3y 



Energy Equation: 



3HT ^ 3HuT . 3HvT . u 3 S?v n . 3T Jrrj 

3— + TIT" + 3"Y~ " ~ + (1 + 0) T? It 

" PT T7^ + B h l Fx (H IT + ly (H T7 )] 



* This equation is vertically integrated to yield a) prognostic equation 
for y, and b) synoptic equation for 8; they are 

at -r l 3 x 3 y ' ao 



« - a 3 i a. 1 ro r 3Hu 3Hv, H ^ 

a -H3i + H/ ( "rr + it 1 da * 



a y 

The latter, upon transformation, yeilds the actual vertical velocity 



for interior points, while ABNR3 is for open boundary points. Note, 
fcr the Anclote Anchorage sample problem, the open boundaries are 
at j = 1 and j = 14, and the imposing tides are applied at points 
immediately outside these open boundaries. Therefore, program 
modif!-ation is needed if different open boundary conditions are 
employed. 

3. Following the inertia terms computations, which may be skipped if 
the Rossby Number is very close to zero, new values of u and v 
at n+1 are computed for all points in the grid. Again, the leap- 
frog and central-space schome is used, but DuFort-Frankel differ- 
encing is applied to the vertical momentum diffusion terms* Two 
subroutines are called here, BVEL for interior points and ASAF3 
for open boundary points. Since velocity at all points is calculated 
without distinction, a subroutine CIVENU is needed to specify the 
given discharge or flowrate at particular points, that is, to replace 



calculated velocities at those points with known values. 

Steps 1, 2 and 3 are calculations for surface elevation, modified 
vertical velocity and horizontal velocity components; they constitute the 
V-calculation. Whether the V-calculation is to be carried out or not 
depends on flagged statement KVEL = 1 or 0. The next group of calcu- 
lations is for thermal transport, or T-calculation, and it involves the 
energy equation only. Similarly, this group is to be flagged by state- 
ment whether KTEMP is unity or zero. 

4. The convective term in the energy equation is calculated next, using 
either a given velocity field, in the case of KVEL = 0, or the pre- 
sently known velocity field at n+1, in the case of KVEL = 1. The 
subroutine for this purpose is CONV. 

5. The energy equation itself is then integrated over time to obtain T 
at n+1. The forward-time, central-space (FTCS) and DuFort-Frankel 
differencing for the vertical diffusion term are used in the calcula- 
tion. The subroutine involved is TCOMPT. Since temperatures at 
all points are computed without distinction as to whether the points 
are with given temperatures or not, subroutine CIVENT is needed to 
respecify the temperature at the given points. 

Clearly, Steps 4 and 5 make up the T-calculation. In either of V- 
or T-calculations, vertical velocity w^is involved. Instead of w, the 
rate of change of surface elevation, — , is used for convection in ver- 
tide direction. 

6. The actual vertical velocity, w, is computed when it is^rteeded foe 
printout. The subroutine is WCAL, and instead of u n and v 
which are defined at half J and half I respectively, the interpreted 
velocities at center of grid cell are used in this calculation. Since 

a space-staggered scheme is used, the water level and vertical velo- 
city are described at the center of grid cell, while the horizontal 
velocities are described at the edges of cells. 

7. The real time (or simulation time) is checked and Steps 1 through 6 
are repeated; that is, the above procedure is repeated for n+2 using 
values at n+1 and n. 

Reference to the flow chart presented in Figure 1 will clarify the 
description of program algorithm. 

FLOW CHART 

Figure 1 shows the main flow chart of the three-dimensional, free- 
surface program applied to the Anclote Anchorage. In the flow chart, 
the subroutines and their functions are described briefly. Table 2 lists 
the subroutines called in the main-program ANCMN. 
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Figure 1. Flow chart for calculating program ANCMN and its subroutines 



TABLE 2. Subroutines Required in Main Calculating Program ANCMN 



No. 



Name 



BAYBOT 



BAYINl 



10 



11 



READT 
IRREAD 

EQTEMP 
BETA 

BNRTIA 

ABNR3 

BVELS 
ASAF3 

GIVENU 



Description 



Reads in bay bottom sounding 
(H) and various grid book- 
keeping matrices, calculates 
depth matrices accordingly. 

Initializes the dependent 
variables such as U, V, W, T, 
ETA, RX, RY also set 1CT to 
values outside the calculation 
region* 

Reads in stored data from 
tape for continuing run. 

If required, it reads in IR 
data as initial temperature; 
thus it replaces T input 
from READT . 

Calculates the equilibrium 
temperature over that hour. 

Calculates surface elevation 
ETA and vertical velocity 
OM (a). 

Calculates inertia terms of 
momentum eq. for boundary 

pts. 

Calculates inertia terms of 
momentum eq. for boundary 
pts. on north and south 
exit of the anchorage. 

Calculates velocity at in- 
terior pts. - main part of 
velocity calculation. 

Calculates velocity at boun- 
dary pts. along north and 
south boundaries. 

Specifies velocity at control 
pts., such as discharge, in- 
take and river head. 



Remark 



Skip if LN=1, depends 
on problem and grid 
work, data file AMATN 
preferred. 

Also depends on probelm 
and run conditions. RX 
RY are inertia terms. 



Skip if LN=1< 



No need if T-calculation 
began with ambient 
temperature. 



Number of hour is NCY. 



Mandatory in 
V-calculation. 



Skip if ROSSBY = 0. 



Skip if ROSSBY = 0. 



TABLE 2. Subroutines Required in Main Calculating Program ANCMN 

(Continued) 



No. 



Name 



Description 



Remark 



12 

13 

14 

15 

16 



CONV 

TCOMPT 
CIVENT 

WCAL 
ANCPR 



17 



18 



TPRLOK 



STORET 



19 



20 



21 



ZZ1 



AMATN 



C2007 



Calculates convective terms of 
energy equilibrium. 

Calculates temperature. 

Specifies temperature at dis- 
charge outlet. 

Converts the OM (n) vertical 
velocity to physical vertical 
velocity W. 

Prints surface elevation, 
velocity (mag. with dir.) and 
temperature at 4 chosen loca- 
tions after completion of 
marching of DT, i.e. printout 
after each time step. 

Prints out the velocity field, 
water depth, surface elevation 
and temperature field after 
one hour of simulation time. 

After TPRLOK for printout, 
the same data and relevant 
parameters are stored onto 
tape, either for later plotting 
or for continued calculation. 

A subroutine called by ANCPR 
and TPRLOK; it is for finding 
velocity dir. based on U and 
V components. 

Is a data file containing all 
grid matrices and bottom depth 
matrix to be read in BAY80T. 

Climatical data used in *he 
run; the data is on hourly 
basis. 



m i nr i ii i ' iii i saacrag 



Depends on run condi- 
tion. 



Locations are to be 
chosen by user. 



Skip if KSTORE = 0. 



Same data could be in 
card form, inserted 
after TINIT card. 

Could also b,e in card 
form . 
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SUBROUTINE DESCRIPTIONS 

This section describes the subroutines used in ANCMN, in order of 
their appearance* 

BAYBOT 

It reads the marker matrix MAR; the elevation matrix ELEV which is 
changed to depth matrix H by adding a constant STAGE; then four more 
marker matrices, MEX, MEY, MX and MY. It interpolates H according 
to marker MEX, MX and MY to find additional depth matrices, called HB, 
HU and HV, respectively* 

BAYINI 

It initializes most of the variable matrices. Since the initial condi- 
tion for velocity is a quiescent condition, U1, U2, U3, VI, V2, V3, RX, 
RY, ETA1, ETA, ETA3, UB, VB, OM, W and TC are set to zero. The 
quiet bay is assumed to have a constant temperature TINIT to begin with, 
so Tl, T2, T3 and TB are set to TINIT. 

READT 

This subroutine reads in input parameters, physical quantities 
and intermediate results stored on tape. 

IRREAD 

It reads in in-situ measured or IR scanned temperature as an alter- 
native initial condition to temperature calculation. This temperature is 
interpolated by hand and stored in the form of matrix Tl. 

EQTEMP 

it calculates the equilibrium temperature T and surface heat ex- 
change coefficient K of a natural water surface 6 . The procedures for 
these calculation j arl as follows: 

1. T . = T - (14.55 + 0.114 T ) (1-f) - [(2. 5 + 0.007 TJ(l-f)] 3 

o a a a 

where T . = dewpoint temperature in °F 
T ~ air temperature in °F 
f a - relative humidity in fraction of unit 

2. S = 0.255 - 0.0085 T ave + 0.000204 T ave 2 

where T = j(T + T .) , and s is an intermediate step 

T* ve = amb?ent surface temperature in °F 

3. f(u) = 70 + 0. 7 u 2 , and u is wind speed in mph 

11 



4. K = 15.7 + (6 + 0.26) f(u) 

where K = surface heat exchange coefficient in BTU/(ft 2 day °F) 
$ H 

5 ' T e= T d + r 

where T = equilibrium temperature in °F 

H e = gross solar radiation in BTU/(ft 2 day) 

Note: The T , f, u, H and T are climatological data; however, care 
must be taker? for the hourly d&a TAIR, HUMID, WIND, SRAD and TSURF 
are in metric units* Therefore, in the above calculations, the basic data 
must first be transformed into English unit, then the final results, T and 
K , must again be transformed back to metric units. T (TEQ) and K 
(^K) are used in subroutine CONV for T-calculation. e 

BETA 

Computes n n+ (ETA3) and Q n (OM) by using central differencing 
from the continuity equation. The vertical integration is done using 
Simpson's rule. The following symbols are used: 

DR = depth at half-integer j point on right edge, i + 1, of the cell. 

DL = depth at half-integer j point on left edge, i, of the cell. 

D2 = depth at half-integer i point on upper edge, j + 1 , of the cell. 
D1 = depth at half-integer i point on lower edge, \, of the cell. 
DHUX = ^ 

dhvy = ?hy 

AH = total depth at the center of the cell, a half-grid point. 

BNRT1A 

It computes the sums of the nonlinear inertia terms, RX and RY, 
in the x and y momentum equation at each interior point of the domain. 
Note that RX (i, j, k) = RY (i, ]. k) = for k = KN. The following 
symbols are used, in their order of appearance: 

AH = depth either at u-point or at v-point 

DET = — 

OHUUX =|iiJJH DHUVX = 1HH1 

o X o X 
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DHUVY = ~^ DHVVY = lH — v 



3 y 3 y 

D2 = depth at forward half-grid point in either x or y direction. 

D1 = depth at backward half-grid point in either x or y direction. 

UBAR2 = average u at forward half-grid point. 

UBAR1 = average u at backward half-grid point. 

VBAR2 = average v at forward half-grid point. 

VBAR1 = average v at backward half-grid point. 

E2 = averaged n at (i, j + 1) or (i + 1, j) 

El = averaged n at (i, j) 

D2 =5 depth at (i, j + 1) or (i + 1 , j) 

D1 = depth at (i, j) 

DUOMS = 1^2- DVOMS = |*5. 

da o a 

DUS « |H DVS = !* 

3a oa 

ABNR3 

It computes RX and RY for points on boundary. The tide heights 
at half-grid points outside the north and south boundary are computed 
first. Since the open boundary is in x-direction, that is, only the v- 
point appears/ RX = and RY is given by 

= |Huv + |Hw 3u£. 3v _|n 

3x 3y 3a 3a3t 

In the computation, it is assumed: 

02 = depth at half-grid point just outside the boundary 
= HV (i, JN) + North Tide Height, if D2 is HB (i, JN) 
= HV (i, 1) + South Tide Height, if D2 is HB (i, 0) 

VBAR2 = v at half-grid point just outside the boundary 
= V2 (i, JN, K) if it is north boundary 
= V2 (i, 1, k) if it is south boundary 
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BVELS 

It calculates u (U3) and v (V3) at interior points by central 
time differencing and DuFort-Frankel scheme. Note *.hat the horizontal 
diffusion of momentum is neglected in the model. 

The equations for u and v are 

RTF _fv 9 F + j3T H RX 

H TT ' fU 9 WZ i TT H RY 

where RX and RY are the nonlinear inertia terms. 

The following symbols are introduced for briefness. 

A2 = Coriolis term 

A4 = g |i or g |- n 
3 3x y 3y 

A5 = J RX or J RY 

A6 = the rest of vertical diffusion term. 

ASAF3 

This calculates v n (V 3) at the v-points on the south (j = 1) 
and north (j =JN) boundaries. It is similar to ABNR3. The tide heights 
at imaginary half-grid points ju$t outside the south and north boundary 
are computed first. The term — is calculated at v-point at both j = 1 
and JN, with tide height at outs&e half-grid point and n at inside half- 
grid point. Since the u velocity on these boundaries are assumed zero, 
the Coriolis force term B2 is set to zero- Symbols B4, B5 and B6 stand 
for pressure, convection and diffusion terms respectively. 

CIVENU 

It specifies velocities at cooling system outlet and intake. These 
velocities are determimjd from power plant flowrate. The river flowrate 
is also simulated by irrposing velocities at river entry point* 

CONV 

It computes the sum of the convective terms in the energy equation 
at each point. Note that T is designated at half-grid points/ The follow- 
ing symbols are used in this subroutine: 

AH - depth at point (i, j) where TC (i, j, k) is to be calculated 

H 



OR = depth at forward u-potnt 
DL = depth at backward u-point 
D2 = depth at forward v-point 

Dl = depth at backward v-point 

3T 
DTZ (i, j) = temperature slope t- at the surface 

a a 
UR = forward u but at T-level 

UL = backward u but at T-Ievel 

DHUTX = 1^1 

VR = forward v but at T-level 

VL = backward v but at T-level 

DHVTY = —^ 

3 y 

The convection term TC is written for 

Tr _ 3HuT ^ 3HvT ^ u 3 flT ^ ri x . 3T 3 n 

TC * FT" + 7y~~ + H 77" + {1 + a) To ft 



TCOMPT 

This subroutine computes temperature T (T3) for each point \r 

the domain. The equation is 

H3t " „2 T^ 4 * B h ( T^ + T^ ) H TC 

where TC is the convection term obtained from subroutine CONV. The 
boundary condiiion on solid boundary is adiabatic, but on open boundary, 
T is^sjumed known. There are three different formulas for computing 
the i-^j- although the same DuFort-Frankel format is used throughout. 
This ft because the surface temperature slope is DTZ (i, j) and the bottom 
is adiabatic. The symbols here are 

D2TX = -2-L- DZTY = ^-~ 

3x^ Dy^ 

CIVENT 

It specifies the temperature of cooling system discharge water and 
river delivery. 
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WCA*. 

This subroutine calculates vertical velocity w (i, j, k) from modified 
vertical velocity a (i, j, k) by using o- trans formation formula. The fol- 
lowing symbols are used: 

DEX , |i DEY . |i 

DAHX = |H DAHY = |£ 

ANCPR 

It prints continuous records of elevation, velocity and temperature 
at certain particular points. 

All data are printed or a single line, with the first item on the 
line being the total simulatbn time TTOT in sec. Items for each point 
are: surface elevation, ref'iltant velocity, direction in which the velocity 
vector points (deg positive clockwise from North), and temperature, in 
that order* 

TPRLOK 

This subroutine performs the major printing tasks. The following 
variables are printed hourly, or controlled by printout interval TPRT: 

AVR = resultant velocity, cm /sec 

ANC = direction in which the velocity points 

W = vertical velocity, cm /sec 

ETA = surface elevation, cm 

T2 - temperature, deg C 

After TPRLOK is executed, the main program ANCMN prints out the total 
depth, H = h + n, at each point. Thus, the hourly printout of the rele- 
vant variables is completed. In this subroutine, there are two flags, 
KUV and KPROF. If both are zero, it skips the printing of velocity 
components u and v either presented in layers fk - 1, KZ) or in cross 
section along x-direction (j = 1, JM). These two flags are flipped only 
by interchanging the statements in the subroutine. 

STORET 

It records all the relevant data and results of the preceeding simula- 
tion hour. STORET puts one EOF on tape after each block, while the 
nuin program ANCMN puts another EOF after the last block of data is 
recorded . 
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ZZ1 

This subroutine finds the direction of the resultant horizontal velo- 
city. 
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SECTION 4 
LIST OF PROGRAM SYMBOLS OF MAIN PROGRAM 



This section presents the program symbols and their definition in 
alphabetical order. In many cases, the symbols are described with the 
aid of diagram to show the definition* 

DESCRIPTION OF MAIN VARIABLES 

The relative position and designation of variables are shown in 
Figure 2* The water depth, h, is described in integer values of i and 
j; the u-component is described at half-integer value of j and integer 
values of i and k; the v-component at half-integer value of i and integer 
values of j and k; the w-component at integer value of k and half- 
integer values of i and j; the surface elevation, n, is described at half- 
integer values of i and j; and the temperature, T, is described at half- 
integer values of i, j and k. The modified vertical velocity, 2, is de- 
scribed at the same place as w-component* Figure 3 shows the space- 
staggered grid system in horizontal projection* 

Table 3 lists all the symbols used for dependent variables appearing 
in the program* Since three levels of time step are used, same variable 
at different level is assigned with different symbols. The rule for 
symbolizing the dependent variables is: for variable F(i^x, j£y, kl a, 
nM) , Fl(i, j, k) is used to denote the value of variable at n-1; F2(i, 
j, k) is the present value; while F3(i, j, k) is the value at n+t thus to 
be computed; and FB(i, j, k) is the interpreted value of F2(i, j, k) at 
a set of grid points differing from where it is designated. 
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v(i,j+l,l) 



u(i,j,l) 




u(i+l, j,l) 



u(i+l, j,2) 



h(i,j) 



Figure 2. Relative position and designation of the variables 
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y(j) 



all points on 
these columns 
have same 
subscripts 



j+ l - + _ | - 




all points in this 
+— square have same 
subscript i,j 

^C ] all points on 
V these rows have 
) same subscript j 



water level (T) 
z veloc: 



- -A- - depth (h) 



evel nn 
ity (w) J 

h) J 



x velocity (u) 
y velocity (v) 



x(i) 



Figure 3. Space-stagaered arid system - plan 
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TABLE 3. Symbols used in t 


:he program 


Symbols 


Argument 


Description 


U1, U2, U3, UB 


(i, j, k) 


u-component 


VI, VI, V3, VB 


(i, j, k) 


v-component 


T1, T2, T3, TB 


U, j, k) 


Temperature T 


ETA1, ETA2, ETA3 


U, j) 


Surface Elevation n 


TIDE1N, TIDE2N, TIDE3N 


(i) 


Tide Outside North 
Boundary 


T1DE1S, TIDE2S, TIDE3S 


(i) 


Tide Outside South 
Boundary 


DTZ 


(i, j) 


Heat Exchange at Water 
Surface 


RX; RY; TC 


(i, j, k) 


Convection Terms in u; 
v; T equation 


W; OM 


(i, j, k) 


w-component ; Non-dimen- 
sional ft 


WAT 


(i, j) 


Total Water Depth or 

h + n 



Note: In the program ANCMN, ETA2 is labeled as ETA and ETX, the 
former is used for calculation while the latter for printout* 

MARKER MATRICES 

The following integer-valued matrices are introduced to describe the 
grid system and to distinguish boundary from interior. 

MAR (i, j) 

MAR (i, j) identifies nodes in the full-grid system, i.e. 

MAR =0, (i, j) outside of boundary, hence no calculation 

MAR = 1, (i, j) inside or on a boundary, as shown in Figure 4. 
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Figure 4. MAR (i, j) matrix 



MEX (i, j) 



MEX (i, j) provides marker to the half-grid system with reference 
to the y-direction boundaries, i.e. 

MEX =0, (i, j) outside of y-boundary, or exterior 

MEX = 1, (i, j) just inside an east boundary 

MEX = 2, (i, j) just inside a west boundary 

MEX = 3, (i, j) nowhere near to y-boundary , or interior, as shown 
in Figure 5. 
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Figure 5. MEX (i, j) matrix 
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MEY (i, j) 

MEY (i, j) provides marker to the half-grid system with reference 
to the x-direction boundaries, as shown in Figure 6. 

MEY =0, (i, j) outside of x-boundary, or exterior 

MEY = 1, (i, j) just inside a south boundary 

MEY a 2, (i, j) just inside a north boundary 

MEY = 3, (k, j) nowhere near to x-boundary, or interior 



r < (& { < <s<. 



I 






/r> > ?> r > >?> 



I 

i 



' {{jtttu 



^u. 



? t r ;", » i rt* *> * * f f j *} » ' 



Figure 6. MEY (i, j) matrix 
MX (i, j) 

MX (i, j) provides marker to u-points, as shown in Figure 6. 

MX =0, (i, j) outside of y-boundary, or exterior 

MX = 1, (L j) on east boundary 

MX = 2, (i, j) on west boundary 

MX - 3, (i, j) nowhere near to y-boundary 
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MY (i, j) 

MY (1, j) provides marker to v-points, as also shown in Figure 6* 

MY =0, (i, j) outside of x-boundary, or exterior 

MY = 1, (i, j) on south boundary 

MY = 2, (i, j) on north boundary 

MY = 3, (i, j) nowhere near to x-boundary 



. MY I 



M -. 



,,- .vridVfy 



*)}?l n>> 



H— 



--1 —3 — 3 — 3 ~3 —3 



—3 



,, ? , ? } } s> > s i > > j / > * />, < >>?>> 



f-2 



-7-0 — 1 — 3 MX 



/r/>\;//> 






' f J >>f/ 



Figure 7. MX (i, j> and MY (i, j) matrices 
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DEPTH MATRIX AND ITS DERIVATIVE 

The bathymetry of the area of imprest Is given by the matrix ELEV 
(i, j) designated at full-grid points. The values in feet are positive for 
take or inland waters as elevation above MSL. However, for coastal 
water, the depths are read from the survey chart and are designated by 
positive values. For certain periods of the year, the water level may 
differ from MSL; a stage (STACE1) is added to the ELEV (i, j) to obtain 
the actual depth matrix H (I, j). Note that the values of H (I, j) are In 
cm. To facilitate the calculation, the matrices HB (i, j) , HU (i, j) and 
HV (i, j) are derived from H (i, j), and are for depths at half-grid 
points, u-points and v -points, respectively. Thus, the following real- 
valued matrices are introduced. 

ELEV (i, j) elevations of the bottom with respect to MSL 

HB (i, j) depth at half-grid points, in accordance with MEX/MEY 

HU (i, j) depth at u-points, in accordance with MX 

HV (i, j) depth at v-points, in accordance with MY 

SIZE OF THE MATRICES, OR DIMENSIONS OF SUBSCRIPTED QUANTITIES 

Let the grid work consist of IN x JN x KN nodes, i.e. IN nodes in 
x(i) direction, JN nodes in y(j) direction, and KN levels in a(k) direction. 
Then there are IM = IN - 1 half-grid points in x-direction, JM = JN - 1 
half-grid points in y-direction, and KZ = KN - 1 layers in o- direction. 

The values IN, JN, KN, IM and JM are the parameters to be specified 
at the beginning of the program and are determined by the grid used. 
Therefore, the dimensions of the matrices are given in terms of these 
parameters. Table 4 shows the size of the matrices already defined. 

TABLE 4. Size of the Matrices 



Symbol 


Least Size 


Given Size 


Ul, U2, U3 


(IN, JM, KN) 


(IN, JN, KN) 


V1, V2, V3 


(IM, JN, KN) 


(IN, JN, KN) 


T1, T2, T3, TC, TB 


(IM, JM, KZ) 


(IM, JM, KN) 


UB, VB 


(IM, JM, KZ) 


(IM, JM, KN) 


W, OM 


(IM, JM, KN) 


(IM, JM, KN) 



25 



TABLE 4. Size of the Matrices 
(Continued) 



Symbol 


Least Size 


Given Size 


RX, RY 




(IM, JM, KZ) 


(IN, JM, KN) 


ETA1, ETA 


. ETX, ETA3 


(IM, JM) 


(IM, JM) 


WAT, DTZ 




(IM, JM) 


(IM, JM) 


ELEV, H 




(IN, JN) 


(IN, JN) 


HB 




(IM, JM) 


(IM, JM) 


HU 




(IN, JM) 


(IN, JN) 


HV 




(IM, JN) 


(IN, JN) 


TIDE1N ... 


TIDE3S 


(IM) 


(IM) 


AVR, ANC 


• 


(JM) 


(JM) 


MAR 




(IN, JN) 


(IN, JN) 


MEX, MEY 




(IM, JM) 


(IM, JM) 


MX 




(IN, JM) 


(IN, JN) 


MY 




(IM, JN) 


(IN, JN) 



AVR and ANC are used to facilitate printout of velocity and angle re- 
spectively. 

OTHER SYMBOLS OCCURRING IN PROGRAM ANCMN 

ALREF: Reference horizontal length L in cm* 

AV: Vertical eddy viscosity, estimated by means of 

A y = 0,0018 H 4/3 , H in cm, A y in cm 2 sec~ 1 . 
BH: Horizontal *ddy diffusivity, estimated by means of 

B h = 0.0018 L 4/3 , L in cm, B h in cm* sec" 1 . 
BV: Vertical eddy diffusivity, estimated by the same formula as A , 
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thus turbulent Prandtl No. is one. 
DHR: Time increment in hour as simulation continues. 
OTX: As a check on when to printout. 
DS: Increment in o-directton, in fraction of unit. 
DT: Time step in second. 
DX: Increment in x-oirection, in cm. 
DY: Increment In y-dlrection, in cm. 
OUMS: 2DS- 
DUMX: 2DX. 
DUMY: 2DY. 
EST: Eastern standard time in the day of simulation. 

FCOR: Coriolis factor = 2W sin(latitude) , sec" 1 . 

e 

W = earth's angular rate of rotation. 
e 

C: Earth's gravitation = 980 cm sec" . 

I: Index for x-axis. 

ICO: Flag. Set as 1 initially; it changes to when the calculation be- 
comes unstable. 

IM : Maximum number of half-grid point in x-direction. 

IN: Maximum number of full-grid point in x-direction. 

J: Index for y-axis. 

JCTR: Index for simulation hour. 

JM: Maximum number of half-grid point in y-direction. 

JN: Maximum number of full-grid point in x-direction. 

K: Index for a-axis. 

KN : Maximum number of full-grid point in c-direction. 

KZ: Maximum number of half-grid point in c-directicn. 
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KSTORE: Flag* Set as 1 to store hourly result on tape* 

Set as if no store is needed* 

KVEL: Flag* Set as 1 if velocities wb to be calculated, otherwise set 
as 0* 

KTEMP: Flag. Set as 1 if temperature is to be calculated, otherwise 
set as 0. 

th 
LN: Set as 1 for 1st run of present case; set as n for subsequent n 

run* 

MBLOK: Data block number which is to compare with data block NBLOK 
which is to be read in. Used only when LN > 1. 

NBLOK: Index for data block* 

NCASE: Case number. 

NCY: Number of hours to be simulated in this run. 

QQ: 57.3, used for changing from deg to rad. 

ROSSBY: Rossby number* 

RR: Water density , = 1. 

RWEX; Number of hours between cllmatological input data* = 1. 

TABN : Ambient water temperature outside of north entrance* 

TABS : Ambient water temperature outside of south entrance* 

THETA: Angle between north and y-axis, clockwise positive* 

TINIT: Water temperature at initial instant before the waste heat dis- 
charge start* 

TPRT: Time between printouts, in sec. 

TTOT: Total simulation time, in sec* 

TZ: Record of time for hourly printout. 

TZERO: EST hour at the beginning of present simulation run* 

The following symbols are used to specify tidal condition* 

AMPL1T: Tide amplitude, in cm. 
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DPHASE: Phase laq per Ax, in hour* 

PERIOD; Tide period, in hour. 

PHASE: Phase difference between tides at north and south entrance* 

STAGE: Difference in cm between daily mean level and short-term 
(weekly) average sea level. 

STACE1: Difference in cm between short-term (weekly) average sea 
level and long term average level (MSL). 

TSHIFT: Time shift for adjusting tide with EST, in hour* 

The following symbols are used to specify the hourly climatoiogical 
conditions* 

TAIR: Ambient air temperature, deg C. 

HUMID: Relative humidity, fraction. 

WIND: Wind speed, cm sec* 1 . 

WDIR: Direction from which wind is coming, deg measured clockwise 
from North. 

SRAD: Cross solar radiation, in BTU/(ft 2 day). 

TSURF: Surface water temperature, deg C. 

The following symbols are related to climatoiogical data and appear 
in the calculation of wind stress, equilibrium ambient temperature and 
heat exchange at surface. 

EPSLON: Direction to which wind blows, in rad. 

WPR: Wind speed in m sec" 1 . 

CTEN: Empirical constant appears in wind stress formula. 

TAU: Wind stress t. 

TAUX: x-componrnt of wind stress x . 

TAUY: y-component of wind stress t . 

TDEW: Dewpoint temperature, deg C. 

TEQ: Equilibrium temperature, deg C* 

SK: Surface heat exchange coefficient in cal/(cm 2 sec °C)* 
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SECTION 5 
PREPARATION OF SIMULATION RUN 

This section describes the preparation work needed for ANCMN run* 
The flow chart and the associated subroutines in Figure 1 and Table 2 
are referred to In the following description. 

1. Specify number of full-grid points, IN, JN, KN and number of half- 
grid points, IM, JM, in PARAMETER statement. Although the domain 
of solution under consideration is usually smaller than the rectangular 
space of IN x JN x KN, the marker matrices will assure that the grid 
points outside of domain skip the calculation* To have a clear print- 
out, the variables at off domain point have been set to 10 . This 
value is beyond the capacity of the computer printout in printing real 
numbers (F format) so that stars will be printed and show the off 
domain area* 

2* Specify run number by input data LN, card #2: 

For LN = 1, i.e* first run, data file or card deck of AMATN 
is needed. 

For LN > 1, i.e* subsequent run, tape with previous result 
is needed. 

Specify flag for storage by KSTORE, card #3: 

For KSTORE = 0, desire no storage. 

For KSTORE = 1, tape must be provided for storing results* 
Specify flag for velocity calculation by KVEL, card #4: 

For KVEL = 0, no V~calculation, thus thermal dispersion only. 

For KVEL * 1, do V-calculation, thus circulation included* 

Specify flag for temperature calculation by KTEMP, card #5: 

For KTEMP = 0, no T -calculation, thus only a hydrodynamic 
model* 

For KTEMP = 1, do T-calculation, a complete hydrothermal model* 
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Specify data block number MBLOK, to make sure the data read in 
from tape is correct, card #6* 

Specify number of hours to be simulated in this run by NCY, card #7* 

Specify the time between successive printouts by TPRT, card #8. 

Specify grid size by input data DX, DY, DS, card #9. 

Specify time step OT and tide data ST ACE, AMPLIT, PHASE, DPHASE, 
PERIOD, TSHIFT by input data, card #10. 

Specify Coriolts factor PCOR and stage STACE1 by data, card #11. 

Specify the angle between North and the y-axis of grid system by 
THETA, card #12. 

Specify reference length ALREF and kossby No. ROSSBY, card #13. 

Specify number of hours between weather observations RWEX, card #14* 

Specify TZERO, the Eastern Standard Time when the simulation starts, 
card #15. 

Specify water density, vertical eddy viscosity, vertical eddy and 
horizontal eddy diffusivity, RR, AV, BV, BH, card #16. 

Specify initial temperature TIN IT, a constant for whole domain, card #17. 

3. In general, the first run of present case has: 

LN = 1, KSTORE = 1, KVEL = 1, KTEMP = 1, MBLOK = 0. 

Then the subroutines BAYBOT and BAYINI are used to initialize 
the calculation. This includes reading matrices, MAR, ELEV, MEX, 
MEY, MX, MY, by BAYBOT from data file AMATN. The same sub- 
routine calculates the derivative height matrices, HB, HU and HV. 
The initialization of various variable matrices is done in subroutine 
BAYINI and in the main program itself. 

4. In general, the continued n run has the same NCASE with: 

LN--n, KSTORE = 1, KVEL • 1, KTEMP = 1. 

MBLOK = index number of the data block which is to be read in; the 
calculation will continue thereafter* In fact, the data being read con- 
tains all the information needed to continue the run. However, to 
allow for the freedom of matching tide of different amplitude, period 
and phase shift, an additional card (#18) specifying NBlOK, TTOT, 
DT, EST, AMPLIT, PHASE, DPHASE, PERIOD and TSHIFT is needed. 
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The data may ba sama as those contained in tha tapa or dl ffarant 
from tham so that tha calculation goas on to follow anothar tida 
format, in accordance with this changa of tide, tha NBLOK, TTOT, 
DT, EST may ba reset* 

5* Tha main loop in the main program ANCMN is the hourly simulation 
loop, which is started with hourly climatological data card containing 
TAIR, HUMID, WIND, WDIR, SRAD and TSURF* The wind stress 
and equilibrium temperature are then computed and held thereafter 
as constants throughout that hour* 

Tha main part of the hourly loop is an internal loop for it Increment, 
in which the main calculation is done in the order of (n, fl), (u, v) , 
T, W, then a printout of elevation, surface velocity and surface 
temperature at certain chosen half-grid points* 

6. The V-calculation controlled by flag KVEL consists of subroutines 

BETA, BNRTIA, ABNR3, BVEL, ASAF3 and CIVENU. BETA computes 
n and fl. BNRTIA and ABNR3 compute the convection terms, RX and 
RY, for the momentum equations; this computation is decided by 
whether ROSSBY is zero or not* The (u, v) calculations are done by 
BVELS and ASAF3. The given velocities at contro! points ^r% re- 
specified by GIVENU. The BNRTIA and BVELS are for interior 
points while ABNR3 and ASAF3 perform tne same purpose except for 
normal velocity points along open boundaries, where the water eleva- 
tion is specified as a function of time* 

The T-calculation controlled by flag KTEMP consists of CONV, TCOMPT 
and CI VENT* CONV computes the convective term TC, then TCOMPT 
computes T, and CI VENT respecines T at discharge points. 

After the completion of marching forward to (n+1).it, the variables are 
relabeled and UB and VB are computed as the horizontal components 
of velocity at centers of (I, J) blocks* Finally, before the printout 
of newly obtained variables at fixed locations to serve as flow develop- 
ment at fixed point, the surface velocity and temperature at a critical 
point are compared with preset values to see whether an instability 
has developed* If instability does occur, the program terminates after 
producing a hard copy of the latest result. 

7* The subroutine ANCPR produces step-by-step records of surface 
elevation, surface velocity and surface temperature at certain chosen 
points* These points are selected because of the variables that are 
believed to undergo the most change, as they are close to open 
boundaries, river exit and discharge outlet. 

The hourly loop included subroutines SOTRET and TPRLOK too; the 
former stores the hourly results as well as all pertinent data onto 
tape for later uses, and the latter produces a printout of resultant 
horizontal velocities, vertical velocities, temperatures at four levels 
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and elevation of free surface. In addition, the main program ANCMN 
itself does the calculation and printout of total water depth, H = h + 
r\, before starting next hourly loop. 

8. ANCMN performs the hourly loop NCY a number of times. 

Therefore, NCY number of climatologicai data cards are needed to 
provide the necessary data. 
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SECTION 6 
INPUT DATA 

The input cards for running ANCMN are given in Table 5 below. 
Note that the data symbols have already been defined in the previous 
section; however, the following remarks should be considered* 

* Free format is used for all data input. 

* Distinction must be made for integer and real number. 

* The order of these cards must be followed* 

TABLE 5. Input Data for ANCMN 



Input 


Card 
Content 


Symbol 


Definition /Value 


#1 


1 


NCASE 


= Case No. 


#2 ; 


1 


LN 


= 1, if it is first run then data file 
or card deck AMATN is needed in 

#18 th 
= n, if it is n* n run 


#3 i 


1 


KSTORE 


= no store, then there are no 

continued runs 
= 1 store intermediate results on tape 

for plotting or next run 


#4 


1 


KVEL 


= no V-calculation, i.e. dispersion 
of T by given (u, v) field 

= 1 do V-calculation , so momentum 
is under dispersion 


#5 


1 


KTEMP 


= no T-calculation, i.e. hydrody- 

namic only 
= 1 i.e. hydrothermal model 


#6 


1 


MBLOK 


= No. of latest hour of last run 
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TABLE 5. Input Data for ANCMN 
(Continued) 



r 



Input 


Card 
Content 


Symbol 


Definition /Value 


#7 


1 


NCY 


= Number of hours intended for 
simulation in this run 


#8 


1 


TPRT 


= 3600 seconds, hourly loop 


#9 


3 


DX 


= x-direction grid size in cm 






DY 


= y-direction grid size in cm 






DS 


= ^direction grid size in nondimen- 
sional unit 


#10 


7 


DT 


J ^max 






STAGE 


= Average level of tide-MWL, in cm 






AMPLIT 


= Amplitude of tide, in cm 






PHASE 


= Phase lag of the north tide behind 
the south tide, in hr 






DPHASE 


= Phase lag in east-west direction, 
in hour per Ax 






PERIOD 


= Period of tide at entrances, in hr 






TSHIFT 


= Time shift (in hr) for tide to agree 
with EST time 


#11 


2 


FCOR 


= Coriolis factor = 2W sin(Iat), in 
sec" T e 












STAGE 1 


= MWL-MSL (datum for sounding), 
if MWL * MSL, in cm 


#12 


1 


THETA 


= Clockwise angle from North to the 
y-axis of grid work, in deg 


#13 


2 


ALREF 


= Horizontal reference length, in cm 
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TABLE 5. Input Data for ANCMN 
(Continued) 



Input 


Card 
Content 


Symbol 


Definition /Value 






ROSSBY 


= Rossby No* which controls whether 
advection is needed to account for 
in the equations of motion, zero or 
nonzero 


#11 


1 


RWEX 


= Number of hours between climatical 
data, generally it agrees with hour- 
ly loop 


#15 


1 


TZERO 


= EST v?hen the simulation run starts 


#16 


4 


RR 


= Density of water = 1.0 






AV 


= 0.002 (AHREF) 4/3 , AHREF = refer- 
ence depth in cm 






BV 


= Same value as A 






BH 


= 0.002 (ALREF} 4/3 


#17 


1 


TINIT 


= Initial temperature, a constant for 
the whole domain 


#18 


A dec 
matric 
quirec 


k of cards or 
es MAR, ELEV 
J only if LN = 


data file, called AMATN, to specify 
, MEX, MEY, MX and MY. It is re- 
1. 


#19 


9 


NBLOK 


= To reset data block number 






TOTT 


= To reset total time in sec if 
necessary 






DT 


= To reset time step AT if necessary 






EST 


= To reset EST 






AMPLIT 


= If tide changes 






PHASE 


= New phase lag 






DPHASE 


^ New phase lag per ^X 






PERIOD 


= If tide has different period 
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TABLE 5. Input Data for ANCMN 
(Continued) 



Input 



Card 
Content 



Symbol 



Definition /Value 



TSHIFT 



= Time shift of the new tide 



#20 



A deck of NCY cards, each card contains six hourly 
weather data: 



TAIR 

HUMID 

WIND 

WDIR 

SRAO 
TSURF 



= Air temperature, deg C 

= Relative humidity in fraction of unit 

= Surface wind speed, cm /sec 

= Wind direction, from which direc- 
tion wind is blowing 

= Cross solar radiation, BTU/(ft 2 day) 

= Ambient surface water temperature, 
deg C 



f 
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SECTION 7 
PLOTTING PROGRAM 



This section presents the descriptions of main plotting program 
PLOTMN and its subroutines. As mentioned earlier, the plotting and 
analyzing of the results constitute Part 2 of the three-dimensional, free- 
surface model* Here, the tape containing the hourly results of simulation 
is the main input* The control data cards help the user to choose the hour, 
the plot and the comparison. The output is in a plot tape which is used 
by a CALCOMP plotter to generate plots. 

DESCRIPTION AND FLOW CHART OF PLOTMN 

The purpose of PLOTMN is to read in the measured and /or calcu- 
lated temperature fields and to plot the isotherms. In addition, the cal- 
culated velocity field is plotted in o-planes (o-axis), in certain x-cross 
sections (y-axis) and in y-cross sections (x-axis) , and the surface 
elevation field is plotted in contour plots. 

Since the main input is the data block from Part 1, the symbols 
and their dimensions agree with those that appear in ANCMN. In order 
to store the IR scanned surface temperature at four tidal stages for later 
comparison with calculated results, a TIR (IM, JM, 4) matrix is added. 
It is to be noted that temperature fields are interpolated at half-grid 
points from the mosaic IR images by hand. Several data cards contain- 
ing the quantities to be used in plot caption are read in as well as con- 
trol cards which assign the data block to be used (NPLOUT) and the plot 
to be done (IPLOT). A flag NSTAND is to assign which measured tem- 
perature field is to be compared with the calculated. The algorithm for 
PLOTMN is simple and straightforward since no complicated calculation 
is involved. The only calculation is to compute average deviation of 
calculated temperature field from measured temperature field at the same 
tide stage. The average deviation is given by 

.1, (TB(i, j, 1) - TIR(i, j)) J 

J,=!J =Tirn 



where TB is the calculated temperature while TIR is the measured tem- 
perature by infrared, and .-.(»# j) is the number of surface half-grid 
points in the domain. l *' 
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The isotherms of T1R and TB are the plots of main concern, as one 
is to be compared with the other in order to assess the accuracy of the 
model in predicting the hydrothermal dispersion of waste heat. Occasion- 
ally, the water surface contour, the surface current and the velocity 
profile are also of interest, as they depict the circulation set up by the 
tide and the wind in conformity with the configuration and bathymetry 
of the waters. The surface elevation contour is done by the subroutine 
ECHKON which is also used for plotting isotherms, but care must be 
taken to assign the contour values, since the surface elevation changes 
with tide; thus, the hard copy printout of surface elevation ETA must 
be consulted in order to choose the right contour values. The surface 
current is done by PLOTUV, while the velocity profiles in j(y)-cross 
sections and i(x)-cross sections are plotted by PLOTUW and PLOTVW, 
respectively. It is noted that the velocity scale for horizontal components 
is different from that for vertical component. The ratio is to remain the 
same as the ratio of horizontal to vertical length scale. Therefore, the 
velocity profiles are exaggerated in vertical direction; however, since 
the horizontal velocities of the top level (a ~ 0) are plotted right on the 
water surface, they show free-surface profiles too. Note that the j-sec- 
tions and i-sections are fixed by given I- J- values. They are chosen 
by the user's concern about the effect of plan-form configuration and 
bathymetry on currents. 

The flow chart for the main program PLOTMN is presented in Figure 
4, and the subroutines are listed in Table 4 for quick reference. It is 
to be noted that several CALCOMP subroutines are also listed. 

SUBROUTINES 

ECHKON 

This subroutine calls subroutines CONLIN and ENDER. This program 
was developed by the National Hurricane center for map contouring using 
CALCOMP or MILCO-type plotter. ECHKON is the entry point for the 
package. It scans the rectangular gridded scalar field, such as surface 
temperature or surface elevation, to determine where to start a new con- 
tour. Each contour is done in a loop. Inside the loop the subroutine 
CONLIN is called to do the interpolation and drawing, and ENDER is 
called by CONLIN to label each contour of the same contour v?lue. The 
exit of contour loop in ECHKON is made when the final contour value 
increased by increment has reached the specified maximum. Here, 
ECHKON is used for contouring TlR(i, j, N = 1, 4), TB(i, j, 1) and 
ETA(i, j). 

PLOTUV 

It computes the horizontal resultant velocity from components u and 
v at each level. In general, there are four levels corresponding to k = 1 
to 4. However, for the present problem, only the surface current is of 
interest; therefore, KPLOT is set to 1. 
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Read IPLOIR, IPLOCM, NIR, KPLOT, NCY 
Read Contour Values for TIR and Caption Data 
Read Control Vectors IPLOT, NPLOUT, NSTAND 
Read Contour Values for ETA and Plot Size XL (in) 



IPLOI 



NIR Loops 



R = 



Yes 



No 



ECHKON: Plot Isotherm for Measured Temperature at Each Stage 



1 



IPLOCM = 



NCY Loops 



No 



Yes 



READT: Read in Data Block from Tape 
Check Whether NPLOUT = and NSTAND = 



No 



Point by Point Comparison of Computed T with Measured T 



IPLOT (1) = 

No {- — - 



IPLOT (2) = 



No 



I PLOT (3) = 



No H- 



I PLOT (4) = 



No 



IPLOT (5) = 



Yes 



«. PLOTUV: Do UV Plot 



T 



Yes 



PLOTUW: Do UW Plot 



Yes, 



PLOTVW: Do VW Plot 



Yes, 



ECHKON : Do n Contours 



No 



Yes J 



ECHKON : Do T Contours 



Stop 



Figure 8. Flow chart for plotting program PLOTMN and its subroutines 
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TABLE 6* Subroutines Required In Main Plotting Program PLOTMN 



No. 


Name 


Description 


Remark 


1 


FACTOR 






2 


PLOTS 


Are CALCOMP Subroutines 


In UCS * ACALCOMP 
of UNIVAC 1100 


3 


PLOT 






4 


ECHKON 


Subroutine for plotting isotherms 
and contour of surface elevation; 
also draws the domain 


Calls subroutines ENDER, 
CONLIN, OUTLIN 


5 


REAOT 


Same as READT in ANCMN, for 
read in-stored data from tape 




6 


PLOTUV 


Plot U, Von different layer, to 
select the layer one can choose 
KPLOT value; normally it is the 
surface layer 


Calls subroutine OUTLIN 


7 


PLOTUW 


Plot U, W on chosen J sections 
west-east across the bay 


J section once chosen 
is fixed 


8 


PLOTVW 


Plot V, W on chosen 1 sections 
south-north across the bay 


1 section once chosen 
is fixed 


9 


CAPTN 1 


Write common heading on each 
diagram 




10 


CAPTN2 


Write title for UV Plot 




11 


CAPTN 3 


Write title for IR-T isotherms 




12 


CAPTN 4 


Write deviation of calculated 
temperature from IR-T 




13 


CAPTN 5 


Write tidal stage on the 
diagram 




14 


CAPTN 6 


Write title for UW plot 




15 


CAPTN 7 


Write title for VW plot 




16 


CAPTN 8 


Write title for surface eleva- 
tion contours 
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TABLE 6« Subroutines Required in Main Plotting Program PLOTMN 

(Continued) 



No. 



Name 



Description 



Remark 



17 


CAPTN9 


18 


ENOER 


19 


CON LIN 


20 


OUT LIN 


21 


FIT 



22 

23 
24 
25 



VECT 

AROHD 
NUMBER 
SYMBOL 



Writ* title for calculated surface 
isotherm 

Write contour value to label 
contour 

Called upon to draw individual 
contour 

Called upon to draw outline of 
the computational domain 

Fit a parabola to three point 
used in PLOTUV to interpolate 
water depth 

Calculates the velocity and 
calls AROHO to draw the vector 



Are CALCOMP subroutines, 
called upon in various sub- 
routines used in PLOTMN 



Subroutine of CONL1N 



Subroutine of ECHKON 



Subroutine of ECHKON 
and FLOTUY 

Subroutine of PLOTUV 



Subroutine of PLOTUV 



In library file UCS * 
ACALCOMP of UNIVAC 
1100 
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PLOTUW 

tt computes the resultant velocity based on components u and w. 
As mentioned, the vertical and horizontal components cannot be made to 
the same scale; therefore, the velocity profiles are distorted. The j-cross 
sections on which the velocity is computed and plotted are preassigned in 
the subroutine itself. For the present problem, plots Mre drawn for ] = 
4, 8, 12. All three plots are done on the same sheet* 

PLOTVW 

It computes the resultant velocity based on components v and w, 
and plots the velocity vectors on the vertical cross section along y- 
direction. The i values are likewise preassigned in the subroutine* 
Here plots are drawn for i » 4, 8, 12 and are on one sheet. 

INPUT DATA 

Table 7 lists the data to for PLOTMN. Free format is used generally. 
The total input consists of data cards and data files, but in the list, 
the input data have been numbered in the order of their appearance, 
regardless of whether card form or file form is used. 
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TABLE 7. Input Data for PLOTMN 



Input 



Card 
Content 



#1 



#2 
*3 

#4 

#5 

#6 

#7 

#8 



#10 

#11 
#12 



NIR 
NIR 
NIR 

NIR 

NIR 
NIR 

NIR 




#9 NIR 



NIR 

1 
5 



IPLOIR 
IPLOCM 

NIR 
KPLOT 

NCY 
TL(I) 
TH(I) 
TI(I) 

TAI(I) 

QKI) 
Q2(l) 

Q3(l) 

Q4(l) 

Q5(l) 

Q6 
I PLOT (I) 



Definition /Value 



* C, If no isotherm of IR obtained tempera- 
ture is desired 

* 1, otherwise 

* 0, If no isotherm of computed temperature 
is to be plotted 

* 1, otherwise 

= Number of IR obtained temperature fields 

s Number of the o-level to be plotted 

= Number of simulation hours 

= Array of tnin. contour values for TIR 

= Array of max. contour values for TIR 

= Array of increment values for TIR con- 
touring 

= Array of ambient temperature, assigned 
to out-domain points 

s Array of EST time, for caption 

= Array of wind speed values, for caption 
on TIR plot 

= Array of wind direction values, for cap- 
tion on TIR Plot 

= Array of air temperature for caption on 
TIR plot 

= Array of discharge temperature for cap- 
tion on TIR plot 

= Oischarge flow rate, for caption in general 

* Array of integers to assign the plot desired 
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TABLE 7. Input Data for PLOTMN 
(Continued) 



Input 


Card 
Content 


Symbol 


Definition /Value 


#13 


NCY 


NPLOUT(l) 


* Arrpy of integers to assign the hour de- 
sired to be plotted 


#11 


NCY 


NSTAND(l) 


» Array of intergers to assign the TIR to 
be compared with 


#15 


NIR 


ETAL(I) 


= Array of min. contour values for ETA 


#16 


NIR 


ETAH(I) 


= Array of max. contour values for ETA 


#17 


NIR 


ETAINT(I) 


= Array of increment values for ETA 
contouring 


#18 


1 


XL 


* Plot size in x-dlrection (in) 


#19 


2 


DX, DY 


» Spacing in x- and y-dlrection (cm) 


#20 


Data File 


1, J 


s The (i, j) value of boundary nodes, used 
for drawing domain's boundary 


#21 


Data File 


TIR(IJ,N) 


= To read in the IR measured temperature 
fields of different tidal stage, flies like 
HDATA, ED AT A, LDATA, FDATA 
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APPENDIX A 
EXAMPLE CASE 



INTRODUCTION 



The present model has been successfully applied to thermal disper- 
sion study at Anclote Anchorage. The Anchorage is located on Florida's 
Culf coast and north of St. Petersburg (Figure 9). It is a relatively 
shallow passage between the mainland and the Anclote Key. A series of 
barrier islands separates the anchorage from the Culf of Mexico. Through 
natural channels to the north and south of the Keys, the Anchorage has 
an unrestricted exchange of water with the Culf. 

The Anclote power plant operated by the Florida Power Corporation 
has two 515 MW, oil-fired, electrical generating units. Once-through 
cooling water is drawn from the Anclote River through a man-made canal* 
The six pumps delivering a total of 1,990,000 gpm (125*6 m 3 /sec) are 
designed to raise the water temperature 2. 8°C above the ambient. The 
heated water is discharged back into the Anchorage through the dis- 
charge canal with a dredged submarine extension. The designed total 
flowrate is approximately 53 times the long-term average flowrate of the 
Anclote River. At present, only Unit 1 is operative while Unit 2 is 
still pending permission. That is, the present flowrate is 62.8 m 3 /sec 
(995,000 gpm) . 

The principal driving mechanism for current circulation is tidal 
flux at the north and south entrances of the Anchorage. The tide is 
predominantly semidiurnal with mean range of 2 feet. Earlier measure- 
ments of temperature and salinity indicated the currents flow in and out 
through both entrances; however, the exchange appears to be stronger 
in the south than in the north, or the currents generally flow north 
during food tide and south during ebb tide. Moreover, the wind plays 
an important part too. The surface current direction depends on wind 
blowing at wind speeds exceeding 15 mph. 

The model as applied to the Anclote Anchorage shows its capacity 
of considering the effects of geometry and bathymetry, spatio-temporal 
variation of the free surface, various boundary conditions, including 
tides of different phase and range, surface heat transfer based on equi- 
librium temperature concept, and changing meteorological conditions. In 
addition, turbulence has been considered by using the eddy transport 
concept, and the effects of baroclinicity have been included. Again, the 
user should refer to Tuann et al. (1979) for the general review, mathe- 
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Figure 9* Anclote Anchorage location in the state of Florida 
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matical formulation, finite difference implementation and numerical method 
of solution. 

The finite difference grid work is three dimensional and is designed 
to cover the area of interest. The grid size Ax x Ay x A a is that the 
(east number of vertical layers is four. That is, a 4-Iayer, 5-level, 
vertical partition is a reasonable choice for the present generation of 
computer, The grid work is allowed to orient away from north-south, 
east-west system, but in general, the x-axis of the grid system aligns 
with west-east, and y-axis with south-north. Thus, the subscript i 
increases eastward, while j increases northward. The z-axls is chosen 
upward from mean water surface, while the subscript k increases down- 
ward from the water surface. That is, the k * 1 level is the free 
surface which is continuously changing, while k = 5 is always the bottom. 

For the study of Anclote Anchorage, the grid is 16 x 14 x 5 with 
five levels, each with 224 nodes for a total of 1120 nodes. The grid size 
used isAx=Ay = 417m. Depths off the natural coastal line are read 
from the Coast and Geodetic Survey chart. The maximum depth is 4 
m at the south end of Anclote Key. It was found that gravity waves 
were the dominating consideration with regard to the maximum allowable 
time step At. A 15 second magnitude of At was found to work well for 
the present grid system. 

Numerical results were obtained with the University of Miami Com- 
puting Center UNIVAC 1100 computer. The time histories of the three 
velocity components, (u, v, w) , the surface elevation, n, and the tem- 
perature, T, for a 24-hour simulation period were obtained with about 90 
minutes of computer time in most cases. This is a time ratio of about 
16:1 (the ratio of real time to computer time). 

PROBLEM STATEMENT 

Florida Power Corporation has a fossil fuel power plant situated at 
Tarpon Springs on the Anclote Anchorage. The discharge rate is 62.8 
m 3 /sec of water at temperatures, in general, 2. 8°C above the ambient 
water. On June 19-20, 1978, a team carried out an in-situ data acquisi- 
tion mission to gather field data on temperature and current. At the 
same time, four flights by NASA/KSC were undertaken to obtain tempera- 
ture by remote sensing method. These four flights were intended to cover 
four different tidal stages in the Anchorage. The remotely sensed data 
were processed into digicolor film. The in-situ measurement of surface 
water temperature at the time when the airborne IR data was undertaken 
provides a reference for IR temperature. With this reference, the iso- 
therms were drawn from the digicolor film. The in-situ current measure- 
ment data were used in plotting current of different depth. The ground 
measured temperatures were used to draw surface and subsurface iso- 
therms. 

Once the model has been verified for its versatility and its capacity, 
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Figure 10. Grid work for the Anclote Anchorage 



and in particular, to the prediction of hydrothermal development in well- 
mixed shallow coastal waters, the model was run with actual tidal and 
meteorological data as input, but the initial temperature condition either 
could be that of uniform state, that is to assume the power plant start 
impulsively, or could be that of an IR temperature field, that is using 
IR data as the initial temperature field. The demonstrative runs were 
carried out to simulate the hydrothermal situation for several days. The 
predicted current fields are verified against the in-situ measured currents, 
and the predicted isotherms are verified against the IR-obtained surface 
isotherms and the in-situ measured subsurface isotherms. 

CALCULATION OF PARAMETERS AND INPUT DATA 

In this section, the specification grid system, reference and physi- 
cal parameters, tidal and meteorological data, discharge and intake velo- 
cities, ambient and discharge temperatures will be presented- The actual 
calculation of some input data quantities is carried out in detail for the 
purpose of demonstration. 

Grid System 

The map indicating the exact locations of power plant, intake and 
discharge outlets, and the sounding of the Anchorage was used to deter- 
mine the size of the domain, the grid system to cover it, and the arrange- 
ment of intake and discharge points in the system* So, a domain of 
about 6 km x 5 km covering most of the Anchorage was used. A grid 
system of 16 x 14 was selected in the horizontal plane* The size of the 
grid cell is Ax = Ay = 416.7 m; this size and th* grid orientation has 
made the intake and discharge outlets to the open water fall in with 
nodes respectively, and the intake and discharge channels have 45° and 
315° orientation respectively. The depth was specified according to 
sounding chart. There are five nodes in the vertical direction. This 
gave a total of 16 x 14 x S nodes. The coordinate system and grid 
work are shown in Figure 10. The MAR matrix, bottom elevation matrix 
and four additional marker matrices are stored in data file AMATN. 

Reference Quantities 



L: Reference length = ALREF = Maximum Length = 6 km 

Horizontal eddy diffusivity = 0.002L 4/3 - 
0.002 x (600,000) 4/J = 100,000 cm 2 /« 



BH 



/sec 



BV: Vertical eddy .diffusivity = 0.002H 4/3 (H = maximum depth) = 



Vertical eddy .diffusivity = 0.( 
0.002 x (360) 4/,a = 6 cm 2 sec 



For shallow well-mijced tidal water about three times the calculated 
value was found suitable. Here, we use BV = AV - 20 cm 2 /sec, i.e. 
Turbulent Prandtl No. = 1 
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RR: Density of water = 1.0 

THETA: 0, since the grid system orients along North-South 

RWEX: 1.0 

ROSSBY; 0. 0, that is, based on test run, it was shown that the nonlinear 
inertia terms can be safely neglected to save computation 

TJNIT: Initial uniform temperature or reference temperature « 20 deg C 

Calculation of Time Step, DT 

In order to determine the time step, DT, the stability criterion has 
to be followed, which is done as follows. 

nT „ DX . 41760 _ - n w 

DT JTg* " JixyHdxiiii " 50 sec 

About 1/3 of this value is reasonably safe to use. Here we use DT = 15 sec. 
Calculation of Intake and Discharge Velocities 

1. Flowrate = 995,000 gpm (from power plant physical data) 

= 62. 8 m 3 /sec 

= 62.8 x 10 cm 3 /sec 

2. Both intake and discharge canals are at 45° from N, therefore 

31.4 x 10 cm 3 /sec is crossing the Ax and Ay at the point of 
intake and discharge. 

3. The average depth at intake and discharge is approximately 4* or 
122 cm, and the width is Ax = Ay = 41760 cm, so the cross-sectional 
area is 41760 x 122 cm 2 . 

4. ~ <* average velocity is: 



U = V 
ave ave 



31,4x10' 



41760x122 
5. The velocity profiles are assumed as shown 



= 6.163 cm /sec 



1 
2 
3 
4 
5 



7.0 cm /sec 
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6. To allow for canal storage during tide change we assume the intake 
and discharge velocities to be sinusoidal, i.e. 

Intake: V3(14, 4, k) = 7 - 3 x cosl^CEST - 7.625)] 

U3(15, 3, k) = V3(14, 4, k) for k = 1, 2, 3, 4 

Discharge: V3(14, 8, k) = 7 - 3 x cos[~g(EST - 7.5)] 

U3(15, 8, k) = -V304, 8, k) fork = l, 2, 3, 4 

where 7.625 and 7*5 are taken to be the phase shift which takes 
into account the time to travel from the south end of Anclote Key 
to the concerned point. 

Calculation of Tide on June 20, 1978 

Simulated diurnal tide is shown in Figure 11, where 

1. Period = 12.5 hr 

2. Stage = short term average sea level - MSL = 48 cm 

3. Amplitude = i short term average tide range = 65 cm 

4. Time shift = 7.125 hr 

i.e« at 7,125 a.m., June 19, 1978, the tide at the south end of 
Anclote Key was zero. 

5. W - E lapse = 0.014 hr/DX 



Wave propagation speed C = J2gf\ = J5c980x360" = 850 cm /sec 
(H = 360 cm is the maximum depth of the Anchorage.) 

The time needed to travel one grid distance is 

™=<J™ = 50sec = 0.0Uhr. 

We use 0.014 hr per DX for phase shift in VV - E direction and the 
imposing tide at the south entrance is 

n$ = 48 + 65 sinl—^CEST - 7.125 - 0.014(1 - 1)] 

I = grid no. in W - E direction. 
6. S - N lapse = 0.15 hr. 

Distance from south entrance to north entrance is about 543,000 cm. 
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V* 




Tide at uouth end of Anctote Key 



,2 m 



Simulated tide for calculation n= 4B ♦ 65 x ain (j-275 (*"♦)) 

/ — -Average Level 




Figure 11. Semidiurnal tide for June 19-20, 1978 at south end of Anclote Key 



Time for wave to travel this distance is iff = 0.18 hr. We take 
0.15 hr as phase difference between the south and the north boun- 
daries; there, the imposing tide at the north entrance is 

n = 48 + 65 sintj^(EST - 7.125 - 0*15 - 0.014(1 - 1}] 

Calculation of Anclote River Flowrate and Temperature 

1. The tstance traveled from South Anclote Key to Tarpon Springs Is 
20 OX. We estimate a time lapse of 0.5 hr to account for the re- 
tardation due to buffering effect of river storage and Anclote River's 
natural outflow. 

2. The average current is estimated to be 20 cm /see, therefore, we take 

U3(16, 1, k) = 20 cosl-y^jlEST - 7.625)] 
V3(15, 1, k) =-20 cost^tEST - 7.625)1 

for k = 1, 2, 3, 4. 

3. The surface elevation at Tarpon Springs is to be calculated. 

4. To be in accordance with given velocities at Tarpon Springs , the 
temperature there is also assigned and its value has a 24 hr period 
instead of 12.5 hr. This temperature is 

T305, 1, k) = 26.9 + 0.5 sin[|j(EST -12)] 

where the 12 hr shift is to make the peak temperature occur at 1800* 
Thus, the water in and out at Tarpon Springs has a temperature 
ranging from 26. 4 (before dawn) to 27. 4 (late afternoon) • 

Discharge Temperature and Gulf Water Temperature 

1. On June 19-20, 1978, the recorded discharge temperature at daytime 
is in the range of 29.3-30.3. To account for the further drop of 
discharge temperature due to cooler ambient temperature at nighttime, 
we a?iume a sinusoidal variation of discharge temperature with diurnal 
period* 

2. Discharge temperature is estimated 

T3(14, 8, k) =29.4 + 0.4 s.n[~(EST - 12)] 
therefore, the highest discharge temperature of 30. 3°C occurs at 
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6 p.m. and the lowest (29. 3°C) at 6 a.m. 

3. The Gulf water outside the Anclote Anchorage as well as the atmos- 
phere is sink to the heat disposal from the power plant; therefore, 
the boundary condition on temperature at the north and south en- 
trance is not considered as adlabatic as in normal case of far-field 
thermal pollution problem. Instead, we specify the outside-anchorage 
ambient temperatures. Again, they are 24 hr periodic and their 
values should be in accordance with the measured temperature in the 
same neighborhood. Here in compliance with measured data, we use 

T ab * 27#0 * °* 2 slnt f5 (EST " 12)1 

for both ambient temperature outside the south and the north boun- 
daries. 

EXECUTION DECKS FOR CALCULATION AND PLOTTING RUNS 

The following execution decks are for use in UN I VAC 1100 computer 
at the University of Miami. These may have to be modified if a different 
computer is used. The programs and subroutines used in these runs are 
all compiled and stored in the file. 

Calculation Run 

First Run~ 

1. 9 ASC, A FILENAME. 

The file 'FILENAME 1 is assigned for the run. 

2. ASC, T 8., 16N, TAPENAME1 

A tape file names '8 1 is being assigned. The tape is 9-track, and 
the reel number is TAPENAME1.' 

3. PRT, S FILENAME. ANCMN 

The main program "ANCMN 1 is printed. 

4. PACK FILENAME. 

'FILENAME 1 is packed together, eliminating the space left by deleted 
elements, and thus, condensing the file. 

5. « PREP FILENAME. 

Prepare an entry point table for the 'FILENAME. 1 

6. MAP, S 

Combines relocatable elements to form an executable absolute element. 

7. IN FILENAME. ANCMN 

8. LIB FILENAME. 
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9. END 

10. 9 XQT 

11. 3 

Case number (NCASE) . 

12. 1 

First run (LN). 

13. 1 

Store the calculation results on to tape TAPENAMET (KSTORE). 

14. 1 

Calculate velocities (KVEL). 

15. 1 

Calculate temperatures (KTEMP). 

16. 

Specify numbers of the latest hour of the last run (MBLOCK). 

17. 13 

Number of hours to be simulated (NCY). 

1 8. 3600 . 

Print the results at each 3600 seconds (TPRT) . 

19. 41760., 41760., 0.25 

Grid sizes in x-, y- and ^direction (DX, DY, DS). 

20. 15.0, 0., 60., 0.08, 0.005, 12.5, 7.125 

Specify time step, the difference of average tidal level and mean 
water level, amplitude of tide, north-south phase lag, east-west 
phase lag per DX, period of tide, and time shift (DT, STAGE, 
AMPLIT, PHASE, DPHASE, PERIOD, TSHIFT). 

21. 0. 66E-4, 48.8 

Coriolis factor and the difference of mean water level and mean sea 
level (FCOR, STAGED. 

22. 0. 

The y-axis coincides with North (THETA). 

23. 8.E5, 0. 

Horizontal reference length and Rossby number (ALREF, ROSSBY) . 

24. 1. 

Number of hours between climatic data (RWEX). 
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25. 7. 

Simulation run starts at EST 0700 (TZERO). 

26. 1.026, 20., 20., 50000. 

Specify wattr density, vertical eddy viscosity, vertical eddy diffu- 
sivity, horizontal eddy diffusivity (RR, AV, BV, BH). 

27. 27. 

Initial water temperature (TINT). 

28. 9 ADD FILENAME. AMATN 

Input data file •AMATN 1 for specifying grid matrices and initial 
water depth. 

29. 9 ADD FILENAME. FDATA 

Input data file 'FDATA* for initializing temperature distribution. 

30. 9 ADD FILENAME. C2007 

Input data file 'C2007' of climatic data. 

31. 9 FIN 

Terminate this calculation run. 

Sebsequent Run— 

1. 9 ASG, A FILENAME. 

2. 9 ASC, T 7., 16N, TAPENAME1 

A tape file named '7' is being assigned; the reel number is 'TAPE- 
NAME I. 1 This tape was used in the first run for storing the hqurly 
calculation results for 13 hours in 13 blocks* 

3. 9 MOVE 7., 12 

Move TAPENAME1 to the 13th block which is the last hour result of 
the first run and is going to be used as input data for this subse- 
quent run. 

4. 9 ASC, T 8., 16N, TAPENAME2 

A new tape named TAPENAME2 1 is assigned to store the calculation 
results of this run. 

5-13. Same as the cards 3-11 of the first run. 

14. 2 

Continuing run (LN). 

15-17. Same as the cards 13-15 of the first run. 

18. 13 

The last hour of the first run is 13. 

IS- 29. Same as the cards 17-27 of the first run. 
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30. Same as the card 20 of tha first run. If new tidal data Is needed, 
this card has to ba changed. 

31. 9 ADD FILENAME. FO AT A 

If diffarant IR tamparatura distribution Is neadad, FDATA has to 
ba changad. 

32. 9 ADO FILENAME. C 2007 

Tha data fila C2007 has to ba changad sinca the weather condition 
will be different from the first run. 

33. 9 FIN 
Plotting R un 

1. 9 ASC, A FILENAME. 

2. ASC, T 7-, 16N, TAPENAME1 

A tape file named *T is being assigned. 'TAPENAME1 1 stored the 
results of the calculation run. 

3. ASC, T 11., 16, TAPENAME2 

A tape file named '11' is being assigned. The tape is 7-track, and 
the reel number is TAPENAME2.' This is used for plotting tape. 

*. 9 PRT. S FILENAME. PLOTMN 

5. 9 PACK FILENAME. 

6. 9 PREP FILENAME. 

7. 9 MAP, S 

8. IN FILENAME. PLOTMN 

9. LIB FILENAME. 

10. LiB UCS'ACALCOMP. 

Call 'CALCOMP* plotter Horary. 

11. ENO 

12. 9 XQT 

13. 1. 1, 1, 1, 6 

Plot IR isotherms and computed isotherms; only one IR temperature 
field is to be plotted. Plot isotherms on the surface lev?l only, 
rnd run for 6 simulation hours (IPLOIS, IPLOCM, NIR, KPLOT, 
NCY). 

14. 27.5 
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Minimum contour value for IR plot (TL). 

15. 30.0 

Maximum contour value for IR Plot (TH). 

16. 0.7?. 

Increment of contour value for IR plot (Tl), 

17. 27.0 

Ambient temperature (TA1). 

18. 13. 

EST for caption (Ql). 

19. 358.0 

Wind speed in cm /sec for caption (Q2). 

20. 110. 

Wind direction for caption (Q3). 

21. 29. 4 

Air temperature for caption (Ql). 

22. 29. 5 

Discharge temperature for caption (Q5). 

23. 62.7 

Discharge flow rate in cm 3 /sec for caption (Q6). 

24. 1, 1, 1, 0, 1 

Plot UV, UW, VW velocities and isotherms (IPLOT). 

25. 0, 0, 0, 0. 0, 1 

Plot the results at the 6th hour (NPLOUT). 

26. 0, 0. 0, 0, 1. 1 

Compare the deviation of computed temperature from IR tempterature 
(NSTAND). 

27. 0.0 

Minimum contour value for surface height (ETAL) . 

28. 0.C 

Maximum contour value for surface height (ETAH). 

29. 0.0 

Increment of contour value for surface height (ETAINT). 

30. 6. 

6" plot size in x-direction (XL). 
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31. 42000*, 42000. 

Grid size (cm) in x- and y-direction (DX, DY). 

32. ADD FILENAME. APER1 

Specify boundary nodes for plotting the boundary. 

33. ADD FILENAME. ED AT A 

Input data file 'EDATA; 1 specify IR temperature distribution at 
ebb tidal stage. 

34. FIN 

Terminate this plotting run. 

The input data file AMATN, FDATA, C2007, APER1 and EDATA 
are listed in the Appendix B. If these data are not stored in the 
'FILENAME/ card decks have to be substituted. 
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SAMPLE OUTPUT (EST 1300, June 20 # 1978) 
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1*415. 
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1*4*5. 
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I 8b IS. 
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ii/i*. 

in mh. 
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Figure 12. Surface velocity, Anclote Anchorage by modeling 
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TItlEC JUNE 20. 19*78 )« 
NINO SPEEOCCH/SEC)i 
MINO 0IRECTI0NC0EO/NH 
H IR TEllPERRTUREJ OEG-C ) * 
OISCHRROE TEHP(OEO-C1» 
OISCH FLOWRRTEC CUM/SEC )» 
LENGTH SCRLEdCfls X CM)) 
VELOCITY SCflLE(CH/SECJ» 



13.0 
356 .0 

1 10. 
29.4 
29 .5 

62.7 

41339. 

52.49 



J« 12 




j* a 




Js 4 



EBB TIDE 



Figure 13. UW velocity, Anclote Anchorage by modeling 
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TIHECJUNE 20. 1978)* 
KINO 3PEE0(Cn/8EC)t 
HINO 0IRECTIQN(0EO/N)i 
AIR TEHPERRTUREIOEO-Ot 
0ISCHRROE TE«P(0EO-C)i 
0I8CH FLQHRRTE( CUM/SEC) * 
LENGTH SCRLEdCfl* X CK)> 
VELOCITY SCRLECCH/SECJi 



13.0 

358 .0 

110. 

29.4 

29.5 

82.7 

41339 

52.49 
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Is 8 




Is 4 
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Figure 14. VW velocity, Anclote Anchorage by modeling 
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TIHEUUNE 20. 1970) t 
UINO SPEE0(Cf1/SECJt 
UINO 0!RECT!QN(0EO/N)t 
RTR TEMPERRTURE(0EO-C)« 
0ISCHAR0E TEf1P(0EO-CJt 
OISCH FL3MRRTEC CUM/SEC)* 
LENGTH SCALE (ICRs X CH)i 
VELOCITY SCRLE( CM/SEC J t 



13-0 

356 .0 

110. 

29.4 

29.5 

62.7 

41339. 

52.49 




OEVIATIQN FROd 1R TEMP* 
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Fiqure 15. Surface temperature, Anciote Anchoraqe by model inq 



JMfc^fcfe^tt, >-> -^— _^_-—. .*mm 



+ * 



APPENDIX B 
FORTRAN SOURCE PROGRAM LISTING 

LIST OF SUBROUTINES OF THE MODEL 
Calculating Part 

1. Main Program 
ANCMN 

2. Subroutines Called (in order) 

BAYBOT: Reads grid matrices and bottom topography 

BAY1NI: Specifies initial conditions 

READT:. Reads data from tape for continuing run 

IRREAD: Reads 1R data as initial temperature distribution 

EQTEMP: Calculates equilibrium temperature 

BETA: Calculates surface elevation and vertical velocity in 
xya coordinates 

BNRTIA: Calculates inertia terms in momentum eq. at interior points 

ABNR3: Calculates inertia terms in momentum eq. on the north and 
south boundaries 

BVELS: Calculates interior velocities 

ASAF3: Calculates north and south boundary velocities 

GIVENU: Specifies velocities at discharge point and river mouth 

CONV: Calculates convective terms in energy eq. 

TCOMPT : Calculates interior temperatures 

GIVENT: Specifies temperature at discharge point 
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WCAL: Converts vertical velocity in xyc coordinates into xyz 
coordinates 

ANCPR: Prints surface height, velocity and temperature at four 
locations at each time step 

TPRLOK: Main printing program 

STORET: Stores calculating results onto the tape 

ZZI: Finds the current direction 
3. Data Files 

AMATN: Specifies marker matrices and elevations 

APER1: Specifies outline of interest area 

C2007: Climates data on June 20, 1978, start at 0700 

HDATA: High tide data from IR 

ED AT A: Ebb tide data from IR 

LDATA: Low tide data from IR 

FDATA: Flood tide data from IR 
Plotting Part 

1. Main Program 
PLOTMN 

2. Subroutines Called (in order) 

PLOTUV: Plots U, V velocities on different levels 

PLOTUW: P!ots U, W velocities at different j sections 

PLOTVW: Plots V, W velocities at different i sections 

ECHKON: Plots surface isotherms and surface height 
2NDER: Subroutine in ECHKON, for labeling 

CONLIN: Subroutine in ECHKON, for contouring 

CAPTN1: Writes captions on the plot 

CAPTN2: Writes captions on the plot 
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CAPTN3: Writes captions on the plot 

CAPTN4: Writes captions on the plot 

CAPTNS: Writes captions on the plot 

CAPTN6: Writes captions on the plot 

CAPTN7: Writes captions on the plot 

f PTM8: Writes captions on the plot 

CAiTN9: Writes captions on the plot 

FIT: Fits a parabolar to three points 

VECT: Establishes the components of a vector 

OUTLIN: Draws the outline of interest area 

* The plotting subroutines PLOTS, PLOT, AROHD, NUMBER, SYMBOL 
are existing in UNIVAC 1100, University of Miami, CALCOMP file. 
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SUBROUTINE LISTINGS 

*ft<iw< 11. ANON POP CREATED ON 1* CCC ?<l AT I0:*l:3* 

2 C MAIN PPfPSAM CF ?-0» FPE£-SU"FiCE ♦ HYCRQtHCPMAL *0QCL. FEATURES: 

7 C CTCS C CUFPFT-FRAKKEL PIF FE^EN CINC ; DIFFUSION CF T CONSIDERED; 

ti C NGNLINE/P TERMS CAN ?E INCLUDEP? A<S0CUT£0 WITH PLOT PPOGP AM -PLO TNN 

s <:*•*«***«• i*« ••»*••»•**+«*** •*.*•••*.*••••**♦**••••*•*•••**••*••*♦••**« 

f PAPA^fTCP I'.'Mb, Jvrt*,KN=5 , TH=15.JKri J 

7 OlfftSLON Ul(lN.JN,*N> t U?(lN t JN,KNl,U3<IN,jN,»<N) t VniN f JN,KN», 

• CV2U K,jN,KNI ,V3Utt»JN,HNI ♦ * UM t JH,KNl ,0M< JM , J" ,X N J f FT A 1<IN,JM1 t 

CE TA(l** f J«<) t ETA3nH,JH) t PX<lN f JN,KNl»RmN t JN f )tN|, MX (IN, JNJ, 

\r. fMYt I* V JN1,MAR<X**JN1,HIXN,JNI •hOClM 9 JH),ELEV(IN«JNI, 

U CUF <I*,J*,KNI fVBCI'-tJHtKNl ,HUI ?N,UNl t HV(lN f JN> ,M£X(IM«J*li 

1? CH£Y( l-.JM) 

13 dimehsIon Tnif.jsuNi , t: 1 i«. jh.kni ,t 3t im, jm»xni » tccxw»j«,kn> 

1 U DIME KS ION AVftt J») * ANGC J*| t 0TZ! X« f JH| t TE(-X M 9 J»,rtM 

IS OIMEMION *TX« XM,JMl f TIDE3M]H),TlCE2N<IMl,TlOCXmi*» 

18 DIMENSION WAT4XW t jMl t TIOE3StlMSTIOE2SlXH),TIOClitlM1 

1* WRXTE«6»2fe7) 

13 C 

1Q C*«1**CASE SO UNCERTAKING, FOR LABELLING PURPOSE 

:p REAO 2, NCASE 

21 C 

2* C**2**LN=1 FIPST RUN OF PRESENT CASE; M SUBSEQUENT RUN 

2? REAO 2, LN 

2* C 

25 C*»3***ST0RE=0 NO STQPEtTEST RUM, =1 STORE ON TAPC 

IH REAO 2. KSTORE 

27 C 

<o C*«*«*KVEL:0 NO V-CALCUIATION. ~l 00 V-CALCULA TION 

2* REAO 2, KVEL 

3C C 

31 C**5«*KTE*F=J NO T-CALCCLATXON, =t 00 T -CALCULATION 

Z2 REAO 2, KT£"P 

33 C 

3« C»*b**T0 ASSURE *8L0« CONTINUE, THUS AS COUNTER OF SIMULATION HOUR 

3* REAO 2, -ttLOK 

3* C 

t? C*«7**THIS PUN WILL 00 NCY*TPRT/3*»00 HOURS OF SIMULATION 

3° REAO 2, NCY 

7 o c 

uO C«*«**FOR HOURLY CYCLE TPRT=36C0; OTHERWISE CLIMOATA 00 LIKEWISE 

m! RFA0 2, TPRT 

M2 C 

n: C**<*«0X t D* GRIT SITE IN Crf; OS SPACING IN SEG^A DIRECT IONI . 251 

uu REAO 2, OX,QY,CS 

«s c 

M* C»»lC*«TIMCSTE e t**wL-TEMPO*AL KWL f 4nPLTT,N-S PHASE. OIFFERE'JCE , 

m7 C**lO**t-w phasl niFF PER DX, TIDEPEWTOO, TIDE SHIFT In HOUR 

<** REAO 2, OT # STAGE ,AM(>LIT f PHASE«OPHASE t PERIOO tTSHlFT 

5C C**ll**FC0erC0PIC1 IS F ACTOR f TEMPORAL MWL-ANNUAL MWLtREF. FOR SOUNDING) 

51 READ 2, FCOP f STAGEl 

5 ? C 

S T C**12**THC TArANGLE BETWEEN NORTH *N0 GRID Y-AXIS, POSITIVE CLOCKWISE 

5" READ 2, THETA 

55 C 

56 C*«l3**ALPeFiPEF'PENCf HORIZONTAL LENGTH IN CM t ROSSBY NO l£S TIM A TED I 
5' READ 2, ALREF,R0SS5Y 

5? C 

*i C*»l****QyE>rso OF HOURS BETWEEN yfATHFR CRStRV* TIONS CXU GENERAL H0URLY1 

h" RcAO 2, R-EX 

61 C 

fc? C«»15**T2ES0-EST OF THE C*Y <*HZH THE SIMULATION RUN STARTS 

6 T REAO 2t T:E o 

b« C 

«,e C*»16**QEN 51TY.VE0T £PQY VI SCOS IT Y , VEP T C HORI EOOY OIFFUSIVITY 

f> REAO 7, RR,AV,BV,QM 

6? C 

fcc C«*LT**IM TTAL TE**P FOR THE WHOLE COMPUTATIONAL DOMAIN 

L<< REAO 2, TIMTT 

7- C 

7; WRlTC(b,321 NCASE 

7? WRITE lb, !) LN 

7? IF IK 5TOSE.E0.0 ) wPJTF < fef 1 3Q l 

7* IF IK 5T01E.3T.C) W<? I T£ C C , 1 3 1 I 

7* UC FCRH/T lit, *OATA ?*nT RCCORflL" ON TiPEM 

76 131 FCR- *T( IX, TATA RrC^RTLD ON Upr.*) 

77 WRlTFfo.tJ'JCY 

7* WRlT£ib,S) ox f OY,r;s,oT 
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7* UftXTCf6 t 24 I STAG£»A«PLIT t PHASE»PPHASE»*-ERIOO»TSHIFT 

tT yftlTElb«7l FCOR 

51 yfllTHb, ICll) TM£TA 

B2 WfflUffrflSl ALPEF 

?? 15 fQ*Htuit t *n*zrz* t r\2.c f * cmm 

S* WRIT?I6,72) T7£PO 

6< URITEf*«73 t RP t AV t BV f 8N 

f* bflIH(6^06l QbfX 

!T 606 F0flPATUX t *RwEXT« t F10.2l 

!' , bflXTECb.SMI TJMt 

«« 3**» FOA»*JT<Il t 'U<i!TS a t F10.?} 

*« 12 FORH*T(U f *T2r^O = * 1 FtQ.2) 

Jl T3 FCR*JTUX t »*P. = » t flS.2 f * AVS*,F10.2»* BV=«,Fl0.2 t * PH=NF1Q*2) 

*! 0C=57l3 

<*» *2xk*-i 

9S IGO=l 

** THCT/rTMETI/CQ 

97 OUHX :2.*ox 

?« OUMt:?.*OY 

S« CUMS :2.*0S 

1CP IFCL*.GE«21 CO TO X 

1C1 EST=17ERQ 

1C2 TAUX:C. 

1C7 T«ur:2. 

1C<» TTOT:3. 

10* N8L0K=Q 

1 -8 CALL r UYBOTriN f JN f KN f IM.JN t U 1 ,U?.U 3t V 1 » V2 , V ?, W ,OM »£TA 1 1 CT A . ETA 3. 

1C rPX.PItMX.MY^lR.HtHPtELtV.UR.VBtHUtMV.MEX.HtY ♦ 

K* CDUMX , nu*Y»OUKS,DX,OY f OS»OT f FCOft t TAUX f f AUY ,G »HC ASE ♦NBLOK , I TOT , 

1C CSTAGC1) 

1 1 C C 

111 C*»18«»A OfC* OF AMATN CARDS OR FILE AMATN IS NEEDED HERE, (ONLY IF LNri)*** 

11? 'CALL PAY TNI rXN,JN f KN, I H , JM, Ul P U2 t U 3 f V 1 , V2 , V ?♦ V f OM t £TA I . ETA »ETA3. 

I 1* CRXtR^MXtMY.^AR.H.HPfELEV.U? , Vb ,HU , H V . Mf X , NFY t 

1 l e COUHX .PUMY 1 0UKS f OX»0Y t DS,Ot t FC0R#TAUX,IAUY t 6»NCASE»KBLOK,TTOT» 

IIS CTl # Ti t T3,tC.TB,TlNlT) 

1 17 00 SCI 1 = 1 ,IH 

1 I* TIDE 15 f T »=0 . 

11° TIDE :S(I):0. 

i:: soi tioe :sm::. 

121 00 5C2 1=1,1* 

122 TIOE !N(I >rr? * 
12 T TICE2N<n=0. 
12* SC2 TIDE 2N(I >:Q« 
12* GO T( 222 

125 C***«*****«* •'**•*•««•«•«*« 

127 1 CCNT jr:uE 

12* CALL 3 EAOT<TN,jN t KN t lM,JHfUl ,U2,Vl t V2,WtETAl,Tl t T2,VB , 

12' CE TA t »X,HY,«AR,H*H« f ua |VP ,HUtHV f f*EX ,HET ,NC A SE • N 6L0 K , TTO T ,0T , 

13? CT10E JN, TTCE'>.,TrDE3N 1 TIOElS,TTO£2S t TIOE3S«STAi;C»tST f 

171 CA-PL TT.PHASE »0PHA5£ t P£RI001 

12: IF(NELOK # NE**PL0KI GO TO 1020 

1 7 T C 

li« C**«**IF AKY OF FOLLOWING PARAMETERS NEED TO BE CHANGED. wrff£ 15 CHANCE 

l! c C*«19*«C0M BLOCK NO, CTNT TOTAL TlHE .0 T f E S T , T I PE 0A1A FOLLOW 

1 It READ ?»KBL0K,TT0T,0T f £5T,AHPLIT tPH ASE # CPNA $£, P£PIO0* TSHIFT 

1 !■ WPIT£tb,3M S T *C E, A -PL I T, PH A SE,CPHASE» PEP 100 »TSHI FT, OT 

l- fe WRITE<b t 33) EST f r*CASE t NPLOR 

1>! C*«**»1F 1B-T IS USED AS INITIAL T£PP, THEN RELIEVE FOLLOWING CALL 

JM2 C*«*»*THEN I DECK OF IR-OATA CARDS IS NEEDED IN 1TRCA0 SUPPOUTirt 

1" T 222 CALL IPREA0ilH f JM t KN t Tl t T2J 

1»» C 

1*S C 222 CONT J\U£ 

in* C •••••*••««*••*•*•••»•*•**•<•♦*♦• •*•**•••*«•••*•*♦* 

mr c*MMMi»»(rs loop in houocy.step. with hourly cli^atical oati 

im.« (:«*****•«•«*«"•*"••«*********••*****«*«*«••«•***« 

!«*« 00 1 CI JCTPtl,NCY 

1ST C 

151 C**2C**4 OtCH PF NCY CLI^ATICAL CATA CASC1 FOLLOWS 

IS"* C**2Q*-E*C> CARL PCCOPOS HP T£*P ♦H»J(- 10 I Ty ,. iPHQ .WPXP.SCPAD »JURF/C£ T£NP 

15? RtAO 2, TAlP t HLHI0 t WIN0tW0IPt5RA0 t T5U«F 

1 5«* C 

1 < r kRITC(6 f :29JTAlR f HlmXC.WtS0,W0rR, SPAOtTSURF 

15* 32* F 0«^ /T ( lY t *TAlR f HUM tO ♦blN0»WCIR«SPAO,T SURF »»6F 10,2) 

15Y C8=2« 
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IS* EP$ieNrrjorSM80. 1/57.3 

IS* 1 yPR=HN0/lC3. 

1*.*> IF <W INQ.LT. ICC. )CTeHr.CCl25*<WPO««.j| 

Ul IF( ( WpiO.&E . ICQ. I .AND. IWIN0.LT.15QC** I JCTCN=.0n05* SCRTI WPft I 

U: IF<V T*Q.Gt.lSC0.)CTENr.0C26 

16 T TAU=CTES«.3CI?:«(wrV0«»2) 

16* T*UX :TAU«SIN(fP*LCN-TMtT4) 

165 TAUT :T4U*C0S CtPSLON-TMtTAI 

166 CALL C0TCMPCTAlP f HUfIO,yi4O f WOIPtSfiAO»I5U»F»TO£tt t SK f TCQ! 
It 7 Txr 1PRT 



2Z1 



2 \ C 
* • » 
2 12 

: ; T 
rn 
2 i f 
rife 
r i t 



2;* 



>!1 



2Jfr 



C***«NOPTK C SOUTH AF8IENT TCNP IS SISE OF 2* HP PfRlOOt OP SPECIFIED 
TAeS :27.r>.2*<J*M.2618*CE$T-12.n) J 
TABNC?7.0*.2*5lW|.2ftia*<£ST-i2.0n 

c 



165 TAUT :TAU*COS CEPSLON-THtTAI 
l6f CALL cgT 

lfc? QTXr 1PRT 

16*» T2 = 0. 

lb« c*«**«««* 

1 7n 7C2 CCNT I-NU£ 

171 TTOT rTTOT»OT 

1 7,? t? = t:»ot 

173 0HRrCT/J6CC* 

1 7i £STrEST*OHP 

1 7* 
1 76 
1 77 
17* 
179 
lc? 

iai 

122 

\*7 
IS** 

ie* 

1*7 
15? 
1*9 
19~ 
191 

i«: 

1 9* 
1*5 
1 *6 
197 
19* 



-* r t 
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217 

V* 

23* 

2«C 
2M1 
2*2 
2*? 

2«m 

2*« 

2M6 
2N7 

2M« 

21C 

* 5 I 
25? 
25? 
2S« 
2S5 
25* 
25? 

2sa 

25« 

26 rt 
261 
26? 
263 
26* 
2b* 
266 
267 
26? 

»*? 

271 

27! 
27tt 
275 
27S 

2 7-5 
2S r 

:?i 

23? 

- Z* J 

2**» 

2«r. 

256 
2t? 

:■* 

2S r 

i gr» 

2C* 

2*7 
2*5 
?$u 
2* c 
2*>> 
2^7 
?•?« 

r*° 

3CC 
JC! 



3C r 
3 C 'j 

T * 7 

!:• 

3C^ 
31- 
311 

3ir 



2CM 
203 



ven 

C0NT 
COM 

DTXr 
CALL 

crx.r 

COUMX 

IF IV 

ifu 

CALL 

CONf 
IF i I 

'•SIT 

r»x ,5 

COUMX 

CXVEL 

CO T 

361 CONT 

IFIT 



,J,K )=(V2(I t J f X)4V2<l,J*l f *!l/2. 



700 



600 
607 



3 »«# 

3*3 



122 



612 



813 
61*» 

sis 

1C3 



CONT 

wan 

FOR* 
00 3 
00 3 
IFCW 

00 3 

irrK 

IF(K 

1 F C I 

COKT 
CONT 
IF |K 

CiLL 
C£TA. 

ruol 

CAMPL 

CONT 

CALL 

CRX,P 

COUMX 

cxvel 

WPIT 

FCRf 

oo a 

00 8 

wam 

1F<* 

CCNT 
00 8 

WPIT 
FCSP 
WPIT 
CONT 





IF(K 




ESQ 




WRIT 


3C8 


FORM 


13* 


COST 


s 


FCRM 


FORM 


« 


Fr.fi * 


5 


FCRM 




» • 


7 


FORM 


z: 


FORM 


33 


FC*" 


3» 


F v .;m 


267 


F*fi- 


7C3 


FOR** 


1C11 


r:^ 


10CP 


STOP 




ENO 



PRT-T? 

«CAL( IN 

•NXtMT, 

CU>*Y f CU 

C3.1.U 

SCTS t 1« 

ANCPRH 

NUE 

0*NC2) 

«6,703l 

TPftLOM 

,HX,MY, 

nUMY,DU 

KTE*P,A 

1C3C 
NUE 
.LT» 



,JN f KN t lM 1 JM f Ui.U2 J ll3,Vl,V2.VJ,W t 0H t ETAl f ET* t rTA3 # 

r*R,htHB 1 ELEV I UP,VB t NU,HV 1 *EX,hEY, 

MS ,OX,0Y t Oi»OT,FC0RtTAUX t TAUt »G »NCAS£ tNCLOK t TTOT) 

• 6C«2':0.) IGO:0 

,2.11 I.6E.5Q.I IG0=O 

M,JM f *(N»ue»Vef£TX|TT0T f R0SSaY»THeTA t T2J 

GO TO 861 

IN.JN.KN, lK i J*« t Ut»ll2tU3 t Vl,V2»V3tM.OH f £TAl ,£TX t CTA3 t 
MAfl t H,H9 t EL£V 1 U»iVb f MU»MW 1 HF,X,MrY t t2 1 
MS|PX»0V v 0S,Of»FC0R t TAUXtlAUV<G^NCAS£ a NBLOK*TT0T t 
Vff f ANG«THCTA ) 



TPRT> GO TO 702 



NUE 
<6,6C 
TUx, 
3 1=1 
3 J-l 

X(I,J 

<• Kri 
C5.II 
EC.KN 
.Nt.l 

NUE 

NUE 

TORE. 

STORE 

t,-Y, 

N t T I 

T,PHA 

NUE 

TP»LO 

t MX,M 

OUHY, 

KTEMP 

(6,U 

TUX, 

3 1 = 1 

3 J = l 

,J)tl 

Yl I ,J 

nue 
• r-i 

(6. SI 
T( lx, 
(6,26 
NUE 



7» EST, 

•EST.TA 

fTM 

tJM 

l.EO«0) 

,KN 

TBd.J 
I T3(I, 
)*ANO*l 



ECO) G 

T(IN,JN 

M4P,H,H 
E2K,TI0 

SE,OPHA 

« (IN.JN 

Y,MAR t H 

OUMS.CX 

,AVQ,iN 

2) 

•WATER 

• IN 
f J* 
^•••9 
I.NE.OJ 

,I M 

5) I, (W 
M=' ,IH 
7) 



TAIR,TEO t TAeN f TABS,TAUX,TAUY 
IRtTEfi f TAN t TAS # TAUX,IAUrc« t SF6.2,2Fl2.61 

GO TO 3«3 

,Kl=T2(i,j»irt«f IH8C I v J>*CT«ritJ»l«OS/2«l«0T?|I«Jl 

J«KI=1?ri|JfK-l) 

U,NE.KN1J TB»l,J,«) = (T2d, J»K>*T2«I,J # «-l))y2. 

TO 132 

♦ KN,lK,J*,in,U2,Vl f V2,W*,ETM,Tl,T2,TB. 

^.Ue.VP.MU.HV.^EX.^ET.NCAjr.veLOKiTTOTtOT, 

?3N, TIPElS.TIOEaS.TlQesStSTAGEtEST* 

SEtPEKlOO) 



SEtPEKlOO) 

,KS f lH.JM f UltU?,U3,V 1,V2,V 3,U,0M f £TAl f FTX f ETA3 t 
,HH t ELEV,U9,VB,HU f HV,MEX t MrY,iB, 
,0r,OS,OT,FC0R t TAUX f UUY v G t NCASK t NOL0K v TT0T t 
G.THCT* ) 

OE°THS # CHM 
WAT(I v JI=Ha<l,J>«CTA(I t J) 



♦ HF7.1/1 



TORE.E 
ILE 6 
(6.319 
M ix t • 
NUE 
T ( > 
T( l>, ' 
T|1T,« 
T< IX, -• 
TIMr s 
T C 1 X f • 
T( IX, • 
T( IX, • 
TUX.* 
Ti MM 
Tl 1K # * 
T(1X,« 



0.0) GO TO 13M 



) 

DOUBLE EOF PLACEa ON TAPE* I 



TEPS'.FlO.at'S* 
r C0R="|€»5.7 f * 



OY^MFIO.C^ CM 



OSSNFIO^, 



LSr' t I6» 

ncy^ . t mi 

Oxr »,F IG.O, • CM 

- - -- - SECM 

PEP SECM 

START"«f"ESTr* f F^.2, # CJiE NCr*,!*,' NSLOk;*,!*)) 

STAGEt AMP f H-$ PH4SE.1-W PMASE»PERf r SHlFT* t 6F7,?| 

COMPUTATIONS BEING STQPPEO BECAUSE OF INSTA6 ILTTT •! 
THETA 1 tFlQ.U* n£C'l 
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♦FICtff 
t 



1>.8AY60T FOR CREATEO ON * PEC 79 AT 09:N*s2« 

C ********* <«•**• ******* ••*•••*•***•**•*•••••**••»*****• *•••***•**«««*«< 

C PCIDS MARKER MATRICES t ElEV MATRIX, CCftfiHIf.ES "ATRICES M« » HU» MV « 



11 
1? 
J! 

10 

1* 
lft 
17 
1? 
1* 

:r 

21 

*, 4. 

2* 



:i 



."»6 

1 7 

? a 

JO 

kO 
« t 

* r 

M* 

«U 

4 C 
M', 

U T 

«n 
m« 
c r 
5 i 

5 T 
5* 
5 C 

* f 

5 ? 

S ° 

= t 

t ! 

t: 

6 C 

fcG 
6' 
h ' 

7" 

7 * 
7? 
7* 
7 '' 
7* 



7J> 



99 

71 

37 



ICC 

29 

52 

53 

5*» 

55 

1 

AZfa 

: 

*C2 
i U 1 

1C2 



109 



SllBR 
CETA, 
CLUMX 
CSTAG 

OIHE 
CV2CI 

tE tac 
cmtc I 
cued 

CMETC 

AG-l 

writ 

00 9 
READ 
yRlT 
CONT 
WRIT 
FORM 
00 3 
REAO 
WP1T 
CONT 
00 1 
OC 1 

ifcm 

ELEV 
ELEV 
H( 1, 
MCI. 

CONT 
W5JT 

cc 2 

WRIT 
CPNT 
WRIT 
FCR" 
OC 5 
RC40 
WRIT 
CONT 
WRIT 
F ?fl* 
00 5 
SEAO 
WRIT 
CONT 
*RIT 
FORH 
00 n 
R£AO 
-PIT 
CCNT 
WRIT 
F ORM 
00 <* 
RCiP 
WPM 
CONT 
00 1 
1 
If O 

M* I I 

CUT 
US 1 
CO 1 

If \» 

ml- < r 

CONT 
00 1 
00 1 

: f i* 

MV ( t 
COST 
-KIT 

f:sm 



CUTINE 9AVPOTtlN t JN t KN t T»SJM»Ul»lK f U3.Vl»V2 t v3iW\OW»ETAlt 
ETAJ.RX.RY^X.MYtMARfM.Mp.EUV.Ue.VBfMUiMVtMEXtrrY # 
iOUHY v QUMS 9 nx«OY f O&«OT f FCQR f TAUX 9 TAUY 9 G t NCASL«NPLOK,TTOT 9 

*SION Ul<lKtJN t KNI 4 U?CtN 1 JN*KNl,U3CTN,UN 9 KNI v VlfTNtJN.KNI. 



1«, JM| t FT«3f 1-rJHI .RXCIN, JN f KN I f R Y C IN , Jfl . H N 1 f MX ( I H f JN 1 

^«JNl|r*APf IK,JNI v HllN,jN| f HPtZh«jHl»ELCVC XNtJNlt 

►♦JM t KN) 9 VBIlM 9 JH 9 fcNl 9 HU{lN 9 4N) 9 HV (IN « JN) f »EX U»» J«> » 

IM ,JM> 

C.**9 

Et6, 51 1 

S 1=1,1* 

?» (MARCI, Jf »J:i 9 JNI 
Fffc*6GI It CMARCTt J) 9 J=1«JN1 
i.NUE 
H6,7n 

ITC !X,*eOTTOM ELEVATIONS OELOM MSU 9 FT*/) 
7 1 = 1, IN 

?, CELEVCI.JI ,J=1 ,JN) 
EC 6. 3d I,(fLEVtI 9 Ul, J=1»JN) 
TNUE 

CO i:i,IN 
C1 J=1,JN 

IRC I t J> •EO.CI GO TO ICC 
II, JI--ELEVC I« J) 
II V J) = 3S,«6*CIEV(X 9 JI 
•I=STIGE-ELEV(I 9 J) 

*ita*axI!C. 9 mci 9 ji ) 

1AUE 

Mb, 31 ) 

< 1=1. IN 

Ff6.3tl X 9 CHCI 9 JI 9 JS1 9 JN) 

i i iue 

Mo,5Z) 

ATI IK, '"EX "ATRTX* > 

2 1=1.1* 

2, f rtx 1 1 f J I .J=1.JM| 
Hb.SO) l» CMEXCI, J) 9 jri 9 jM) 
INUE 
EC6.5M 

/T( IV, **EY **ATR!X« 1 
! 1=1. 1? 

2, CrEYiI 9 j» ,J=1,J*» 
[(6.6C) I, CH£Y(T f J I t J = l,JMl 
1NUC 

Ue ♦ n 

M( ix f •** MATRIX* I 

;* i = lin 

2 , <^XlX,Jt f J=l,JKl 

£(fa.bO i. ipxr x«j) 9 j-i 9 jm) 

INUE 

f (6,3) 

M ( IX, ,M Y MATRI*»| 

C2 1=1, IM 

», («YC 1,J» ,j:I,jki 
t (o ,bC) It CMYCI.Jl .JH.JNJ 
I'.UE 

CI I-l.TK 
CI jri.^K 
Excitji.EO.r) 5 n TO 101 

,J|r(M(I 9 J)*H(I*l 9 j)*HCT*l v J*l)«H(I 9 J«S)S/M« 

iNCfE 

c: r = i f p: 

c? J = l ,JH 

JlltJI .tCC? 00 TO 1C2 

,JI = CMI1 «J)*K(I «J* II J/2. 

l^uE 

C3 I r LL- 
CS JM ,JN 

M t.JI .EO.CI r, To 103 
,JI ZIHCI ,J J ♦H(I*1 9 Jl 1/2. 

/TUX, *HE MATRIX'S 
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uc 
111 

112 

in 

2 

:c 

31 

*0 
M 

u 

61 

250 



00 1 

wan 

WP1T 

f OftH 
00 1 
WSIT 

yRlT 
roa* 

00 1 

W9IT 

fork 

roan 
roi*H 
roan 

FC»P 
FORf* 
FCRM 

FORH 

forh 

RETU 

£N0 



10 1 = 1 

E(6. 11 
*T( IE. 
12 1 = 1 
M6,3t 

Mlfx, 

1*» 1 = 1 

EU.30 
M < ) 
lTIJX t «I=».l<i.l5Ka.ll 
MtfX t «WM£« nCPTMS, CK'/» 

MUX, 

MMJ, 

MUX. 
MCI 1 



» X,tHB(X.J>.J=WU*) 

L ^MU HATffTXM 
• IN 
) Xf CMU(X«Jl«JSltJM1 

•MV MTftltM 

.1* 

1 I«(MVCX,J|.J=l,JNl 



~X=*.!«». 1<E6.2> 
•POTtOM SLOPtS, CM/CH*/I 
•MAR HATRTX»/| 
•l=*tl«t 10X V 1S141 
•STAGE=*.F10.2) 
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•'lowm.pAYfu fop ;»E*Teo on « nee 79 at 09*«i5::m 

2 C INITIALIZES MOST OT THE OCPCNO'NT VARIABLES 



T 



<* 5LSBCUTINC *AYlWXIIN,Jfc t KN,XM f jM t ui,U2,UJ»Vl,Virrf3tM»0»»CT*! f 

t CCOMJ^UJJY.OUhStrx.DY.OS^T.rcOR.TAUX.fAUY.GtNCASetkeLOKtTTOT, 

* CTl f T<»T3fTC f T8»TI^IT) 

t*£f$ICN Ut Clh«JN t NKl fU2eiN f jr.KNl,USCIM,wN f KNI t ViaN f JN,KN| f 
CV:iXSJN t NNI t Vi(IN t Jh v nni f MriH 9 jHtKNI ( ONfXHtJP* 9 IINI t FTAlClH t JNI t 



CETAf l-,JH) f CTA3<I%JHI t»m*fJfc.Kfc| t ftY<IN,jN t KN) t Kj(<lH f JN> t 
rHYIJK, jM.WlPtlN.JNI.HUN^NI.HedHfJHI.eLlVUN.JNI, 



U2tf ,J«K >:J, 

u3»t ;j v m >ro. 





w 

1? CWCYM^.JMI t T?U* t JP f KN) 

1JJ OI*E*SiON Tl(lM t J»l t «NltT2II» t JM»K»iiJ # I*llM f j« t KNI,TC(lHtJH t HNI 

tfa QC*!Cn fri.tw 

p 00 ICO Jri.JN 

!• 00 ICO K=1,KN 

2: VUI ,J»K)ZC. 

2! V2II ,J,KI-t. 

?• V?II ,J f * ire. 

i«i Rx<I ,J,k )ra. 

24 »Y(I ,J f K 1=0. 
2* ICC CONTINUE 

1* 00 ?c: in .1* 

:« 00 2i? jri.jK 

5: crAi ft f ji-o. 

31 CTA( T.JISO. 

3? ETA3 (t,J»r:. 

3' 00 3C0 « = 1 ,KN 

3* UfiCl •JtNi::. 

3 s V&ll ,J f K )=C. 

3S u Ht »,m irG. 

3T 0**ll ,J,K »r?. 

'■ Till ,j,ki::tmt 

5« T2C1 ,J»N ISTINIT 

** T3(I ,J f KITTINIT 

*»l TCtI ,J»K >=a. 

*C TMI ,J»K| = TIKIT 
*« ' JCO CONT J.HUC 
«** iOO CCNT ]SUC 

»* ASjr 13.«»9 

*6 CO **C0 1 = 1 tH 

M7 oc **CC jri »JH 



*•« £ Til l( f JJ = 19J 

ST £ TA( I t J)rA«J 

SI £TA3 (t,J >I1PJ 

5-' CO 5 CO k;1 ,hn 

J- T5(I ,J.Kli*3J 

5« b f Ij^KI-ABJ 

; c SCw COSTIMJC 

;J «Cu CTNTINUE 

* p £S0 



IF (?£*<!, J> .*E.01 GC TO *G0 
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♦riOWtn.fltAOT TO* CREATED ON 7 Q?C >1* AT lC:C«:Ct 
2 € ftC'OS X> HOURLY RESULT lURCO IN TAPE FOR CONTINUE RUN Off 'LOTTING 



• SUB!* CtJTTNC PEAITI IN, JN 9 KS 9 I«» JM 9 U1 9 U2 f Vt. V2 9 W« ETA UJ!. 22*76, 

6 CTIOE lN»n0C2N t TI0eiN t TlOClStTtO€2S,TI0C3StSTAGEteSTt 

7 CAMPL 1T,PMASE f rrHASt*PtflIOOI 

1 CW(I* .JH,KN) 9 ETAi( !* f JH| ,CTAf !**,JH| f -X (IN, JN * 9 *VC IN, JNI ,MA* tXN, JN1 , 

10 CHUN ,JNI ,HMIN,JHI ,Utf Xf*«JH,KNt ,tft (Irt,Jtt,KN| 9 

U CHUlI K 9 JNI t HV(TN*jri) • *CXU« 9 J»> v *£Y<tN,J*)«TP< TP.JP^XNl 

1? ... .... -._ .... 

IT 

!• 

1« REAO C7 f £NOril NSLOK 

1* RE»0 (7) ( CIU1(X«J«K1 >K:i,KN) t j:i, JN»,I=I,IN! 



(.null * 9 J'l I |KV1 [Ntjrii t"t*i l"iJni ,nti \ l.T 9 *J-» • i r» l*.wr,^i 
OtMpSION Tl (lK,JN ( irNI t TXPClNf X«l t T!CL2N|lf«l 9 TTOe3tllTN I 
OI*E*SION T2(T*,J» t KN) ♦U6ClSU u l .UDcIStlf »,U0t Ji<l*) 
X CCNT INUE 

REAO C7,ENOril NSLOK 
If REAO (7) ( (IUl(l,J*K1,K:i t KN) t j:i, JN»,I=I,1 

IT CI((U;u f J ( Kl t K:i,KM,Jil,JNI,I:l,IM ( 

1» CH(VlU 9 J,MI 9 l<:it«Nl 9 jrl,jNt 9 X = t 9 XN), 

1« Clllv;(ItJ f Kl 9 K:x,KNI,j:] v JN) v X=l«XNIf 

2 2 Ct I MX (i,J),J=l, JN> ,1:1, IN) . 

2? CKMY (f,J) , J=l, JN) ,1-UXNI, 

2? C( (HU II ,J) ,Jtl,UN) # Iri 9 INI , 

2«t CMHV IX,J| , jrUJNI V I:;,XN» 9 

; c C< IHH.J) tJ=l tJNlt irXtXM 

It. REAO 17) (11 Ml, J, A J ,Kn,KNI 9 JZ1 9 JHltX=I#XNli 

2 7 C((IUPrX v J«K| 9 K=l v KN| f J=l t JHI t I=t f XP) v 

Z* Clt(V<?(I,J,KI,K = l,KNI,j:i t jH» f iri,lM, 

2« ril(TlUtJtKI»N = ItKN} 9 J=ltJHt v X=l 9 irl 9 

2 * CCUT;fI 9 J*K|»K:i««Nl t J = I, % jMl,I = liXM» 

31 C(<<TE(I*J 9 K>,K=l 9 f<N},J=l.JHJ«X = l 9 iK>< 

32 -- --- - -- ...--- 



c< i£tmi 9 ji 9 jfi;jH) ,|=i 9 !ni; 



VI l t i C * 1 « W 9 * I ,n-i 9 "nf,*j-4 ijn , 

C( (CT ftiit ji *jri 9 j*4i l i:i l !*<i 1 
MI, J) 9 jri t jM) 9 f 
Ml, J) ,J = l ,J«) ,1 

HI,J) ,J=l,J*»,1 

if 9 JI 9 J=U j^I ,i: 

r lN( I I ,1-1 ,1- > , 
E?N« II ,1:1, JM) t 

[ 3nii i ,r = i,iH) , 



14 C( IKC XI, J) ,J = J 9 J«| ,t = l 9 I*>, 

: r . ci («c x i, j i ,j = i,j*»,:n f i*i, 

!f CI (MR II.J), j=l, jk> .irXtXHI, 

! * CMWrixin, ' 

3" ci not: - 

jo come: 

«n ci riocisix j ♦ r = i , i« » » 

«*? CI TIO€!S( XI. 1=1 f I«l 

M* KEAP (7) TT0T,KCASC*0r t STAGCtCST 9 AHPtIT 9 PH4SC 9 0PHASC t 

m* cPE^tcn 

<•« »CAC (7,CfiO = 50CI AAA 

Mfc W3XTJ(6.2I 

** 2 ^»h;Tiix.'vo EOF AT £N0 OF OATAM 

*e NSLOr = *iro 

*« GC TC 10C0 

sr SCC CCNT iMjC 

5t UflITCfb.li;> TTOT.NrASE.^et^X 

«? 112 FOR^;T( fx, T r>ATA P?AP FfcOH T»PC, TTOT= t 9 riC.0 9 

5^ C CUE N*fr = %15,* etOK NBPS^Ib) 

5* ICOO PETUfN 

55 ENO 
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•ruowt n*i«»e*o sy* c*caud on * oec t* at r«:«7 ; s* 



1 c < 



H SUM CUT INC rPPtlOU*»J>UKN t Tl f T?> 

* OtHt»510<4 TMlH|J*t«Nl 1 TiM* t J* 1 «N> 

6 C*«**«0*T* «tAO r*ON OHt Of !h€ 4F *H»C *U lCAf A 8* Iff *•••*•*•*•*•••** 

7 00 5!0 j:l.Jf 

* SSO fiCAO J, fTltX**JtlltX:i»X*t 
« 2 fOflHITI I 

OC }?5 un t g* 



00 *!5 1=1.1* 

00 5«J Kr?,** 

55J TUT *J«K t = Tl f T ,J, II 

o: s«* k = i,»(n 



* 555 CONf INUC 

7 RCTUfiN 

9 E*0 
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•rtOWUKCCTCf'P SY> CffCArfO ON 12 OCC 7* *T 2Cs3»iS»» 

* c ********* «♦••*••••••••••••••••••••••••♦♦•••••••••••••••••••••••••••••• 

J C COMUTC! tOUlLl6«IU" T ANO SUflfACC HCAT CXCM«KfiC Cat? MC ttN ^ s J ##t# 

c , TAIR l = TAlk**./$.M2. 

» TSU(MUTSu9r«*./J*«T2. 

* WlNOUWtN0/*«*7 

* S»AO U*PI0*1**C 

i fo£whu«.5S»,i.**T;tRii«M,-Mu;roi , 

Ift lCtwJ:(<2.S*.tJfl7*T*tRll*U**HUHl6ll*«3 



n #f !• :.25*-.2C{5*TrUH#*a0023*»TrilH 

{; FCNU:70»*«7»WtKPl*»J 

it $t(il *. 7 ♦liCT *♦♦:•» *FCNIJ 

li FCfUFN 

2C CNO 



••i 
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• TLOKULHTA FC* C?t*Ttn ON • OCC 79 AT rilMtl* 

? C CQ**UU? tT» ISO 0»*tSA Of AOVANCt T I Ht ♦ g$CS St"»SCM*i UUlC f Ct TNTG 

? C« 

SI 



« 

11 
1? 
IT 
I* 
1« 
1* 
17 
t« 
1* 
2C 

ii 

ii 

2« 
3? 
31 
IT 

J* 



3? 
3* 

MP 



CM 

C 

cc 

C" 



100 



£00 



900 



UlffCUMNt *CT*U*t JWf*N« I* ♦ J* .U* ♦ V2 » 0M f C Til f C T* f t U3 ♦ 

B.Mi.HV.ffFXf-et.ci.o:* tOTtO* tOt i 

lf*€^;OK U2(fN*JN«KNI,V3 lT'l 9 Jfe 9 ftfil 9 0ft(lH,J»SKNI «CTM<X«tJP>i 
TA< ]« 9 JMI 9 CTA3U««Jni ,H8U«tJft? 9 MU<lfi«JN) 9 HVUK f JNt 9 HCXUM t J*l 9 

ItfftSlOM Ptftl 

!0=*.*C 
75=*. 7S 



«! 



• c- 

9 

*«* 

L=M 

r«* 

2 = H 

M 

H 
I 

1 

hux 

MVf 
(HI 
0*T 

u*r 
CT: 
TA3 

ru 

TA3 

CTr 

CM 
H = H 
* (I 

•mi 

: ( { 

tTU 
NO 



(••0 
►•I 

ca i 

CO J 

Mil 
L<I« 
IH? 
Ml, 

mi 

Wit 

mi 

C3 « 
:fOR 

:(0? 
:pmu 

iur 

Ml 

•?*• 

IT 9 J 
N§< ! 
► »*i 

• T.J 

ICT* 

l*\Jt 
HI, 
•J. 2 

«J«1 

isuc 

FN 



II 



T 

= ltTr 

ri ,j* 

•JtttC.QI Cr TO 930 

l 9 Jt*f£fMl*l«J>*CU<ItJn//. 

•JI.CC.2) G»:HU<lM 9 JI*£T*ll, Jl 

JI*<FTA<:*l,Jl*CTA(f # JM/?. 

• Jt.19.1) Cl = MUI l.JMUIIhJI 
J*ll*l£fMX,JMI«CTM f f J I I / 2 • 

iji.tr.?> o?=Mva.j*i )*CTi(t # ji 

JI*(f7A(l 9 J*i)*CTHT 9 jlJ/2. 
tJl.fQ.ll 01±MV(X 9 J>*CTAU,J) 

*u* t !♦!• j»*I-0l*u2<I*J»«* )/0r 
•V-U.J*l 9 KI-0l*V2<I*J f ^n/0T 

X.CHVT 

<tl* 

I=£U1 (T,Jt«AlC»CCT 

• J 1»CTA3< T.J) 
f . C • J GO tj TOO 
lr(lC.*»C -6>)-H§< I, Jl 
3(T*J)*CTil<I 9 Jti/AeC 

JMCTAI T, Jl 

)r(-5?5«CCT*05*«r(ll*r? ?)!/•> /AM 

» = (-r5C«o?T.o$*4rtiwR*K:t*ft3)/*n/»H 

iri»S75*0tt*0S««Ftl )/*♦*! 2 J*f (2 I •*(«!/« I) /am 



tf8 



LOW* 1J.9NRTIA FO 

c cc m pute 
c *•««**•*« 

SL6R 
CCTA t 
COUPX 

DIME 
CV2II 
TETM 
CHYlx 
CUStI 
CMETC 

ABC = 

ACE = 

C****«CCMP 
00 1 
00 1 
IF(M 
Ah:h 
AH = A 
OETr 
00 1 
02=H 
OUH 
U°AP 
UKAR 
OHUU 

in* 

£2 = 1 

EUf 
02-H 
OUH 
OUA 

o:=a 

UPAP 

ubap 

VOAR 
V3AR 
CHUV 
GO T 
1G* CONT 
C 2= f 

02=H 

02 = A 
U3AP 

ve*R 

OHUV 

106 CONT 

IF1K 

AUU 

AUU 

ouc« 

DUSr 

GO T 

1C7 A3rO 

a2=u 

A UU 

OtCK 

Otir 

ica SIG: 

R)MI 

in ccmt 
ico cost 



C*****CCMP 

00 2 
00 2 

Ah:h 
amt 

oet= 

00 2 
02rH 
OUH 
V»AP 

VPAP . 
OHVVt 



CPEATEO £N * DEC 79 AT 09:*8:<»Q 



NONLlfiEA* TEP«S RX/RY AT INTERIOR HALF-GRID U-/V- POINTS 



I1T11E aNPTTAClN.JN.KN.IH.JK.UNUZfUJtVl.VZ,* JiM.OP^ETAl, 

TAj^t.RY^ItMTtflAP.M.HB.ELEV.US.Ve^UtHVtMEX.MeY, 

OUMTtOUPStOKffOrtOStOTtFCORtTAOXi tAUY ,G »NCAS£ »N«10K ,TTOT r CB) 

SION Ul<lN t JN t KN! ,U2 t I'i . JN f KN I t U 3< IN . JN, KN » , V 1 CIN t JN t KN I , 

,JN,KNI,Vi(IN f JN f HNK«rlM,J» t KM f OM«lH f J« f KNI,rfAHIW f JMi f 

H »JNl v ETA2f X«»Jftl ,RX<I* f J*,MJ> ,RY(lN t JN,KN> ,HX(IN,JN), 

,jN| t *AP(lN.JN! f H(TN 9 JNI ,HR t I N v JH I ,£IE V < TN» JN > , 

• JH^NltVBtl^tJHfKNltHUdN^NltHVllHtUNltnEAIlMtJHlt 

"»J*I 

-I 

8»0T 

TE5*RX 6 AT TNTCP10R U-POINTS***** 

n 1=2,1* 

o J=1#JP 

<I,J>.N£.3l GO TO 103 

Cl!jlMETAtl,JMETA(I-l f J|l/?. 

AX1(AH»*CEI 

rTA:(I,J)-ETAiri t jl*ETA3<l-l»JI-ETAltI-l t Jll/A8Cy2* 

KT1.K2 

{I-T t Jl«ETA(t-l f J| 

= tU2ITtJ«H»*U2U»l»JtKll/2. . 

= CU2CI f J t Kl«U2tI-l»JtKM/2. 

T(ar»cuB«R?*«2i-oi*t(jaARi**2i»/ox 

IXvJl.E9.il RO TO 104 

TACI-|.J#ll«ETAtX »J*n»ETAfI-l»JJ*ETAlI,jn/»l. 

TA»I,J>«£TA(I t J-l J*eTAII-l»J-ll*£TA(I-l t jn/* # 

r f j«n«e2 

I,JI«*1 

AXKCltUl 
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C QIMi^iOM 9 U2?tN«V^* ffN >* v 2tXM 9 JI« t NMI t U»CfN t U» 9 «Nt»V«fZ«f t Jfl t lltl| 9 

DlPTftSXON HtXlXH ( Ulll # HCf <X«tJltt»OTZtlH»JHI 
ABCrf«»OT 

K2:k»-*1 



uc 

100 



RR.RYI 
00 |H X-ltlrt 
00 ICC J^itJ* 

OCTr KTA3fI,J)-f IAlfI f j| l/ARC 
Q*SMHX*ltJt*l£TAU*lt JM£TA<I 9 JII/2. 
IM>*£Y(l,J1.C0.2l 0»:MU(lM.JMCTAtf t JJ 
QLSMirXiJlMCTAf f-l,J)*CTM t,J) 1/2, 
I? C»*f wt I , J). £0.1 5 OL-HUlitJt*CTAUtJI 
02tMUl.JM)MtTAtI f J«t)*eTMt f jii/:. 

WfBlMlSlEJ'iM'fiHSHiMirtM:-! 

DT2< l,J|C-HK-(TCO*T?U t J t ll) 
OC I JC KSl 9 tfZ 

TR = T<rl*JfM 

JM'ETt UJKSC.2I Tft=«T2tI,J,«MT2<I«l f J,M>/2, 

lKJfJfljJUfce.il TL=n2II,J f «I*T2<I*UJ*H)l/2* 

ONUT>:(0ff«U««rR*OL*Ut*1LI/0« 

Vfl=(t2Cl*J*l»«)4V?UtJ*ltK*m/?« 

vi:<t?irtj*Ki*vr<x 9 jtK #111/2. 

T R Z T I ( 1 , J t K I 

lF4»CY<ftJl.»C.2) TR:iT2<l,J,M*T2(XtJ«l*KM/2« 

TlsT2f X»Jt«l 

IF (HeTli.Jl.MC.il TLtlTJCIiJtHMTZIXfJ-l.Kll/?. 

0MVT>;J02»V*»TR-Ol**L*lLI/O* 

0*R=C*<ttJ»K> 

M L:C<X 9 J 9 R*1I 

IF IN .*C. II TR:(T2(X»JtH»*T2tI t J.K-lll/2« 

inn .eu.n TR=T2<r»JtK >*(os/2.i*ot:i uji 

If IN ,WC#K2> Us I T2(tiJ*K>*t2<Xt<i 9 X*n 1/2. 

0T0^5t|0HP«TR-O^t*ft»/0S 

OTS= |TR«Tll/C5 

SX6s*'LeAT|N-ll*0S-OS/S« 

I C 1 1 t J,«IS0HUTX*OHVTT*AM*0TOHS*CCT*OTS*«l«*SXGl 

CONf !\UC 

CONT XVUC 

RCTU«N 

tNO 
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•rto«f iktcohm for cacateo cn i am so at i*t2't3£ 

2 C CQPfUU! T3 BY CTCS*OUFORT-F»AN*ei SCHEME TO INTC6KATC T-EON 

T C«»*«*« •••<***•••••*•»••*****••**«•••••••**••••••••••••••••••*•••*••«*• 

«i CMi.Te.Tl.lr.TS.KCXtHCTttTA.iBlTAlS.tTAltnAitOT^.jV.IlMlTASNjTAMl 

h ClMEMION ?CtXKtJ*tftfcl t Tlt t^,JM,HN» f T2t tM t J-,KNI,t3U» f JH t HNJ f 

* OIMCMtON OTZfl^ c JM| 9 HCTtilt v JNl 

• K2:*ft*i 

1* AecrCT 

ii oe uc izi f tp 



00 ICO J=1»JP 
lFi»fvtI t ji.rc,«i tc 



I* ANSMffXtJ»*CT*U 9 jf 

1< AMirM(I,JM€TAMI,J) 

if AAcsev/AH/tos*«n 

1! ACC=IM»AH 

1« OEFr ] •/!!.♦ AeC«ev/ANi/ AH/ I0S*«2 II 

:0 0E6 = L/<U*ABC*9V/2.SAH3/AM/C0$**2II 

i\ If (rtx(irjf.fo*3l02TX = a2(I*l,J,K)»T2U-l t J,N ) »2**T2f X 9 J V K)) 

:» lF<PFY(I t J>«E0.1l?2TX;tT2f X«lfJ 9 K>-T2<X 9 J 9 fl)>/(0X*«?| 

a IF<MY<i,J>.?0.3)02TY:<r2fI 9 J«l 9 «t*T2(l 9 J~l 9 KI*2.*r2fr 9 J 9 K|l 

27 C/ tCt««^l 

2« IF t«EY< I, Jl iEO.tl02TTriT2lI,J*l f lO«TAbS-2.«T2CX.J f K 1 1/ <OY««2 I 

2' lr(K!V(ItJt.CQ«2in2TY5(T2ll»J«UKl«rA8N«2.*TltItJ 9 KU/lOY«*2) 

3? ITCH .CQ,1I G<3 10 5w 

3t IF l« .C0.K2I GO TO SI 
1*51!'?" 



12 SAA1ST2I It JtR*ll-rilI,j f iil«r2(X«JtK a ll 

33 BAflrt-TCII.J.It )*AiC*8ARl*ACC*<02T»*Q2TYI 

3* Til! ,J l K»:CCf*<*H«T2tI f J f «)*A8C-*Afi2l/AH3 

!< . 6C TC 230 

35 SO CONTINUE 

37 SARI :T2f 1, J,*«l l-M U 9 J.« l/?.«OS*QT?C X 9 Jl 

!« T 311 .J,** ):OCG«(*H*T:U •WtK)*AflC*eAB2j/AH3 

*? GO TC 200 

«l 51 CONT 1N0C 

»? SA9l:-TUX« J,KJ/2.*T2f X,J,K-1' 

<? flAP2:»TC(I 1 J 1 Kl*AAC»BA«i»ACC»(02T)(*P2TYI 

M fStl tJ,M l:6cd*(AH*T?(I t J f KJ*A8CMAR2l/AH3 

*« 2C3 CONTINUE 

46 100 CONTINUE 

#7 »eTU*N 

** t*o 
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•fUOWUl.GIVCNT $Y> CUATtO ON ♦ OCC 7f AT lQ;22fU 
t C*»**«JPtC irt fC"»tff|tU*C AT 0MCMM6C **0 ftTVtff u "t'Q*"*** •••••••*••• 

C SU{M10N VltlNtVt 9 KN»tV2ttN t JI««MNI,VlftN v JII*Kflt« 
CtUiUN 9 JI$ 9 Mf«l 9 U2iXlI«jii*l(llltU)ilN,4litllNl 



200 



00 ?C3 HSttKN 



*ctu** 

CNO 
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«novui«tfCAL ro# cu*tto o*« • nee ?? at iux«mi 



i V 



! 



CALCULAfCS »CTUU fC'TlOl VCUOCITT V 



* SI 

: m 

T 01 

• CV 



11 CUia^*J»tK^l«VtClM«JII 9 liMltMUCtK»«lftl 9 MVltl«*JMI*l*CX(lNfJIII« 

1? ClJltf )*.J*»I 

i* «?=«*-! 

1* 00 tCO ISltttt 

U tniitifi|ji*co«oi oo ?o too 

is tasflijJiiitirt 3i ,,tJ,, ' , * /af 

i« xmhems ♦Jf.eo.n to to is 

sr: XMftva *ji.cc*ii go to i 

h fri»ly<i,Ji.ce.?i on to 2 

2 5 * 0?5sieYaf!«ltJI-fTMI-l 9 jii/ouni 

2« 0AMi:tH0 4t«lfJMCTAU#l t UI*M0a«it«lt-CTAtI*l f Jll/OUNl 

JC So T( 3 

U ilscni!-UjMi.!*iffIii f ji-CT*n-itJM 

3? Cl = tCTACT-l,Jl»CTi;t,JM/2. 

3* 0AHi:iKUII#JI*C2««UU»ltJt-'XI/0I 

3« 00 TC J 

1? » tfWfcfl.JI.ft.il <C K 11 

♦• irtncttitjueo.it •<> ro * 

3« iM*c*n.jKCQ.:» go ro & 

5? l * S?r= <?¥*«!, J* n-Ct*q t j-in/OU*r 

*: o*My:iNdaIj*n*tTAij f j*ii-«ca*J*-CTAii,jj!/ou«t 

*»T co tc * 

ft* 4i OCT: ICWU ,j*U-CT*n.JM/OV 

ft« r 2S f f TM I , J* 1 1 *t Tt < TiJIJ/2» 

** tin uu.3*if-wsMttA<f»j«n-CTMi # jii 

•I 60 TC ft 

SS C«sC1«U*J*tl*l*S«c'TAlit Jt-CTAfXttl«lll 

SI t|TCCTMI.J*n»CT*«T,JM/2. 

St CC TC <i 
10MTI4UC 



5? " 6C ••« «rt»K2 



T CQ»Ht l*0CT*U*<X*«Jf*l*0€X«Vt4ltatftt*0Ct 

«*% C3*t nut 

• tec CCNU*UC 



*: »ctu*h 

41 €K0 
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vflQHClt.AfcCM '*• C»t*TtO OK J OCC 7« *T JOtUiM ^^ 

¥*¥• Ittlttt 
XUSC»fttt*«*«V«*7t 

II BfoSHHC 



7 

to 



1* StSCUtll 
H 



ZCC I 

340 < 



TMStm.t.ll 

viKtlli:S:!i 



11 

j iutti&f" »«•»■• 

27 CALL 7{HtUV.*(DI 

31 f»j:l2tltl2tll 

J| u=ut ntitttit 



3? *t¥t IMfttll 



>H0< 



14:5C'T<U*«2*V+*l> 

:« citL zzitu.v.;toi 

3f |«iso«Um*2COI 

3? in?«.lT.Q.» 2%r2«*! 

3* I* 5 ! 1 !*!* 1 

« ! * wPXT|li,ll*?T0Tt$»*JfU/l,twtfS2tl7»22#l«2»Slt«lf2 3,TWj t 

J) i C ?o4JinJi!r7.o f 2r*.i f r».o,ar4 # i t r$ # CtV*«UM.c^r4«i« 

%« tf.O 
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•rioitt 

1 

5 

< 

IV 
11 
X? 

IT 
!• 



*7 

: 

21 

22 

!? 

26 
27 
2* 

29 
30 
31 
22 

t T 

!• 

25 

36 
J7 
3» 

1Q 

ic 

*! 
*! 

46 

*7 

•• •* 

MO 
c r 

Si 

s? 

5? 
59 

I? 

57 
< « 

59 

6" 
61 
62 
6J 
ft* 
65 
66 
67 
6" 
b« 

7t 
T2 
7! 
76 
7* 
7A 
77 
7* 



11.TPRL0K FOC CSEATEO ON « OCC 79 IT 11*31:18 

C**?5lMTS CUT MOU«LY**CSULTS*OF HORIZONTAL RESULTANT VCLt W-VCL» ANO T 
t AT * LEVELS* »NO THC SURFACE ELEVATION ETA 

SUSP CUT I C| TPRLCH C IK t JK 1 RN,lM t J« 1 U t.U2 t U3, ¥!♦¥*,¥ J^y.ftN- 
CETA 9 ETA3 t 5* 9 RY 9 *K 9 Ht 9 MAR ' "" " *" " *" *~ 



COUrXtOU^T.OUfStCXtOt.O^ 
CKVEL ,KTE«P.VR,ANG 9 THETA1 



KfKNtlMtJv.Ul.U^tUStVltVl^JtU.AN.ETAt, 
AR.H^R.rLiv.Ue.VttMU^V.H^Xt^EY,!}. m 
S*0T t rZeA v TAUX 9 tAUY 9 G»NCA$S»llBL0l( ( TfOT t 



OXflCtSION Ul(lN*JN 9 Kh>,U2(XN 9 JN*KNltUi<IN 1 J* 9 KN) 9 Vl(T^ t .JN 9 KN» f 
CV2tI>,JN f KNl 1 V3aN t JN f KNI t tfCtW,JM t »Ni,0ff<lH.J»',KNI,ETAl(ir< t JHJt 
CCUI lM t jtt| t ETA3(I M f JNl f RXlXN. JN,KN1 t RY(IN 9 JN 9 KNt 9 NXllN v JNI t 
<*Htm.jN>,"AP( IN, JNttHIIN v J*iitH0fTK«JH} 9 ELEV(IN 9 JNI 9 



CU?(I^JH f Whl t VB(IH,jH,KNl ,HU(IN 9 JN»,M%tIN f aMt#*CX<IM t JH» f 
C*EY( l*.Jtt| 
OtPOSlO* VffUHttANGiJHttUCXNtJHtHKl 

NPROrrQ 

00:51.3 

' K2=KH-l 

IFtKVrL«eC.O) 60 TO 88 

IFIKLV.E9.ai GO TO 200 

00 2 K:i,K? 

WRITH6 f 5CCI K 

00 3 Itltl* 

W PIT E 16, 50 1 1 If (U« Cl 9 J f K> 9 jrl t jK) 
3 URITC(6.5C21 t VMXt J»K > t J=l t J*> 
2 CONTINUE 
200 CONTINUE 

00 1J K=l t K2 

U*ITE(t»7CCl K 

00 12 irl.lH 

00 1 J=ltJP 

VRCJUl0.»»o 
1 AkU( *)=1Q.*«9 

00 l! J=l,J* 

IFIHEXU ,J» .EO.O) GO TO 13 

UrUB II 9 J f Kt 

v=ve (I.j.kj 

Vflf J1=SCRTtU*«2*V*«?l 

CALL ?ZHUtV«2E0l 

AUCI «>r9C*ITHETA«2En» 

1F(A*6<J>.LT.0.1 ANGJJ>-360.*AhG<J) 
13 CONTINUE 

WRlTEf6 v 701) 1,(V9(JJ ,J21 t J«| 

WRIT 1(0,7021 (iNG(J> 9 J = l»JHl 
12 CONTINUE 
H CCNT 1NUE 

00 * RSl.R? 

WR1TC(6«SC3I K 

00 5 I=1»I* 

5 wPITff6,50«l I 9 C«fI 9 J 9 K) 9 JSl 9 JKl 
* CCKTINUE 

WRlTf (6,5051 
00 b I:l.I - 

6 UtllTF(6 1 SCSI I 9 (ETA (I. J| ,jrl 9 JHl 
IF(RF30F.f).01 GO TO ICC 
fcaiTCf6 9 «C6l 

CO 1C JritJH 
««ITC(b,6C^t J 

co ic K=ijnr 

W9ltE(S.5C71 
10 CONTINUE 

WfilTf I6,5Q9> • 

00 2C 1= I. TH m 

00 2C K?t,K2 

WRIT C(b, 507 » K 9 (Vftt T 9 U 9 K> 9 U?lt«J»> 
2C CCNT TNUF 
ICC COM Ihu£ 
66 CCNTlNuS 

IF(K KNP.EO.G! GO TO 89 
00 3 1 nri,KN 



K 9 (US(T t JtK1 9 l=ltl^) 



dRiri(b 9 55?l K 
S5C FCB>Mux 9 »TE«PfRATlWE$f 

CO 3 7 l:i,! M 



OEG r t 



RS%I«I 



IFIK.Ec.l) UP!T£l6 9 4$2) I . I T2 C 1 f J 9 « I f J = ! , JP I 
IFtK.r.T.l) WflITt(6 ? 55ll I 9 lT2(I»J f M 9 U=l 9 J»M 
32 CONTINUE 
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7« 552 FCRHTUX, *I= # tT9t« TC«^=% 1**7.2// 1 

I? $9 CCMTIWUC 

45 S03 FCR*irux 9 *U 9 V VELOCITIES, CH/SfC* KS* V X*/J 

«■ sci FCR>*Tu<;Tr«.u.« us'.Isn.n 

af 902 F04«JT<iX,ax 9 *%:*,15F?*l/> 



K:'«X<i/I 



5C5 rc«tt#riu t «jo»f*ce cuvaiioss. rr** c**/i 
Sot fo»**Tiix;*f?* 9 i<i f ' f TA:*,i*t7«i//i 

lin FCR«iTux;«ve»nclL prcMlc! or u velocities*) 



aa 
a* 

«c 

91 *—.•„.. , , fc -..*— ^-« . M * « - ^ fc w^.. ftfc * , 

91 509 FQft»jTUX, 9 VE*TlC4l PROFILES OF V VELOCITIES*) 

95 *CC FORHITIMM 

9* 700 F0RHlTClX t *R£SULT4M VCLOCITICS ANO OIRS* K:*<19/1 

9* 7C1 FCR*iTUX 9 *If f I* 9 # VRCSS* 9 1SFT.1 ) 

9J 702 FORWlTlU t 12X»l5rKo/l 

ICO CNO 
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**L0ttM>.STO»CT FOR CREATED ON • OCC 79 AT XitZtxZk 

"? C *S TORES* ►CURL* RESUIT ONTO TAPE P OR*LA ?ER*RUN OR PLOTTINfi 




i 

« 

1C CHUN,JN» f H«il*,JM»,UeUK.JH,KNi,VbUM,J* t KNl t 

1! DIVISION T2IIM f J«.KI«l t T10ElSf:^>tlIOE2SaHl t T10e5$a»l 

1« *9L0*rNBLOK«l 

1? WRIT? Sll mUHI t J*K> f K = t t KNI t J=l,JN) f I=l,IN! t 

1« CI((Vl(t v J,KI 9 K7t 9 KNltJ=l 9 JNI,iri v Ihl v 

1« CI C<V<(X t J f N| f KZl t KKI f J-I,JNI f ISI f IM f 

2? CI l»AMI.j)«J±l 9 JN) 9 T:lilNt 9 

21 CIIHX CT t Jl v jrl v ^M*|ri v ffilt 

22 CI IHT IT, JI»J=UJN| tfriiXNI, 
2 3 Ct lHOlT t Jl t Jri»jN|»iri t lNl, 
2« CI IMVII.JI, Jrl,JNi t in,IN» f 
2? Ct l HI ],jl 9 j:i t jNl»I = UIM 

26 WRITE Iflf (If yf I 9 J 9 m 9 N=l 9 ^N> 9 j=l 9 JNI 9 t:i 9 THSt 

27 C(liUPCl 9 J 9 Ki f »<rl f KN) 9 jri 9 JMi 9 t = l 9 jr-i f 

2« C( l<TilI 9 J 9 Kt 9 K:t 9 «*> 9 J = l 9 *iKl 9 iri 9 lni 9 

3C Ct(<T;tI 9 J 9 K),K = l 9 KN) 9 j:i 9 JH) 9 X:i v X>*| 9 

3! CI UTflI.J v KI 9 K = l t KM v J = l 9 Jh I 9 I=1 9 IM1 1 

J? Cll£T*l<t,Jl 9 J=X 9 J-l 9 I=lifHl 9 

1J CHET«Cl 9 JI 9 J=l 9 JH) 9 t-l f fMl 9 

3» Cl CNCXlt Jl tJ:itJMI ( i:ifZ!lIff 

:« CII^CilI 9 JI 9 J=l,JH» f lri ! lNi 9 

3a ciiMeUrJi.jruuM tX7i. t irtt 

il SiHBUS'.il'.illtfl: 

3Q ClTIDF3Ml>.r-l,lW>t 

• i Ct TtO€3S <ii.i=i»1h> ■ 

*! UR1TE HI TT0T 9 NCASC 9 0T 9 STAGC 9 EST 9 A*PLIT 9 PH*SE 9 DPHAS£ 9 

*» CPERtCD 

*»* £NQ FILE 

** WRITE«6 9 112I TT0T,NCASE 9 N8L0K 

*7 112 FGftniT<lX v *0AT4 RECOROEO ON TAP*, TTOTz* t riO ♦O 9 

M C» CASE fc8P=*tTS«* 6LOK NBR*» 9 I5> 

*•* RCTUfH 

SC ENO 



104 ORIGINAL PAGE IS 

'^ PWR QUALITY 



•flCW(l).22l roR CPCATCO ON 4 OCC Tf AT UtJVrQS 
] C««tt**Oet£fiMt'4C$ TMC INCLF OF ftCSULTANT WCLOCITY 

sueftcuriNt ?ri(UtV«?cot 

IF(MS(V»*LT. .QOOII GO TO 10 



iQ 



IOC 



CCNTlhuf 

tFfu.GT.o.) zrori.sn 

iFtU.LC.O.I ?£Oi*.TlJ 

fftTUt* 

C*0 
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•PtOUflUAMTK ST* CREATED ON 2? IP* 79 AT 10:28 i*5 


1 


0,0,0,1,1, 1, 
1,1,1,1.1, 1, 


a, 0,0,1, 
► 1*1,1,1, 


:i:i:8:8 


T 


liliMtltti 


,1,1,1,1 


,1,1*1*1 


• 


1,1,1,1.1, 1, 


► It 1,1,1 


►1,1,1,1 


r 


1,1,1, 1.1, 1, 


► l.uull 


►1*1,1*1 


i 


1.1, 1. 1.1, 1, 


► uuui. 


►1,1,1,1 


■, 


1,1,1,1,1, 1< 


uuuu 


►1,1*1,1 


• 




1,1.1,1*1, 1, 
uuuuu j< 


:i:l:l:i: 


► UUU1 
►1,1*1,1 


1C 


1,1,1,1,1. 1, 


► ltl,l,li 


►1,1*1,1 


I! 


1,1,1,1,1,1, 


► 1.1, Ul, 


► UltUl 


ir 


UUUltUl, 


a, 1,1,1, 


►1,1,1,1 


u 


?:S:l:l:i:J: 


.1:1:1:1: 


:1:1:1:1 


15 


UUUltO,C, 


► O.UUlf 


►1,1*1*1 


H 


UUUQtO.C, 


► 0,0, 1,1, 


►1,1,1,1 


17 
1' 
1« 


O.O,C.O,0.Ci 
UCl.Otl.Ci 
12. ,11. »8.C, 


►0.0,0. c, 

:l:8:i:8: 


^i:?:2:f:8:?:8;i l :8:i:8;i:8;?:?:S:§ 

► 3.C, 5 .0,5.0, 2.0, 2.0,2.3, 2*0*2 .0*4.0 


ac 


8.G,8.C,9 # C, 


► 10., 13., 


► S.C, 8.0,8.0, 6. 3*3.3, 3.0, *».0, 5. 0,7.0 


21 


10.,S.0*d.C< 


► 4.Q,1C, 

► 8.0,9.0, 

► 7.0,7.0, 
►5. 0,3.0, 


>9 .0*9 .0*9. 0,8 .3*6 .0**1.0 ,6.3 ,7.0 ,8.0 


22 

23 


9. 0,9. 0,8. C, 


►8.0.*«3,8.O,«.0*7.g,6.C,7.3,8.0,9.U 

► 6. 0,8 .0,8.0, 9 .0,8.0, 7.2, 7.0, 8 •0,9.0 

► 3 .0,5 .0,9.0,8. 0,8 .0,7 .0,7 .0,7.0,8 .J 


21 


2« 


3.w,*.0*3.C, 


► 5.0*2.0, 


► 2.0, 2.0,7.0, 8.0, 8. 0,7.0, 7.0, 7 .0,8.0 


2* 


2.0*3.0 ,2. C 
2*C,2.0,1.C, 

l.o.uoa.Ci 


»H.0,*.0, 

► 0.0. 4*0, 

► 3*0, *.C, 


:i:§:5:S:l:S:?.1:?:8:7 7 :8:?:8:?:g:?:§ 

►:.C,2.0,2.0,6*0,7.0t7.0t6.0,7.0,6.0 


2° 


C,ti,C.C,3.C, 


► 3.0,3.0, 


► 1.0*1 .Q,2.0,3.C, 6. 0,5.0,5. 0,5*0, 5. u 


3C 


UC,3.3,3.t 


►3.0,3*0 


»O.S,U0,2.3,6.C,3'3,3.3,i».0,(t,C,(».C 
► O.C*C.gt6*Q,2*0,2.0,3.G,H.g,*.JJ,3.0 

0. G.O «C, U.O, 1.0,1.0, 2 .0,1. 0,2.0, 2. U 


» 


l:Z:i±{:h 


'oroide! 


3? 


0,0,0*1,1, li 


►2*0,0,1, 


► uo, c 


3* 


1,1,1,3,3,2 


► UUU2 


►3,u,C 


i: 


i:!:i:!:!:i: 


•J: ?:!:*: 


\UM 


27 


3, 3,3, 3, J,:, 


J:i:J:§: 


3,3*2 


3» 


3,3,3,3,3,3, 


► 3*3*3 


2° 


3,3,3,3,3,:, 


3,3,2,2, 


►3,3,3 


HP 


3,3,3,3,3,3, 


► 3,3,3,3, 


►3,3,3 




J:i:?:5:J:S: 


i:i:J:i: 


-t:f:i 


*: 


2.*. 2. 3, 3, 2, 


3,3,3,3, 


3,3*3 


tit 


G, 0,0,3,2,:, 
CO, 1,2 .5, C, 


3,3,2,3, 


3,3,3 


«« 


2,3,3,3, 


3,3*3 


*»fr 


1,1,2,0,0,C, 


0,2,2,3, 


3,3,2 


*7 


;,2,c,c,c,c, 


► 0,3,2,2, 


►2,2,2 


«l* 


0,0,0,1,3,4, 


0,0,0,1, 


2,CC 


H« 


1,3,3,3,3,:, 


3,3,3,3, 


2,C,C 


5C 


1,3,3,3,3,2, 


3,3,3,3, 


3,3,2 


51 


1*3,3,3,3, !, 


3*3,2*3, 


3,3,2 


5? 


1*2,3,3,3, ?, 


► 3,3,3,2, 


►3,3,2 


53 


1,3, 3, 3, J, :, 


3,3,3,3, 


►3*3*2 


*• 


1,3,3*3,3,!, 


3,3,3,3, 


fZ*l 9 Z 


55 


1,3,3,3,3,:, 

1,3,3,3,3,2, 


3.3,7,3, 


3,3,2 


56 


3,3,3,3, 


3,3,2 

3,3,2 
3,3*2 


57 


1,3,3,3,3,2, 


► 3,3,3,3, 


5? 


1*3,2*3,3, :, 


► 2,3,3*3, 


re 


c,c,c,i,3, :, 


► 3*3,3,3, 


►3,3,2 


61 


C*L,U2*C,C, 


►1,2,3.3, 


>l f 3*2 


61 


i,i,:*o*c,c, 


3,1,?, 3, 

► CO. U2, 


3,3,2 


6: 


Ui.CC.CC 


►3,3.2 


6 T 


CtCtCtl ,1, 1, 


► 3,0,0,14 


►1,0*0 


6* 


1,1,1,3,3,: 


► 1,1,1.7, 


► 3,0, C 


6« 


3,3,2*3*3,: 


► 3,2,3,2 


►3,1,1 


66 


2,3,3*3,2, 2, 


►1*3,3,3 


►3,3*3 


67 


2,3,3,3,3, :. 


► 3,3,3,2, 


►3*3,3 


6* 

6* 


i:S:i:!:i:i, 


► 3,3,2,3, 

3.2.2,3, 


*S*Zt2 
►3,3*3 


7: 


3,3,3,3,3, : 


► 2*3,3*3, 


► 3,3*2 


71 


3,3,3,3,3, : 


►3,3,3,3 


►3,3,3 


7? 


3,2.3*3,3, 2, 


► 3,3,3,3, 


►3,3,3 


7? 
7* 


1:1:1:1:1:1, 


'l:i:\:h 


'1:1:1 


7« 


0,G,U3,2f <, 


►2.3,3.3, 


►2.3,3 


7* 


ui*2,2, o; c, 


2.2,3,3, 


►3,3,3 


77 


3*3,:*0,C,C 


,3,:, 2.2, 


►3,3,3 


7« 


2,i,S,C,C,C 


rC,0*2.2, 


► 2*2.2 
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7« 

i? 

• t 

«? 
9* 
t« 

9P 
ft 
92 
93 



C,0*0,l*3, 


!:!:§ 


1,3*3*2, J, 


:«3*3 


1*3 » 3,3 *3* 


m 


l,!*3-,5*3* 


%3,3 


1,3,3,5,3, 


2t3,3 


1,3*3,3,3, 


.\3,3 


t;];j;i;l; 


J*3* J 


Q,C,3,1,3# 
0,0,1,3,2* 


*:!:* 


1,3*3,2*3* 
1,J*2,0,C. 


t:M 




•Hi 
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•fLO* 


>tlMPC»l St« C«€ATCO ON 19 JUL 7? AT 11:16:2* 


1 


UU 


3 


1*,* 


* 


13, «» 


S 


1 1,3 

1«,3 
l*,l 
16,1 


* 


s 


9 


16,3 


10 


1?,3 


{} 


15.* 


!*.*• 


li 


11:1 


is 


13, 1 


16 


1*,* 


n 


16,6 


19 


IS,* 


i? 


15,9 


i*,9 


a 


16,10 


:: 


3,1H 


§2 


;.ij 

1 ♦ 12 


25 


i,:o 


lb 


:,n 


:t 


2,7 


ZS 


l.T 


2« 


1 .* 


JC 
21 


$•5 
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• rLOWUI.CJQC? }TW C»CATfO ON 23 JUN *0 AT llxX3|£6 
I 2S.O,.e2,J*7.6,?!3.C t .?$,26.6 



6 
7 
P 
o 
10 
11 

ir 

i* 
ir 
is 



26.1 
27.2 
28*1 
21.1 

29. «• 

30.0 

38:! 

30.0 
29.9 

la. 3 

ft* 




90.0,.*U,26.7 
110., .30,26.1 

110.. ,;u,26,« 

110. ,.93.27.0 
110. ,.90*27.0 
90.0, .60. 27.Q 
90.0, .50, 27. 3 

M8:::i!:i!:? 

ICO.,.30.27,2 

110., .20.27,1 

:3:8::&!:3?:8 
iS:S::K:i!:8 
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*rio»u i. 



H'J 
7. 

• • 

t 7. 

<7. 

I ' • 
i7. 
£7. 
i7. 
i7. 



U I * 

lit* 



7 • y , 

?«:♦ 

7)wi 
7 tut 
7.L, 
7*Ct 
7.0, 

!•*• 

7tUt 
7*Of 
7*C t 
7,J« 



17. J. 

:7,c, 

4 7 f wf 

47. C, 

;7,c« 

:7«u» 

;7.j, 

27, C. 

t 7. Uf 



TtO 0»l I 7 way 
* 7.C»* 7 



r . u f 
f * w , 

r.n. 



;7«3t«7 

*7.Ct^ 7 

27.3,27 

27.^,27 
*7«J«i>7 



id 

• Ct 

•a 



7. 
7. 

* I • 
?7. 
27. 
27. 

!!. 

n. 

27. 

27 

27 



7* k 
y t < 7. 

ti,27. 

U,27, 
J, 27. 
L v 21« 

u«*7« 



T ]i 

w ,4 7 

Ct{7 

S:? 
§:!] 
S:3S 
8:1} 



:i«tS9 

• b |. 7 i 
«u »*7 . 

• i ♦ *. I • 

IW|4I * 

.C*47. 

«J»27. 

• #*7. 

• U fit 7 4 

.0,27, 

• t*7 • 

• %i 7 . 

• ii 7 • 

• vi«4 7. 



w»2 7 *w < 

i t«7«C< 
L,27.i, , 

Wf 47.0 ( 
Lt<7.Ci 
U2?«(m 
Li*7 «C 

w»27.C 
w,27.C 
0,27.0 



: t n.c,2 

» #4 7 • 3 ,2 

[ »4 7«Ot< 

It?*. 3. 2 
t^8.2,2 

I::!:*:" 



:i?: 5 




7,C,27 
7.0,27 
7.*, 27 
7. £,27 
7.C.27 
T.Cti7 
*«C,29 

7.7,|i 
7.t,2l 

7«« f 2« 

7.0t|7 
7.C27 




.5,27*3 



no 



4«f LOWU t. 




i?.w*t?7 
t7*L«27 




«7«««4 




SO -T : 3:! Cl*C 

?:!i T :?:5? 






i 

iiii 

9*41 







1.27,3 



I f*4«7 f 



;7.o,j 
iT.o.a 

:S:J 
H:|: 
B:: : 

24»lt 

*7*Z# 




111 



Z J « C $Zl «Ut 4 7 *Cf 



I.e.??.? 



ICO 0« 1 
i / «0 1 2 7 * 



*»i 



<'iWtiTibf«ritf 

<7«u f X7tw»27«C| 
*/ • Lf ?7 * Jt< 7«U 9 




13 




,7.Ct27. 
i7«*.47 



9 *AY 

ww 

C«?7«y 



7* AT jj*5 
«27*C»:7«0 




W7.C.17.2 



ftt'7 

• 4 /• ifZ 7 

•*Wi? 

tc*«u»29 

: :t:il 

•*7«fct*' 



C, 27.0,27 



.c.It. 



0,2J.U,iT.3 

\it±h:l 

,27. a, $7.0 
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«ru«(l Ur^lTt 



to 




;:«-,*7.u»<7o. 
;t«u*<i*Ci«*7«w* 




:I:!:JJ:J 

7»C t*T« 
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•nottUt«'LOT»'* ro« c»catco on n occ if • ? utscto* 

? C MfN MC^ttA* '0* AWAt?*XN« ANO PtOTTlN* fMC »I5UIT1 OSTAXNCO It ANCH* 

« P**kPt1Z* tN:U,J«Sl% # KMSS«THSl!( 9 JllS|3*NCMSSt*HCTSlS 

5 OIKCMICN tfiyni0C0t 9 fLMI«TN(4l 9 riMI«0lMCYI 9 e2fltCY4t0JCIlCYtt 

? CO««*CYI,Q$f«tYl f fAll«| _ ^ . _, — * 



OXtttMfe* mOTfS) 9 *»lOUT<*CYt 9 «STA*lOtffCYI 
I! OI"CM?0* TtllH 9 J* 9 Mlff2«l»tJM»«NltU«X*t«IWtKNI 

1? Bl^SiflS HMISlISt;l|gfi3}lC{-?I8(»HS:;USiKiiCsSS 

\h OlnthilOH XtCDICNItfYiNCNT.CTALCttCTItCTAHf^CYI.CUXNfCffCYl 

iS 1KB i:ilfc!ll:i51:Hil 

21 HCAO ZtCTKTt.lrl.NlRl 

22 ftCAC 2 t lTAl(Il t ls| t Nl»| 

8* «<*o r v cc;m,t:i»Mt»i 

2* dCAO 2tf *?3II It Xrf 9 NIAt 

27 «2*0 2»<OS<Xf«I?X«NtiU 

$• «?AO 2«<XPLOim f t:i,S) 

|t ft?AQ I.tNJtANOlfiarijKCtl 

|r «ca5 3,ccT*uxi«|:i*»itfi 

ij hcao 8*cctahui.x:i»m!iii 

;• ft'AO 2,lfT4lNT(Xl.i:i f MftJ 

*\l$ h 0». 0Y 
OC ♦ NCsXtNCft 
tCAO 2» IfJ 

xrf*ci:riCAni-! I 9 

TY|NC1ircOATU*ij 



X7 OC t NCsXtNCft 

*• " • 

2f t tOMTTNUC 

*: 2 fount ii 

*? *t 7s« i^-i 

*• MCT:l, 

K* CICt rACIO»«FACTJ 

*4 iNK:fCOAT(ltN*n 

«7 OCUMCC* 

4t APMMr.C* 

«« AffHA >S.1S 

«? o«=.;s 

5! AiatftAOCCCO. 

«? PUsC. 

«r coir :ros uhii 

«< lSCAl£iKtyOV/rLOAT(TN*ll 

56 Y3CAt?ilSC*tC 

57 2SCAlt:,2G* 
S> OjCAlfi.OCT^ 
«* VSCAl'-USCALt 
4*" bSCAl*::. 

ti i$c=K/r.s*/*scAtr 

HZ YJCr t./J.{«»/Y}C»tr 

*! 2JC5 j«s2.i*/*fc*i* 

4« 5C- I./2.5«/U5':ac* 

4) V5C5 !./:.t*/V!C*LC 

ie Mse;i*/2«s*/i»S'*i > . 

j; ice rS«fUfi;^*?l*:?l?:!J« t«:«,ru.:,' rso'.nz.zi 

4 oaro« , .«r?{5i. , e $ .5:?fe:Si5 iK«..rw^.' niosru.fi 

7j call •toTsiiiur 9 I:c%tii 

77 CALL »L9Yn.,l # ,*1l 

7! C«»»*^LO*ritCTMfI*S*?0«*7HC*X<P*oS^ t%A* IriT' »*»OLATCO T«riCtO 



74 OC a *:! .M* 

77 OC j JM.J* 

7« *£A$ 2*CU9<X*«J*M 9 X:i 9 Z»l 
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?• i comhui 

ic xnx'toxt.co.ot co to ? 

«* 00 U M:l«*Xft 



12 ;i:Ojt*l 

• S MlOMHI 

fC Titortrimi 



ill 



:ali •lot<io««q« 9 «i) 



«? 

«« CAtt CA*TN*«* # HI0C*l 

i;; 12 call rAUHJtuo 

l8l _ CILC »lOT<0. 9 0«,Jt 



CltL CA*TNf *-- 

m •• IV& SttWi 

I8f ^...JKi!«[j;« 

!;i c «•* *C*QLi N "WW"! 1 **° "*" c * c cw0,c * or tfM * CH Tc •' * t0TTC0 

1C« CAt. H*OT«lStJII«Miltt9«WM.Ut»(l2#SI«V2fVtCTfltTl.Utr|. 

JC« CtTAiM^T.^MtH.MS.IIt.^OjHUfHV.-Cl.^CttMCASt.WlCOIltTTOTtOT* 

J 'f;R«K»Wrt!^ m to 1000 • 

11! Nrft.$lA*o<*CTH 

il« 00 jc5* t*i,X" 

Jl* 00 ?C?k JSl.J* . 



ii* ir<»*xit*^i;co*Q> T9(X«j«ii:tai(ni 

iij 30 c* continue 

M=OMnGI 

M£?*f«6l 

* 4 ; MsC'fKQI 

12* JU**C. 

K! «:M10T 

Uf 00 t<^ Jsl*J» 

12* 00 t«* I-ltl* 

!{• Iuw::uHMMu t J,M-T;«u f j f *it«*2 

i!i 



!?{a^i?t^T^!«2:SV^C0 TO 1900 



tit . ... . .. 

n* CAtt •toni9«tS.«*si 

15« M?Ai:?« 

fst.X* 

i?« XM*»t«ti,jt.ro.ci *e to 



IN 83 1 



!!• AMrMM| 9 JI*CTAIt { J| 

(*% H"Ai:A*AlltAH*H«'Af I 

1*1 C*«««*»tOT< HfSUtTA*! VCt OM H I 0* AfcAOtT HORIZONTAL $XC** *tA«C 

1*4 Cll,T«.AMt« # irrt7!2^^#i^iN i Atf*Al.C0lr 4 &lM/ # ftO 9 M» t UOU«l* 

1»* CJCL**»r4tU»»*k#0> t Ot f M«*»#MKT,«rtOT« 

isc CAtt ei»rfi4cM t i«tBeuf 

l«t CAtt CI»T»*?INi 

is* c«tt *tcMc..o«,3» 

i«? call •tof ui.#c..»!i 

{«* C*»«OtQM «(«L'tr2Sf Vft cr u u * o* t-* VC'UCM, scctxans 

15« 10C1 XMXkPHJf.lO.OI fti» 10 130? 

1$a CAtt ttOrutffXtt.JN. «*'.lf. J*»«*.-.*2tC*tO»«OS« --( ^ 4- 

15 f C^L0r»T«MltAf f y^AX 9 W«*AA«KCt v HU«A« tf XMtAllKAy 9 XL*C0Sr»SXMr»H«fftA«« 
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\% ^aL Ti ^VTHj;i^j* i ,;«»«3s:iii:s^'"« c *««"*«««»«"'i» 

16C CALL CAPTNSr« f HlOCMl 

161 CALL CAPTM6IM1 

call n 



16? CALL 'L0T<a, v O.,3> 

16* CALL PLOT(lO.,Q. t -3l 

16* C«**9«PLQTS PSSULTANT V*L OF V C 



16? 1002 IFUFLOTUI.CO.C) 60 TO 130? 
" "-v«(i 

AX 9 U 
LC.r 
NltP 
HtiH 
H7tK 
17? CALL PL0TCC.,0. V J1 



W ON N-S VERTICAL SECTIONS 



65 CALL PLOTVV(IN,JK*KK»lH»J^ t VetU f K2 f Cir t OV t 0S t 

16* CNPL1 



#LGUT .MWAX.UMix^HAX^iXtMV.ifiliujAPM^x'xLtCOSriSIiir.HB^Aa* 



16" CALL CAPTNlCPl f P2. . 

17? CALL CAPTNSlNtHlBeul 

171 CALL CAPTN7CM! 



\\l C lSoHwhtS?5l!io5K^ $#f T0U * f "MOM VALUES HOST SC GIVE* INOXV 



17* 2LIT:tTALCN6| 

177 2FIG=CTAH<NGI 

17* SAMCCNrZLlT 

170 CINf reTAINTCNR) 

18? CALL eCHKCNCETA t I* t JH t l f I* t LJH,<i.a t $.6t«0«#SAMC0NteiNT f 

iai SBS*!'* 1 ?^!^!*?!!?*!^?:! *^^ 

182 C 

1*6 C 

18« CALL CAPThaCNI 

186 CALL PL0Tn. f 0.,3l 



CN,tXtOY f XSCALrfYSCAt£ f doSr f SINr f XX»rtl 

ALL CAPThaCNI 
ALL PL0Tn. f 0.,3l 

1«" KIKPIOT 

191 OC 9*8 I=l f lP 

19J go 9%6 J=l.JM 

19! 99« £TAC T,J|rfflCl f J t Kl 

i9» mrmfNci 

19« 2?IG:TWINGI 

19* $A«CC*=2LJT 

1?7 TlNTlfKNGl 

ill gitl: g^«l« , i:i;!g£Si- p, " P5 ' *- J,sc ^ scl 

2C? CALL eAPfMiSUK) 

2C» CALL CAPTK9CNI 

205 CALL PlOT<C.»0.,3l 



2C6 CALL PL0Tng*.C.t-3) 

2C7 991 rcfcHlTUX, 15M.2I 

2 0? 992 TOP- AT| r8. 2,' OC VI A TIC 



,C9 iOQO CONTINUE tON H...T. I. T£*P AT- .IS.- »<OCG-C. %Fl ?.Sl 

2iC 1« CALL PL0T(0. 9 0«*9991 

!W 10i 11V 
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**XO«U).PLOW FOP CREATCO ON 7 0£C 79 AT 1C:05;19 

C •••••♦•••«•••••*♦•••••*••••*•••*••••**••*• ••*•****** 
C TO PLOT RESULTANT OF U C V ON SIGMA PLAME K2KPL0T 



I 



7* 



• SUBRCI'TINC PtOTUVCIM, J*.KN t I«, J*,U6 •V**fHEX,P£ Y t NCN f XX,YY , AKN , 

« Cr€L7|7PPf»*AffHlN t AaMAX t C0SPtSlNFtKQtH8«UOUHl,U0UH 9 XttUAAX t 0X*0Y t 

6 CH*«X ,'4PLT,KPLOT t 

7 CXSCAlf.VSCILExZSCtlFtU&CAU.VSCAlC.WSCAU.STafieflTAI 

CXXC»t*>,YYINCMl,HPhH»4Wl.£TlTl*.*llU 

IP MPITfl6 t 32» 

11 32 FC«M#T(U,«UV PLOTS*! 

17 A l=CT*X$CALE/2« 

13 MJUPIOT 

1* IT(H .GT.ll 60 TO 20 

IS DEPT»>:Q. 

lfi 00 3C isltl« 

17 00 3C J^ltJH 

1» I?t*£*<I f JI.EO.O* 60 TO 30 

10 AIrFlOATi!*ll*OX«XSCALC 

2T AJ=F10AT(J-1I*0Y*Y$CALE 

21 AA?.:*I*UB(T»Jt»l*U$CAL£ 

2? AAj:*J*V9ti.J,Kl*V$CAl£ 

27 CALL VECTtAl t AAl f AJ f AAJ f ARKIN f A«HAX*COSf tSl*F» 

2* 30 CONT INUC 

25 60 TC 8 

?6 C«*«« «•••*•»**•*•*************••*•• 

27 20 CCNT 1»:UC 

2« OEPTl-r0£L2*FL0ATlM«ll 

Z9 00 «C lrt t I- 

2C go. *c j=i,j* 

21 IFt«CX(I,j).CC.Cl 60 TO 40 

2? AN=MErX«J»»CTA(X.J} 

3? IFlO£PTH.G£.AHl GO TO «C 

3* 0077 JH/AKN 

J? Ll^l «INT(Q£PTH/0021 

36 298 L2IL 1*1 

27 L2rL««l 

3* IF(L2.lC#KMl GO TO 299 

29 L1=LW 

«C 6C TC 298 

*1 J99 CONTINUE 

*? Z1=OC?«FLO*T(L1*11 * 

*«3 22-21*007 

*»*« 23tz;*oo2 

m< euuciltJtLii 

it E2=ue<I,J,L2l 

M CALL Flt{h,72.23,Cl 1 E2 1 E3 1 4 f B,Cl 

no UO£PlMrA*0£PTI- *--*..-* 

50 Ci:VHX»J t Ll> 

51 £2?VH 1,J V L2> 

I! 

S* AI=FIOAT(I-1)«OX»XSCALE 



£3£UFIX £ J.L?I 

CALL Flt(h,72.23,Cl 1 E2iE3 1 4 

U0EPlMrA*0ePTH*»2«e*0£PtM«4 

EUVMl t J f Ll> 

£2?VH 1,J V L2> 

^lLL e rIfy5V!:2,73 f El,E2.£3.A f B,C) 

VCCPlMrA»0EPTH»»2*S*0EPTH*C 

AI=FIOAT(I-1)«OX»XSCALE 
Si AJ=FcCAT(j-l)#OV*YSrALE 

57 AAI=II*unEPTH«USCAL^ 

5? A»jr*j*VO£PTH»VSCALE 

*9 CALL VLCT<AI t AAX t A4*AAJ f APMlNtARHAX f CQ$F,$I«FJ 

6? *c c:ntinuc 

6 1 8 CONTINUE 

62 C**»»*0PAW5 OUTLINE OR «Ol'NOAPY OF THf INTEREST AP£A 

6' C*LL PLOTt-ll,-Al,-T>» . 

f* CALL OUTLIMtIN, J*.tKN f NCN»OX t OY.X$CAL£tYSCALC t COSF,SINF,XX»YYl 

*5 CALL »10TI Al f Al,*31 

6A NPLT:?lPLT«t 

t? W^ITE(b*la)^PLTfOCPTH 

6« J* FCH-IT( it, -otCT ?*«« = •, I6»* COMPLETED* OtPTM5 t , FlO.Qt • C"M 

*« 13 CONTINUE 
SOO CCfcTlluE 



1\ Z FOPPITf | 

7.: RETURN 

n £no- 
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II.PCOTUH FOR CREATED OK 19 OCC 79 *T J?:UM* 



C PLOTS RESULTANT OF U I % ON SELECTED E-« CROSS SECTIONS 

SUBROUTINE Pt07UwnW»J*i.HS,tM,J^ i Ue,H,«?»Ol,OY i OS. 
CPtOT^T f H^AX v UMAX t W««k v HeX«HU 9 i(l v XN f iftn*X«XLtCOSFffilHr v HRffHAII # 
CNPLT tXSCAt€»TSCH€t?SCAL£»USCALC 1 VSCAC£»«SCAU^»STAG€teTi»l»Xl 

OI*E*StON "AftUN, JN> f H6<lH,J*),ETAU»«*JH> t»rtlN,JN> 

OXHC^SIOM U8(lM f J- f KN) tWtlMtJM.KHJ.HtxaH.J^JfHUflNtJNl 

WRITE 1**321 
32 f OR^ATltX^UW PtOTSM 

K?tK*<»l 

AUCMXSCALE/?. 

do i erg xj:i.3 

IFtI*.E0.3l JSJ2 
CALL PL0Tt0.,M.5 t *3l 

oo us r=iti« 

lFfHfXII,J).£C.01 CO TC 70S 

Al=f tOATM-1 J«OX*X$CAL£ 

AHSMttZ, Jl^TMl, J> 

OC TC5 *n,K2 

AK=<<TAGE«£TA<I»JI-FL0AT<tt«ll*0S*AH**4SCAL£ 

AAl=II«u8(I,U f KI*USCAL£ 

1AK:|K«UIT >J 9 K I*WSCAL£ 

YU?.2«SQPTMAAI*AII**2«(AAK*AIC1**2I 

YWrA^AXltAP«IN t AMINnTW f AfWAXIl 

CALL AP0HOlAX 9 AK t AAI v AAK # YU«Cr«C«l2> 
7C6 CONTINUE 
705 CCNT ]NUt 
C*****ORAW< BOUNDARY OR BOTTOM OF TH£ CROSS SECTION 

CAtL PL0T«-Al,C #> -3! 

00 710 I=lfIN 

XF(KMX»J)«EO«C) 60 TO 710 

NN^NMl 

JFtNt.GT.I I 60 TO 711 - 

Ai:ri4ATCX*ll«0X«VSCALE 

AK=irAGE*ETAtI,J>>*2SCAL>e 

CALL ?L0TtAl,AK,3» 

AKS-ZSCALE^CHUtt.JI-STAGCl 

CALL PL0TtAl # AK # 2) 

GC TC 712 

711 C0NT1VUE 
AXrFLOATfl*l)*OX»tSCALE 

AK = -<SCAL£*<MUU.Jl-STAG£l 
CALL PLOTtAl*AK ,21 

l\l\n 

712 CONTINUE 

713 CONTINUE 

1F< 1.LT.IN? fiO TO 707 
AK? (STAGE +0.Q)*2SCALC 
GO TO 708 

AKriSTAGE*ETA(I»JM*2SCAC£ 
CALl PLOTtAlO,AK*21 
WPITEU.T7I J 

77 rc&vtrtlt* •jr',191 

FJ=F10ATIJI 

CALL SYMBOL <S.« » -n • 7 f .l,;Mjr f e.t2> 

CALL VU-IEJm.T,-!}.?,*!,?^,:?.*-!) 

CALL PL3T(AUC.»-5) 

CCMJNUC 

CALL PLCTC0, f -*.5»-3> * 



7C7 
708 



13C 



MPLT:*jfH.T*i 
-■HIT"- "■ 



walTC(*,UI NPLT 
13 f CR?»ATUX t ^LOT NBRSIl** 
RETURN 
£H0 



COMPLETED*) 
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U 



L0UI11.PC0TVW FO* CREATEO ON 19 OEC 79 AT 12:1C:*1 



PLOTS RCSULUNT OF V C^U ON SELECTED N-S CROSS SECTIONS 

SUBPCUTINt PLOTVUrtNtJN.KNtlHtJ^.VBtUtKr^OX.OTtOS* 

CNPLT.xSCALE.TSCILCt-'SCALE.USCALj.VSCAUE.USCiL^sfACetiTAtNYl 
DI*E*SlON M*R<lN,JNI 9 H8?Ifl 9 J''I 9 *fAUH 9 J»> 9 KV(lN 9 UN) 
OIl^SION Ve(IHtUM t KN) 9 tfaH f Jtt t irNl v PexCZN«UHttHVClN«JNl 
U51TeC6.32> 
32 FQRt-ITux 9 «VU PLOTS*! 



a!io>* 



XSCALC/2* 



DC 1CCO P = l f 3 
IFir " 



«02 



_ iH.ecii u* 

IF(H.£0.2> t:« 

CALL PLOTO. v *l»$t*3) 

00 8C? J=l tJ* 

IFlHCXII.Jl.rO. 0| GO TO 

AJ=rLOATCj*ll«OY*YSCALC 

AMrHP(l,J)*CT*(I f J) 

00 8C3 K=l 9 «2 

AK:t<TAnE*ETAtI,J)<-FL0ATtK>lJ«O$*AHI*2SCALC 

AA % S*J*VB<I.J»*UvSCALE 

AAH=*K*W(j t J t H )»WSC*Lt 

YW;.:*S0RT<(AAJ-AJ>**2«fAAK-AK)**2I 

YW:AMXt(A?t<IN t AHlNmw 9 ARHAX> I 

802 CONTINUE 
:«**«*ORAU< BOUNDARY OR BOTTOM OF THE CROSS SECTION 
CALL PL0TC-Al t 0« f -3! 

00*8 10 J51.JN 
IMMtCI»Jl.EQ.O> 60 TO 810 

IfTna.GUII 60 TC 811 
AjrrtOATlJ-ll*CY«YSCALE 
AK=| *TAGE*£TA<I 9 J) )*2$CALE 
CALL PL0TfAJ t AK t 3) 
AK=-7SCALE*tHVt!,J|-STAGe) 
CALL PLOTCAJtAKtil ■ 
GO TC 812 

811 CONTINUE 
Aj=FLOAf(J-t J*OY«YSCALE 
AHr-asCALE^ihVCl.Jl-STAGCJ 
CALL °LCT(AU v AKt2) 

JO=J 
AUO-'J 

812 CONTINUE 
810 CONTINUE 

IF<J .LT.JN) 60 TO 8T7 
AK:C !TA0E«3*Q1*?SCALE 
60 TC 806 



837 

aaa 



AK=(!TA6E«ETA(1.J1)*ZSCALE 
CALL PLOtiAJO, Alt r 2l 



URXTE<6.77) I 
77 FCRHITI 1T 9 *ir* ,T« I 

rimoATm 

CALL SYHBCL(S.5 f -g.7 f# l t ?HI: t 0..2l 
CALL NU«9ER«S.7,-0.7 f .l t f I t O.»*tl 
C*LL RLaTCAl t 0-t-31 
1000 CCNT INUE * * 

C*LL PLCT(C. f -*.5 f -JI 
NPLT:NPLT*1 
bRXTCCtt!8> NPLT 
18 FORHlTilX t »PLOT N?R* V 1* 9 
RtTUFN 
IHO 



COKPLETEOM 
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*nowm.ouU!N rot cae*Tro on to oec 79 *t 11:22:23 

} C**o5»Ss*OtT*l5t/n3NoI»Y*OF*THC*rNTc5nT IkEA 



SUEftCUTINC 0UTLlN(IK t Jh 9 KN 9 NCN v 0X t 0Y v XSCAtC«TSCAt£tC0Sr t 

csrwr .xx.yyi 

OIMCfSlOM XXCNCNltYYCNCNl 

00 1! NCzitNC* 

A=XX l*Cl*OX«XSC*L£ 

0:YY<NC>*OY*YSCAL* 

MVS2 

IMNC.EQ.ll KV = 3 

x5A»cosF*e*si^r 
Yx-A*«isr«fl*cosr 

15 CALL f*L0UX t Y,ft<M 
HCTUCN 
CM) 
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ECHKON ELT CREATCO ON R 0£C 79 AT 09x59:18 

SUBROUTINE £CH*ON 

1HIS IS ENTRY SUBROUTINE F0» NHC CONTOURING PROGRAM 
CCALCCHP 0» MILGO TYPE PLOTTER? 



ELBERT HILL. NHCt MIAMI, FLA. 



FALL* 1970 



Tt>£ COMPLETE PACKAGf CONSISTS Or 3 SUBROUTINES, ECHRON, CONLIN, ANO ENDE* 
ACL I ARE CATALOGUED TOGETHER IN THE UM 360/65 UNOER MOOULE NAME ECHKON* 
AM) OECKS ARE NOT NLEOEO. 

ANY FECTAN6UUR CRIDCEO SCAtAR FIELD CAN BE CONTOURED ON HUGO 
OR OLCOMP TYRE PLOTTER BY SETTING UP PR0PE CALLING APGUHCNTS ANO 
PROCEDURES AS INOICATEO BELOW ANO lH£N CALLING ECHKON* 



-CALLING STATEMENT IS AS FOLLOWS- 



CALL ECHK0NCHM f INl r TN2,NCXUNEX2,NETl,NCY2 t MI t ai0fPLTTNC»SAHC0N f 
2C0NlK,RGRIOtlN3,IN«,2LIT,ZBIG,ANORTH,ASOUTH,AEAST,AW£ST ff NOASHO« 
SNCASl-U,XLAeELtSP00TH;iREceT1 



—DESCRIPTION OF CALLING ARGUMENTS — 



HH I! ARRAY COKTAINING GRID OATA TO BE CONTOURED. ITS DIMENSIONS 
---■■--■-_ JI^CNSION HHCIN1.IN23. POINT 1,1 IS LOWER L£f 
INI IS OTHENSIOK IN X OIPECTION ANO IN2 IS 



ACE INI ANO IN2, 
CCRNER OF GRID. 



DIMENSION TN V DIRECTION* 
tX INCREASES FROM w£ST- TO EAST ANO T INCREASES FROM SOUTH TO NORTH} 

NEX1, N£X2* NEYU ANO NEY2 0ET£*MIN£ THE PC*TION OF HH GRIC TO 

BE USED* KEX1 ANO NEX2 A*E THE FIRSTCLEFTH0ST3 ANO LASTCRIGHTH0ST3 
COLUMNS TO BE USED* NEY1 ANO N£Y2 ARE THE FIR&TC90T T0H3 ANO LASTCT0P3 
RCWS TO BE USED. CTHUS ANY SECTION OF MH CAN BE USE03 

FOR FULL GRIP 

N£X1 > 1 
K£X2 > INI 
N£Y1 > 1 
N£Y2 5 IN2 

HI I! HEIGHT IN INCHES OF CONTOUR MAP BETWEEN LIMITS NEY1 ANO N£Y2 
WIO IS WIDTH IN INCHES OF CONTOUR MAP B£TW£EN LIMITS NEX1 ANO N£X2 

PLTUC IS STRAIGHT L*N£ PLOT INCREMENT IN INCHES TO BE USED 

ALONCi CONTOUR. 6000 VALUE IS .OH, BUT CA*' <?£ VARIED UP 0» DOWN. 
SINCE LARGER VALUES CAUSE PROGRAM TO RUN A LITTLE FASTER, IDEAL VALU£ 
I! LARGEST THAT WILL STILL GIVE SMOOTH LOOKING CURVES. m mm mm 

CO SOME EXPERIMENTING MITH TT. START WITH ,03 OR .06 ANO INCREASE. 

SAMCCN IS ANY SAMPLE CONTOUR VALUE, IT IS USED AS A START1N6 POINT 
FCR COUNTING UP ANO OOWK TO GET OTHER CONTOUR VALUES. 

CCNlftT IS CONTOUR INTERVAL TO BF USED. 

RGRIC IS AN !NTC6fR«2 STORAGE ARRAY USEO INTERNALLY IN PROGRAM 

If? AS ARGUMENT IN . _ 
DECLARE AS INTEG£R*2 



A HO NEED NOT BE INITIALI7ED. IT IS INCLUOCC AS ARGUMENT IN OROER - 

TC TAK£ ADVANTAGE OF VARIABLE DIMENSIONS. 

BEFORE CALLIKG. 



XKJ JNO IN* 



ARE X 
BE AT 



AKO Y 
LFAST 



OIMEKSIONS 
AS LA^G" 



IM *US> BE AT LFAST AS LA9GE AS NEX2-NCXK1 
IM ^UST H AT LFAjJ AS LA9GE AS N£Y j-*E Y 1 < 1 

CTHUS RGRIO MUST BE AS LARGE A$ PORTION OF 



OF Rf.BIC OIMENSION RGRI0CIN3,lN<i3 



OATA ARRAY HH BEING USE01 

2L1T 1N0 2eiG ARE LflW*CR ANO UPP r R C0N1CUR CHECK LIMITS. NO CONTOUR 
WILL BE DRAWN RELOW VALUf OF 2LIT OR ABOVE VALUE OF 2BIG. 
CtlStFUL TO PREVENT QRAWtNG FOR ANY COMPLETELY WlLO 0ATA3 



ANC3TM, ASOUTH, AEAST, ANO *WEST CAN BE USED TO ELIMINATE ANY 
NL^WER OF INCHES FRO* ANY S10E OF FINAL DRAWING. 

FC» FULL DRAWING WITH HEIGHT > Hi ANO WIDTH > WlO, 
TNITIALI7L ALL * OF ABOVE ARGUMENTS TO 2ER0. 



FCR EACH OF THE ABOVE WITH POSITIVE VALUEt THIS "ANY INCHES 
VTLL 3^ EIIMIHA'EO ON SIOC TO WHICH I? '(PLIES. .. . . 

iHtS ALLOWS US TO FIT ANY RECTANGULAR G*IO TO ANY *£RCATOR 
OR OTHER MAP LIMITS WITHOUT ACTUALLT ADJUSTING THE GRID. 
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7« C NCAS^P ANO NOASNU CONTROL TYPE Of CONTOURS TSOLIO OR OASHEC LINES3 

JO C IF EITHER OR BOjH ARE ZERO OR LESS« CONTOURS ARC SOCIO LINES* 

$2 C IF lltM APE POSITIVE, CONTOURS WILL BE OASNEO AS FOLLOWS- — 

!• C PIN DOWN SECTION LEN6TM > NO«SHO«PL TINC CPLTINC IS INCRCHCNT LENGTH? 

F5 C * PEN UP SCCTION LENGTH > NOASHU*PLTlNC 

66 C CTHUS LENGTH OF OASHCS ANO SKIPS IS FULLY VARIABLE? 



S« C XlAdEL CONTROLS LABELING OP CONTOURS* LINES ARC LABELEO 

S« C ONLY IF ILA*CL GREATER THAN 2ER0. VALUE Of XLA6CL 

«T C IS HEIGHT IN INCHES OF LABEL NUMBERS. LINES ARE LABELED 
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OTN TS A CONTROL FOR VARYING CONTOUR SMOOTHING. 
INITIALIZE SMOOTH TO SOME VALUE BETWEEN 0.2$ AND 7.5, 
NV VALUE OUTSIDE THIS RANGE IS SET INTERNALLY TO 1.03 



«1 C WITH NEAREST WHCIE NU*PER VALUL OF C05ICUR . IF, SPEC 1*1 

92 C LABELING TO INCLUDE ONLY PART Zf NUN*C D OR TO INCLUDE 

43 C DECIMALS IS OESIREO. SUBROUTINE ENOER MUST BE CHANGED* 

<• C 

5« C SMOOTH 

*fr C I"" - 

07 C tAN 

9R C LARGER VALUES GIVF SMOOTHER CHART WITH LESS DETAIL, WfclLE 

• • C SHALLCR VALUES GIVE LfSS SMOOTHING AND PORE 0CTA1L. 

1C" C CMRMAL VALUE FOR MOST RUNS SHUD BE ABOUT j.SJ m mm 

Kl C AfYfHlKG LESS Tm»N ABOUT 0.«« CR LARGER THAN ABOUT 3* IS 

1ST C FROBABLY NO GOOD* 6CGIN WITH US AKO EXPERIMENT UP OR OOWN 

1C 1 C TC OCTER M INE MOST 0ES1RABLE VALUE FOR YOUR NEEDS. 

1C* C CI*PUT GRIO DATA VALUES ARE NOT ALTEREO IN THIS SMOOTHING! 

1C5 C 

1C* C IPCCCT IS PLOT TAPE RECORD COUNTER. INITIALIZE TO NUMBER 

107 C OF PLOT RECORDS WRITTEN BEFORE FIRST CALL TO CONTOUR SUBROUTINE* 

ALL OF THC ABOVE ARGUMENTS EXCEPT ARRAY RGRIO MUST BE 0EF1NE0. 
ARGUMENTS ARE NOT ALTERED WITHIN PROGRAM, ANO RETURN INTACT* 

BLOTTER BUFFER SPACE MUST BE SET UP AKD CALL TO PLOTS 
MADE BEFORE FIRST CALL TO THIS SUBROUTINE. 

PLOT TAPE MUST BE CLOSED OUT AFTER FINAL CALL. 

<*Y NUMPCR Of SUCCESSIVE CALLS CAN BE MAOE TO CONTOUR 

•UBROUmC rCNKOu. EACH map becomes a separate plot rccoro. 

*0 INTE#?4AL MAP SPACING IS °ROVIO£0» WITH PEN RCTURNIN6 TO 
CRIGIML ORIGINfLOWER LEFT CORNER} AT COMPLETION OF MAP. 
CTHU5 IT IS STKPLf TO PUT MORE THAN ONE SET OP CONTOURS ON SAME NAP} 

A* v J P ECTAL ^ARKI'NCS OR LAPELS THAT ARF CESIRCO »UST BE OCNE 
fUTSlOE THIS SUBROUTINE. THIS SUBROUTfNf DRAWS CONTOURS ONLY 

WITH INCOMING ORIGIN BEING LOWER LEFT COPNER OF CONTOUR CHART. 

SUBROUTINE ECHKON (HH, I M f IN?, NEV 1 ,NEI2 .NEY1 ,N^Y2 , HI ,WIO,PL TINC » 
?SAMCCN t C0NlNT,Sf;RinaN3,lN«i ^LITjCelGtANORTHtASOUTHfAEAST.AWEST* 
?NCASKn t MOASHU,XLA»EL f SHOOTM 1 IRECCY,lN t JN t KN f 
CNCK,CX,OY,XSCALE,Y$CALE»CQSF,SiNF,XX,TYl 

SEE ABOVE COMMENTS FOR DESCRIPTION ANO USE OF ABOVE ARGUMENTS 

•CN /5TR< 
.CGJG.U.* 
jwiJ ♦ LCLV.P 
*N*Xt I.^YIl, 
50CL»6S,0llt! 
CCMMCN /0£NnEO/HlXEN f WC0E f HOGM,XXLAST,YYLAST 
LCGI(»L PASHEP,nOLA»i,OUTS 
OlMEtSJCN MH(IM f lN2l t fttftI0CIN3«lN*) 9 HCNt*},XXCNCNl*YY(NCN1 

intec;r jGRro 

fi«TA 1ma*»/Q/ 

KHC S13R0UIINE FOR CPMGURING 5C«LAP TIELO ON CALCOMP OR MILGO TYPE PLOTU 
CLPERT HILL —FALL, 1'Oa — 

ItI-#PM 

bfiITM6 9 !C)I 
10 FCRf /H///.2X, 11HCCNT0LR -AP , Id , U ,2*M INCOMING ARGUMtNTS FOLLOW.* 
WPITrU,l2HM,IN2,^EXl,KCX?,MEvi t uET^#Hl,WlO f PLTIhC,SAMCON t CONlMT 
2.IN3 , l f *<» 
12 FCRPITI/ t «« v 63HINl TN? ^E»l NEX? NEYl NFT/ «I WIO PLTINC :*MCON CO 
2N1NT IN3 IN« c //,2x,frI5,:nc«R,!P12.!.^Ili:> 
WAI7Ett>,l«li;L:T,29iri,ANORTH,ASOUTH,ACASr 9 AwrST,NOASHO t KOASHU,XLAB£ 
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1S.« ?LtSWCftTH.IR€CCY 

1«9 14 FCR>"T</.SX t 7CH?lIT 2BIG AfcORTN ASCUTH AEAST AWEST *1)A$H0 NOASMU X 

UC 2*.*«5i SMOOTH XRECCV,// t 2X 9 2rX2«3*9FlQ.9«2I6t2FlQ«3 f r6»/l 

161 MAXC F0S3 

162 bHAT:-99. 

163 COM^CSCOWTNT 
16« IFICCNINC.^E.O.IGO tO 3 

us _ w»itci6»2i: 



tf»ITCl6t 
FCRHATt/ 



66 2 FCR*'rt/,2X,3H*AP t l3,X4H 2£«0 INTERVAL* 



ib« J IREClrlPSCCYU 



r writ *iu | 

16° LCASfUNOASMO 

17T LOAS^INOASHU 

17! OASMfSr. FALSE* 

17? IF<LCASH2.GT.O.ANO*L0ASHl*GT«O19ASHER£«TRUe. 

173 HINU»?KLABCL 

\l* ?FlHlNU«!GV!EriOOLA?$r.TPO€. 

176 XF<CCNINC.LT.0.IC0«IKC=*C0*INC 

177 PVCL:.0Q5*C0NXNC 
XT HOSUCra 

179 ViLL 1*1-989*9*, 969 

18C NUVI:NCX2-NCX1*1 

l«i NUVT tNEY2-wrri*l 

18? IF!NLVX.GT.J.AfcO.NUVX.LE.IK3.ANn.NUVY.GT.3.ANO.HUVY.LE.IN9lGO TO I 

18« 60 TC 120 

lCA C SKIP IF MUVX OR NUVY LESS THAN 4 

127 8 YCtfThrHl-ANORTM 

18* SCUTl-rASOUTH 

18^ £AST:wiO-AEAST 

19? WCSTZAUEST 

19? iF(Eisr.GT.wJniCASTTfc«ro 

19* IF<$CUTH.LT.C.I$OUTMSG. 

1«9 IFircRTH.GT,HI)YOfflH=Hl 

1«S HIxEMWjNU* 

196 WCOErCAST-WEST 

197 HOCMrYCRTH-SQUTH 
19» XXLA J7i99. 

199 YYLA5T=99. 

ZCr OlNCtPLTIIWC 

211 CL1T:01NC/U99 

2C2 CBXG:eiNC/2. 

20! T*AX:*.0«IY0PTN-S3U*H*EAST«>WEST> 

JC« XCRIC=wrO/FLOATUUVX-ll 

2C< YGftlCrHX/FLOATlNUVY-il 

2C6 HlNCtXG^IO 

2C7 IFIYCPXO,LT.XGRIOIHIWC=YCRIO 

?:• XrSHCOTH 

?:• X M X, L T., 25. OP. X ,57.7.51X2 1.0 

21? HlfcC:y*MlfcC 

2 11 CUIOFrS3»TIXGRI0*XGRXO*YGR10«^GPl01*.GX 

21' CL05 !T=.Q« 

2M C CL0S1T IS VALUE FOR CLOSED COKTflUR CH£CK 

219 NfXl:NEXl 

2t c NPY1 :N£Y1 

216 NHX1 1SMWX1-I 

217 XPY1 1TNHY1-1 

2ii *;x4:*exim 

210 n£Y4:*CY1-1 

22? c 

22! c nr it nrTCPPimc hax^ako "X* values Ih scalar fielo 

22? 2»AX :HM(N£Xl f N£Yn 

22! 2^IN:2max 

22* 03 3C lZNfXl t N£X2 

22« 00 JC J=N£Y1.MV2 

22f XF(M^rr t j)«f ? f.;^Axi7HAx:iiMcr f ji 

22 7 If <HMI.J|.lT.r*IN>2*r*=HHir»J» 
2:« 30 CONTINUE 



IF<2*4X.GT.20I< l2*AYr2eiG 
Xri?MK.LT.nJTI2''I«:2l>IT 
N?*T ft£TKP*!IM£ BOTfOH ST 

pval:sa«cos 

32 IF|Pm.Gt.2PlNIG0 TO J* 

PVAL :»V/a*COKl*C 

GO TC 1? 
39 XFIPWL~CG'J1NC.LT«2 W IN)G0 TO 35 



, . ..ARTING V*LUE FOR CONTOUR 10C» 

pval:sa«c5s 
XF|P%AL.Gt.2P!WIG0 TO J* 



»23 



n 



22? PVAL^PVAL-CONINC 

u c G0 TC * 3 * 

9C * 35 XPPrC. 

2*il YPP=C 

2*»J NUNIVX-1 

f*3 N2SNLVT-1 

:*• c 

215 C CCNTCU* LOOP STARTS BELOU AT STATEMENT 36 

2*)6 C THIS LOOP DETERMINES WHCPt TO START A NEW CONTOUR. THEN CALLS 



:*y 



C THIS LOOP DETERMINES WHCPt TO START A NEW CONTOUR. THEN CALLS 

t all u ?Sn°V^I HSMM*."" " CM C0,iTOU *' em li *' h UH * H 



}*« c 

25? c tkre arc 2 scams for each contour value* ftrst with va»xa*le outs as 

251 C FALSE SELECTS ONLY CONTOURS ENTERING GRID FRO* OUTSXOE EDGES* 



25*1 C 

255 C 

25ft 36 XF(PtAL.GE.?*AX160 TO 110 

257 OUTS:.FALSC. 

25* 00 37 1=1, Nl 

259 00 37 Jri»N2 

26r 37 RGRICft.JI^O 

261 3* 02=9*9999. 

262 00 ICO 1 = 1, Ml 
26? 00 ICO J=1»N2 

26*1 IFIOITSIGO TO 6?0 

265 irtI.EQ.l,aR.J,eO*1.0R.r.E0.N1.0R.J*EO«N2>GO TO 600 

267 60C inft£*l£fItJl.tQ.lIGO 10 ICC 

26! ir«RCRI0t:»j!,GT.UJN0.OUTSlGO TO 100 

269 II=Ne*«*I 

27" JJ=NEY9*J 

271 HEN( n=MMfIT f JJI 

27? MSN! «l=MMf XI,JJ«1> 

27? HEM !lSHH(XX*ltJJ*ll 

27* H£M*>=HHtIIM ,JJI 

27* 00 *C0 *C = 1 ,M 

27* XF(A85(MEN(KI*PVAL).GC.PV0L>G0 TO 900 

2 77 IFCHCNfKI .GE.PVAL IHCNtK I =*»VIL *PV0L 

2 7* XFtHCNItf >*LT.PVAL!HeMK>:PVAL<*PVOL 

27« 900 CONTINUE 

Ztr IF401TSIGO TO 250 

261 NCNNU 

2*2 IF(X .e0.1.AS0.HEN(ll«GT.PVAL.AN0*NEN(2)*LT.»VALlG0 TO 6*U 

2*1 NENN * T 

2M XFlX^0.Nl.AN0«MEN(3).GT.PVAL«AN0.hEN(*ll.LT.PVAL»G0 TO 601 

235 NENN:9 

266 XFtJ,fQ.l,AN0.HENr*)>,GT.PVAL«AN0*HEN(lJ«LT.PVALI60 TO 601 
2S7 HtHHzZ 
2*f IFIJ.EO.N 
28* GO TC 602 
2* n 25C 00 910 KTl,«j 

2*? li:nO 

297 iriK .£C.9tX2=l 

29* IftHeNdl ) .CT.PVAL.AN0.H£N1121 .LT.PVALJGO TO «G« 

2*« 910 CONTINUE 

2*ft GO TC 622 

2*7 t»rs NENNCK 

29f oOl IFIR(PI0(1»J).EQ*3«0R«CUTS)G0 TO 6*0 

21* I t=RCOX0(I .JIMC 

3oe I2=rcho<i.j»-io*xi * . 

JC1 iFIJWEO.f.CN'M.Oft.IJ.EO.NJ^NIGO TO ICC 

K* 6*0 GO Hi 3«1,l*l2 f 3*9 1 3*6> ,NENN 

JC 602 IFIRC*IO!I,Jl.CC.O>°GRID<ItJI=l 

3C9 GO TC ICO 

XC r 390 r:YGCIO*|FL0«T(J-ll«4PVAL*HENIlll/(HCN(2)*NCNtlIll 



ICM CASE. 



2"? XFlJ.f0.N2«AN0*HEN(2l # 6T,PVAL«AN0.HEN(3).LT.PVAl)G0 TO 601 



Kft xrxjtlO^FLOAtll-ll 



2C7 GO TC 95 

3C» 36 2 irXGRT0MrLCATU«lJ#<PVAL-HENI2lf/lHEN<3)««ENI2lll 

K9 Y:rGC70*FLOAT(J| 

Z\~ GO TC H5 

111 399 Y:YGFTC*(FL0ATfJ«lM<PVAC*HENr*lM/(NCN<3)-NCNl9lll 

31? X=XGCXC«F10AT(1 J 

Jl! GO TC *«5 

319 396 XIXGPIU* irL'fATiT-l 1 *iP MAL'HENt 1 1 1/ iHENt*! *HENt 1 1 1 I 

315 YIYGSXO'FLOATC j-11 
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316 63 oz<%**Ppfin«KPP)0tr-*p*fir»*P*i 

317 IMO.nt.DJIGO TO ICO 

3n NfNT!»SMCNN 

3J1 LCLXTl 

m mi 4 

IV TTTst 

$?• 103 CONflNUC 

32« iri02.GC.6M960.>«0 TO 109 _ 

324 IF<RCRIOUCLX 9 lCLV»«£0*0>R€fflO<LCLX t LCLY)Sl 

*!: c hwi 

]i: U*ITrt6»101IPVALtOOLAl$tOUTS 

H? c lal ***cxl*{ALL*!fiMouJifcC co*lw to actually oraw contour with value *val 

33; CALL CO*LIMHH,lNl,l*itRG*lO»t*JfIM) 

33* C NOV 30 «AC« TO IHkt* LOOR TO SEE IF THERE ARE OTHER P¥AL CONTOURS 

$5? C TO «E ORAWN. 

33? C 

32* GC TC 30 

3*« 1QS IFCOITSIGO TO 612 

3«* OUTS:*TRUE. 

3«2 GO TC 39 

30 612 PVAL:RVAL*CONTNC 

3«* € INCREMENT CONTOUR AM) GO TO TOR Of LOOR RO* NEXT CONTOUR 

3M GC T( 36 

3*6 11G CALL RL0T(0..0.**3t 

5*»T c 

3«« € NOfc TO ORAM THE OUTLINE OF* INTEREST AREA BT CALLING OUTtlN 

3*»« C 

J! - ! AUD»**SCALC/*. 

351 CALL RL0T|-A1,-AI,-3I 

35? CALL OUTLl«UNtJN,KN f NCNfOXtOY«XSCALC*VSCALE» 

:<? ccosr t snir t tx t tTi 

35* CALL PLOT* Al,Alt-3) 

3S« . iRARMHtR*! 

356 I»CCCYSt»ECCtM 

357 yPITf ft,ll5)I^*P,rPFCl »I»ECCY 

35* 115 FCR*lTC/,iax,nMCONTOUR NAR.I3 t 26M BEGINS MTTM PLOT RCCGRO ff 13vl6H 

Si! 2 SS?Tf?2h?2J^!iVc,miiH. N *»c fi r.wM*i ■ 

J»T 120 fiCIUCN 

J6« CM) 






*«*& 

** ^ 
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»riO«Hl! 



.conu* en e»c»tco o« « »?« T ?„! T Jii2*!2t 



IWS8 f lS»ii»itin£S!;i;ii;!fii!^8liKsiz«Sifi 

EdSftj" 



r:^i8:r64!6:cuTor.s«Ni.soMi.TMi.^» 

WWawiHC'r. 

LOGICAL I*C$«OOtA»$»OA$MC»tCLO$tOUT$ t OA$H!X 



,XX»L0Tf2?St2ltMAXC^I f UXC(41t 



THIS SUMOUM* IS CALICO TO OftAU CAC" XNOmWAL CONTOUR 
ir OOLAJS CNU«$ A$ TllUt, tAOCL COHTOU«S WITH MtUHT HXMI" 



0«SHlX£OASHtH 



AOLXTSQ 



YMX:-f • 
xr*x:»f« 



WttoSUMCMTC* 
lOPLCTrJ 

*HA*C=LOASMt 



MS0ri=tOA$M2 
NUGGSO 

xxrx 

YYSY 



xnic-xx 
tftcm 
LZX: 



I»ZX = ICLX 

tctsc. 

NCO*C=G 
CCOS^MLSC. 
60 TC 100 



CNO SCW . «0IH LOO* THAT PICKS CXACT STRIXGMT LUC SEt'CMCO T*AVC«SE 



250 IM*(0«O.IT.»COI60 TO 252 
6C TC T30„ . 



COITOUft 



252 NC&»tTKCO«C»l 

MTPTtTrWYPTflf*MYP* 
Ml^!«»lC050iSMYt{ 
CORO tNCOjOa UXXSO 
CC«0(NCO»0,Jl=TtJO 



701 



2*CUT!.AS0 I*CO»O.OT.J|CO TO TCI 

GO TC ?«53. 



m«XO*YO«YOI.tC*CLO$IT»A«0 



mm*- 



oka9l;ts*i 



700 \\WViY.lSl?,H ( HU***tXtl.t 



»flc:!I?0 



K SS!:k!:i<JK oiir»c 5 to n 
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iMon.mo.iotuo. 

If (OM»OT,ro»IO»OTi:TO«IO 

f:i?»*««xii 

StAUT CUT tOimt 100* 
MAXf tUMHII, Jl 

NAXt* 

H 




%C1 



ax< !>:MHCT*l t J*l! 

AI( MIMM! f ♦$♦ Jt 
OC *Ct IUIiI 
irtAfS<HAXtXX)-'»VAlLQC«PVOltGO TC 401 

in* #itm.Gc.pv*uNAK<:ii:»v«i**vcL 

inHIVCfll*LT*MACIMAXiillS^VAt-MCL 
CONT 1MUC 



■Jiciw 



00 4j« Mist,* 

XMXIX*C0.4IX2:1 
C STIUKCNT ICIOU SCICCTS SIOCS THAT HAVC CUT JOINTS 

XnMJXtTlKLr.*m«A*0**AXtT2t«GT«MAUC0 TO 420 

60 TC 4JS 
4^C *UH0lTtr4UK0UT«l 

ir<Ni-out*tc,t »w*i=Tii 

in*l»OUT.C0.2l«N2:TXX 
*22 i?2TJ|JvA(!*MjMUf/*MiiJiJ-HAllin>*TCtT0 



c 



) X*2 f * 
>0 TC *TC 
0X2t HPViL«MAX4ZH/<HAK<3>*tUXt2lt1 

CO TC *t!S 



424 OT2:<fPVtl-MAX<*M/<HAX<3l-MAXt*t»l«VGfllO 

co tc *:a 

421 2x2;f(^VAC*MAXflll/fHAXt«I*MAX(lt!|«XCAX0 

Ov2sC* 
*3C CIOC <!IX,L1:ex2 

cioi 

435 CCK 



CIOC «TXXt2lsCTZ 

C Pk T IMLc 

UUCiS WC HAVC MILL POINT SOU Aft C, NUAOUT IHUO OC 1 HXTH OUT AT 012* Of 2 

xn«i«our,NC*i)eo to «22 



412 ?MNt»OUf .C0.2IGO TO 4J1 

4X1 tFtTeU^tlirxaZY.NUMOUT.PVAl.XOlGjTMG 

*3* roaM«r</.2x t icHMo mat ouT t sx t 3iio«3rid*2t/i 

GC TC 503 

C tCGI* SCCTION IHAT CCH»«IN?$ P»0'C» MT'» TH»U 0«XQ $00#»C 

C COM'lftlMG HYPCMOtXC CCHf IGU«4T2C«U C2 CMTfff AJlO 2 C»XT JI0CS3 

439 Xrf»C*lOUrx f L2YI.GT*t)eO TO 442 
XlO:CTOCt*M, 11-0X1 

YlDSCtOC<K*1.2l«0*i 

AA7!9ArurO«XXC«vlD4TXOl 

tiozcfoctN*2»2i-o*i * . 

C»»r!0«T<xl*«xxr«*ic«f XOt 
iriOI9.LT*OAAIGO TO ««C 
43* 0*2rC:DC<*i«!*l I 

GC TC «1« 
44C 0»}I(lOr (NN7.1 I 
Ot2tftoeiS»«2 f 2» 
NfU Ut*2 

414 *G»ita7x.tZYt:io«Nrft»TeiM«icxtT 

GC TC *»** 

4*2 it:fff«xc«i;f 9 i{Yi/ir 

ir»ij.of.c.»ho.i2.cT.a.*NC.n.M.ii.*»>o.»i?»itt«uT.cico to *it 
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W 



20 



ii; 



:« HSKmjMlttBrtJ'di'HiKJMiiii 

•0 TC $30 
417 LCir MISO 

lilt 1*1 *3 

text tmsi 

lFtLfuni}.^o.e.ANo.n.MC.NCMTft.ANo»cioc<iitii.ttT.<-«$j.A«o*ctoc 
411 coNtr»uc_ 
414 



•9!SUUIH:1I 



:«o myikmouc emo scume see no* 



ir<ffC»I0U?X t l2Y)«C0«0 1*6*10 Cl?X,LZY I si 



XlfcO:OX?«»GfHp«;Eo*T(LZ*<»1> 
TlN0:ftrJ*TC»f5»rL0*T(LiT-{l 

iriffc»io(L?x t c2r}«co«o 

ixso3»ino-i4|c 

HY#?:S0*t7xxS§*XXSQ*YY_ 

M^ITJU.SfTttJx.L^t.ldCNTln.^CXCtttllo.TirOtXIIOfTtlllOfOtl.Oflt 
?OTi 9 HSO«rV$0 9 MV#r 9 NCN$T.L(LX t CCLY 9 MCOR0tXXtVV 

344 GO TC |i§ 

ClNC SCG*C*TCO CONTOU* T9AV?*SC no* CO'FteTf, NCXT. OTVfOf TMTS 
fftAtfftSC T*T0 MI1C CCUAL StGrCHTl. THIS M«4C* IS rUHCTfOH NOT 



Oftiv Or ir«kOT 
(OKTtOli OC 



l»Al !5*!C*fr»_*uT *r Xkconiitft aroumCNT smooth, which 



rate or shoctnjho 01 



ICO* 



730 IHNCO*0.l!.M*XC*OI«0 TO 732 
l«*r*C04O 
t*VAL 



HMC 1 



XH.»0UCl0$.'*O*4CeXO*0T.l.AM0,N»lftC*6T.l)CO TO 734 
lF4CtGS.A*0**COftO«ftT«3.A«0«HffI*C«GT*2>€0 TO 774 
00 TC SO? 

ncxt, sct up CHttY »*o tin sto*e cata ro* ri»$T sc«»CNr icro»c 

CfTCdlHO nAlH CUfttlllfltCA* XHTCft^OCATC A*0 4tOT COOP* 

T34 x«{$:ix 



Villi 

Xr(4rTMC«HANC*LT ## 75lOASHlX;*rAtSC« 

HSO. 

9{Hi 

00 ?< 



IX 

•0 Y lrl t *lC040 
«H!P»SC:J«GC«MAWCIG0 TO 742 

nrr^s'M 



743 

74: 





iriH 

ro re 7*5 

x;im "" ■ * - 

X'.M) 



:v£fcO*'**C**C<Tftl 



74$ X 



in** 

7*7 W»IT 

?*• ro«* 

02 r 

?sa cmo 

SA4C 

HYP* 

M 

? 

xt:i 

TOST 



tYtio*t*e« 



751 



11 



ictrorxsa*fx;Q«YY:o*rr9Qi 

4C.0C*.Cbwi M.0 TO 753 
tU.7*»«|f4tG,f»r0.PYAL f HY»C 
*T</.jX f 1«MHY4£ TOO S*All t 2X»4rU # 4) 

:XX5C/MYr€ 
tVYSC/HY^e 

tOSlSO TO 75? 
iCrCA;*r» 

C 7*4 

> 



t28 



DO*??! Isl.ttCOftO 

j:*COO-t*l 

ir<M*Ml»SU».CC»«A«CUO TO 7$$ 

ip=M-co*ou f n 

7$J V*:tP-C0ftOtJ 9 2t 
00 TC 7S7 

xt:xF~x«cOff?u v h 

T^=YP-*»COtO«J t }| 

>YT50t 



757 



XKS0:K9C6*M 

ccsbic'h: 

j:.}«aT*Sj 

SASf INIM 

OsCCSMt 
$ 

i 



lI*l*ClTTSC/WY*t 
CC$BIC"I«SO/MT^r 
7S9 A:/}«*riS2<S4ft6«C0$IAC*CA«6*$TftMC»CA*0*C0SSAC«UNa«Sl*tAC) 



SEftT:$X*8AC»C4«C0tt«C*SA 

~* B!5 " IM -tttt5tfaaW M 



KTItSAtAttJ 
CC£*1:C£*T 

c ewttn "aim cu*mi*c** mrff*oiATc a*o no? loop 

00 IC3 tUPtitt*M*e 
XMU»£*NC*M*TftC>00 TO 7*2 
IMClollOO TO 7*0 

)OUT:$a*G 

ccur:CAK6 

00 TC 209 
74Q SOU?:iSC*T 

CCUTCCCMT 

CO T( iJo 
7*2 XI*0:XA 

f l*0:rr 

2A*C :HA*C«riOA T ICUHM I 

00 Tt* Irl,^CO»0 
lMM«»l»'ittl«Gr*2ANCI60 TO 74* 



"""" " T 1 2 I 



7*« vtHQ:viNn*cotO<t( 

OO TC TM 
7** X:i2**C-Ml/MirP5(TI 

X|«0:»l*0*«*COffO<T t I! 

TlM6:tpo«»«COOOlX«?l 
7*1 X«SO:*l»iO-»{*e 

tysc:yimd-tC*o 

H**t :$C»T|XXSO«XXSC«Yf $0»*YSOt 

IM»n»C*lT.*CC01IGO TO 7*7 

S!fcrt»rr¥$o/Mv*j 
cosrct sxx$3/my»C 

MT^f(»SHY^t 

A:.S<lTA«2C$l*r0ft«CA*O*COSF0ft«SA««fCQSfOft«CA«O«Sl*t:^«S«*G] 
$*=$ t*«AI 

$CuT:{a*G*CA«CA'C«$a 
CCuT:CA«G«:A-SAtG*S' 
CKT :iT*N2IS0UltC0uTl 

200 mpr*YH 



Irixc^(5.Lf*«csT«or 



S**XCNO.0T«rAST.Off*YCMO«tT.SOUTH«0t,rCNO«0T*rAftrHl 
?€0 TC 7*0 
I'ftKSKO TO *** 



2*LL »CoV^etC»WCST # t|CC-SCUTHt!l 

*** Isj#e» 

!'<OMHtX.A«0,X«CQ«?ll:tOPL9T 
CHI 'L3M"„ft0«WC$1*TC*C«S0tJTM,X> 



00 TC 7*3 

1C6Z* $***C X*TC**CUTXOK fOO S€G»t*T 
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316 


*• 


3* 7 


C 


I* 


c 


c 


IK 


c 


c 


32? 


c 


in 


I 


32^ 




316 




?!- 




32* 




330 




331 






33? 




333 




33« 




33« 




336 
33? 






32* 




32* 




SUP 




K! 




3«3 




J** 




3*«5 




3*»6 




3*7 


c 


3»» 


c 


I*'' 


c 


"JSC 


c 


351 




25" 




25* 




35*i 


c 


355 


c 


356 


c 


J57 


c 


3?o 




36C 




3b! 




362 




363 


c 


3fc« 


c 


365 


c 


366 




36' 




263 
$6* 






37" 




371 


c 


H? 


c 


c 


37** 


c 


37« 




276 




i 77 
37* 






370 




38^ 




381 




i*3 




3M 




3a« 




336 




3«7 




J 2 ' 

18* 






39? 




39! 




39? 




3*! 




39* 





ENTIR 4*0 ExET ARE TRUE ENTRY AM) EXIT DIRECTIONS AT ENOS Or SEGMENT 
CSMKELIKE INTERPOLATION IS H'OC FOR FIT BETWEEN 2 END 0IRCCTI0NS3 

SANG A NO CANG ARE SINE ANO COSINE Of STRAIGHT LIME ENTRY/EXIT CONNECTOR 
SCGXft SECTION TMAT INTERPOLATES ANO PLOTS THRU SEGMENT 

*t%9 SrSE»T*CAN6-CENT*SANG 

C:C£»T*CANG*$EkT«SAMG 

Cl-C J.»ATAN2CS,C1)>HTP 

SlN r i:SANR«C-CANG«S 

CTNi>:CASG«C*SANG*S 

SsS£ir«CXNG-COUT«$XKfi 

t;cCtT«CXNG*S0UT«SlN6 

C2 = < <«*ATAN2tS f Cl WHYP 

fxrXKC 

TYxYFCG 

NINC:C 

TYPr.CBIG 

HYP*IXrMYP-CLIT # 

H2S=.25*MYP 
*5Q TYP=1YP*0INC 

OisTYP 

iriQl.GT«H25I0i:H2S 

02-T1P-M25 

IF<0<.LT.C.>02=0* 

SlN6: r NTlR-Cl*TYP*C2*(02-0l* 

TX=T>«OlNC«COSCSI*G» 

TY=TY*QINC«SIN(SING> 

EKO SNAKE INTERPOLATION SECTION — TRY ANO FIGURE IT OUT ANO GO NUTS- 
NEXT STORE POINTS THRU THIS SEGMENT FOR FINAL AOJUST ANO PLOT 

IFINIMC.LT.275 JGO TO »53 

WfiITC(b,*»5*i)PVAL.XBCG»rBECt*EN0tYEN0 
*5* FCRH*T|/ 1 2Y ? i?hkj\x SHUTCFF t 2**5F 12 »3 I 

IF SETOFF MESSAGE °tCEIVED HERE, ARRAY XXPLOT IS TOO SHALL. 

FOR l^FRECUENT KCSSAQES* OOMT WORRY ABOUT IT, SINCE LACK OF 

CLOSURE IS ACJUSTEO OUT* IF MESSAGE PERSISTS. EITHER INCREASE 
•LOT INCREMENT LENGTH OR SIZE OF XXPLOT. 

GC TC H$5 
*IS3 NlNC:NI>lC« t 

XXPLCTCNlNCf 1I7TX 

XXPLCTC*IIMC,2>-TY 

IF«TTP.LE.HYP*AXJOO TO *50 

AOJIST F0» CLOSURE ERROR- YHFN PLOT CURVE ALONG THIS SEGMENT. 

H55 X?R=U£NO-XXPL0TCNI*C f 1 1) /FLOA T (NIKC > 
Y5fi = (YENO*XXPLOT(NINC««n/FLOAT(NlNC) 

U = u. 

V=(i. 

NUNC=C ♦ • 



^t^lH SEGMENT PLOY LCOP -OASHEO OR SOLIO CURVES— •• 

SUFPOUTINE EN0E9 IS CALLED TO LABEL LINES 

DC 610 HUNINC 
UrU*>ER ■ 

v=v* >rR 

X:XXFLOT(I t l |*U 
YtxxFL0T(l,2)»V 
!F<X .LT.WtST.0R.X«GT»EAST#OR.Y.LT,S0UTH.OR.Y.GT.Y0RTHJG0 TO 60* 

IFtl^fSIGO TO 603 

INCS:«rRUE. 

1F(IFCR.C0.2)CALL PLOT (XBEG-UEST , Y8EG-S0UTH , 3 1 
603 XPPrx 

YFPCY 

XUXr XPP-WEST 

YUYT>0P-$OUTM 

IF CL /"tTT.GT.Oir.O TO 60* 

IF(L J?Lf T.E0.C1CALL ENGER I XUX t YUY f PVAL # 1 J 

LA3L IT x I 
bG*» IF(0'<HXX.A^O.IPE°.EO.«!GO TO 615 

C*LL "LOT(XUX f Yl/Y,I«>ER> 

GC TC 6^9 
blS CALL PLOTfXUX,YUY,ICPLOT) 
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39« IMICPLOT. £0.3)60 TO 430 

39* n*apc"nma RO-l 

397 IFtNMRO. fit, 0100 TO 40* 

39» ICPLCTrJ 

399 NmaJ»C*LOASMI 

ICO GO TC 439 

901 430 N$ori=NSorr«i 

952 ir(N«orT,ttT,otco to 409 

9C3 fCPLCrr* 

9C* NSCM=L0ASM2 

9C5 4C« IP£»:2 

*!? in.»0t.0QlAes«Oft«YUV,L£*YfUXlGO TO 41C 

9C! YHAXtYUV 

9C« X«AX:*UX 

Mi: GO TC 410 

911 404 IP£PM 

91? iFlL fPLlLNCDGO TO 410 

9 1? CALL £NOCRtX-*€ST»Y*SOUTH,PVALtl) 

919 LA3L lf=0 

91? 410 CQfcMNUE 



916 
91? 



C E>0 AOJUST ANO PLOT SECTION 

91! C 

* 19 NUGG:NUGS«*UNC 

9 2? T0ItT0T«01NC»FL0AT<MJNC) 

921 XF(TCT.LT.TWAXlCO TO 790 

9 2? MfilTtf6t«4?>PVU»TPJX 

92* 942 FCRW/T</,2X t 14Mff£ACHC0 TMAX ON t2Fl?.9» 

92* GO TC SCO 

925 79Q IFULPt.CO.NfttNOGO TO 0QQ 

9Cf MYPH:my?F0I? 

927 StNT :SOUT 

9Z« CENT:COUT 

929 CNTIKzCXfT 

9 3? $AfcG:$lNrC9 

93! CANGreOSFG* 

9!? Xfl£G:X£N0 

92! YC£G:Y£NO 

93* X£*.0:XIN0 

935 Y£\0:riNO 
800 CONTINUE 

93! C C*0 NAIN CURVILIHAR INT£»POLAT£ ANO PLOT LOOP 

*» J« C 

9«? SCO IF(. ►or.CLOS t Oft..MOT.OCLABS#OP,THAX.Lr.,atiCO TO 501 

«*«1 %PPz >*4X *WE$T 

MM"? YPPT>"*X«SCUTM 

«*9! CALL FN0CR<Y-AX 1 Y'«AY t PVALf2l 

9«ii 501 IFt.*OT,ClOS.ANO.LA«LlT.£0.1>CALL £N0£P tXUX t YUY»PVAL» 1 1 

*»9? IF(NlSG.L£.*0SINC)6O TO 502 

9** HOSI^CtNUGG 

9*7 VALL JNSPVAL 

94« 5C2 »£TU*N 

999 £N0 



9jfr 
93? 
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*'t»OUUI.£NO£P ELT CRCATEO ON 8 **Y 79 AT 10 1 17 151 

t SlBftCUTINE EK0C*<X,Y t PVAL t IC01 

-» CCM^CN /CENOEO/MONU*,MCOE t HOGM»*LAS f YLAS 

7 OIMCKSIO* 0131 

• C TUS SU8»0UTTNE IS CALLEO TO LABEL CONTOURS 

7 oh:sc?ttcx«xlas>«fx*xla$)*(y*yla$>*(y*ylas>) 

* ifio*.lt.2.«monuh.o*«honu«»le.g»1go to 25 
« if<ico«cq.?)go 10 1* 

K 0'i> uescY-MonHi 

11 Om:l§SlX«UOOE) 

It QtJI=A8$<Y> 

l? a = i 

I« 0*=AfSCXI 

1? 03 1C U1.3 

14 <1F<OU>.G£»OMGO TO 10 

17 0"=0!T) 

1« H:P1 

19 1C CONTIMUE 

2? GO TCtl2,l*.16 f ia> t K 

2! 12 YACr-N0NU«/2. 

22 1 = 1 

23 60 TC ICO 
2* 1* YA0=.92 
2* 1=2 

;* GO TC 100 

27 16 XAOS .02 

2« YA0r*HaNUH/2. 

I* GO TC 20 

:c ia y*o=*monuh-.o2 

31 i=: 

3? GO TC ICC 

3T 20 X10=>*XA0 

3« YAOst*YAO 

3 C CALL NU»8CP<X*C»YA0«MOMiM t PVALf3*t*21 

36 CALL PLOTlX t Y,;i 

37 XLA$:» 
f* VLAS:y 
JO 25 »CTUPN 

«? 10Q XA0=*.75«HCNUM 

• 1 inPVAL.Ge.9.5.0R.PVAL.LE.C-.5IIXAGrXA0-M0MIH 
*•; IFI?VAL.GC*99.S.0R*PVAL*LE.I*9.SI I XAO =XAO-HONU* 

• T lF(PVAL.G£.999.5.0R.(»VAL.L£.i»99.SnXAO=XAO-HONUI1 
«« IF 1 1 .gC»2)XA0?,5*XA0 

*»5 GO TC 20 

M6 EhO 
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•rtOtftn. r 

! 

* C* 



1' 
11 



IT *g* CPfATCD ON * 0€C 79 AT 10l02»3*l 

FITS A PARABOLA TO TH0?C POINTS* USED IN »AKlNf? PLOT. 

SOBffCUTTNC riT(21»5? t Z3.Ei 1 !T2 1 £T 1 A,e l CI 

= l2l*»2>M22-2Jl-21«W2**2-2 3»*2>»2$»l22*«2l-22*<Z3««2) 
ASIC i»<22-2S)-!l«tE2-C3MC:*2T-^5*/21/0 

BM<21»«2l«<e2-rj>-ei*t22*«WJ«*2l*ej»CZ2»«2l-C2»l?3»*2M/0 
CM(n**2IMZ2»Ci-Z3*e£Wl«<E3*(22«*Z)-£?*<Z3**2ll«Cl*tZ3*<Z2**2) 
C-J2«J?3"2lll/0 

£NO 
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t £•*•** TO OPAb VCLOCTTT VECTOR dV CALLING CALCOM^ $UpPOUTIN£ A*OwO»***« 

SU0RCUTI*C VCCT|AT»AAX 9 AJ,AAJtAf>NlN»AAH*X»CQSFiSINF> 
tW = . ;*SCftT* f II 1-4 II ••2MA1J-4 J J •»£) 

bi=a i*c5sf«aj*Stnf 

BJ = -«T«<;fNF*AJ«CO$F 
BBI:'AI*CO$F«AAJ*SI»F 
8eJ=**Al*SUF»AAj*CnSF 
CALL A*0HQ<BI ff BJ,9BI.86U»YU,C.»12> 
IP ACTUffN 
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♦ruoaiii.CAPTiu tl^ 



\ cli c»£at*d o* i« sc* i* »t 10:2**2* 

SUBROUTINE CAPT*1 IP|, PJ.P J, P«.P?.P6,C 7 # OBI 



? CMM*«wAttC CQ**QH MCAPt^C FOR CACM »UOT yHCTMCP TCHP Oft VC10CXTT««**«» 

• CALL SrMPCL(a.*,6.7CO.l0.2?HTI»CIJUNC 20. l*78»t .G.0.251 
■ CALL NUMB£PM.C,6.7r t 0.ia.P1.0.r>.*j> tw»u#«i 

• CiCL SYPJOLCO. 8,6.5^.0. 10, 25HWJ20 SPCtO<C^/$tCls .0.0,25 

• CALL Sfrtp3LCG.«,fc.So.a.l3»25auX«0 0Xft£CTXON(0Cfi/N|{ ,0.0.25 
" CALL ttUH9C*f««a f b.3%O.lO.P?,0.*.«O) ' ' 

}5 CA <-L 5Y*BQLl3.8,b.in»0.13t2**AIP T£»PCRATUPt fOCG-C) t .0.0,25 

)i s;u *u*3CPc«.c,6.io,o.i3tP».o.'> 1 M* %«•*•« 

}*, CALL SrwgOLCO.I,S.4o,C.lO,2 c .HfifsCHA»G£ T£PP <0CG-C 1 1 .0.0.25 

\l CALL *Uwe£P«ii.e»5.*O.0.ia.P5 t O.O.M! • f 

}; CALL SYMeOL<0.*,5.7C.0.lC.2SMOISCM TLOWPATC (CUH/StCl * .O.C.25 

}\ CALL 5UM8fftf* # 0,5.70 f O.J0,P6,O.O.*n • *.« 

}* CALL S7H58L(C,« t 5.50 f 0.10 f 2^NL£NdTH SCALtfK*: X C Mlt .0.0*25 

}J CALL VU«B£ftHi. 0,5.50, 0.10, 07,0, P.*CI ' '" 

17 CALL «UHB£ft<**Cf5.30.0.lO.G*»0.*i#*2l 

ji 8M UM 
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TLOWm.CAPTMZ SY > C&EATfO ON 11 SEP 79 AT 2?:29*52 



3 
4 

c 

6 
7 

« 

o 

1C 

11 



C*L 
C*L 
C«L 

CAL 
CAL 
CAL 



SUBROUTINE C*PtN2tNl 



t PLOT TITLE ON BOTTOH OF QTAGRAH FOR CAICULATEO UV***« 
SYMBOL' l.C.-l.A..l«»t2*»MfIG SURFACE VELOCITY, .0 ,2<t 1 
SY*flOLU.O t -?.2,.l*t29HANCLOTt ANCHORAGE BY MOOEL lNG f •0 f 29| 
PL0Tt-LC,-2*S,3> 



*ETUEN 
ENO 



PL0T*-UO f *7.5,2) 
PL0H«*.S,«7«$ 9 2) 
P10TI«6*S»-2.S»2) 



PL0T<-l*O*<»2«5t2> 
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•rLOfcUlKCAPT^J %VP C»CAT£0 0* 11 SCP 7* IT 2?t31rSS 
$U8ffCUTlN£ CAPTWJCKl 
C****ttftlTC PLOT TITLC ON PQTTOH OF OUGflAM ANO 9L0CK IT**** 

CALL $y«80LC1.0t*?.efl* t 29MFIO TCHPCPATURC f»OH Ift,C*0»29J 



IQ 



CALL PLOT(.UO,-2.S»J» 
CALL »LnTt*L0 9 «?.S f 2l 
CALL f>LCTl*6.5 9 *7.S 9 21 
CALL »L0T(*6«$ 9 *2*$t2> 

$tvh■^ of, • | • 8 ••*• , • l, 

CNO 
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•flOUtn.CAPU'tt SYC CffCATCO ON 11 StP 79 AT 22tZbtZ9 

1 StSftCUTtNt CA'TMISUHI 

2 €•••• WRITE OCVUnOK VAlUCtFffO* Iff) ON 1SOTHC** PLOT •*«•#•* 

♦ CILL SYMIOL <0*8,~0.5fl3,2JH0tVlATI0N FffOM Iff U*** 9 ««23l 

• CALL *Utti£*M«Ct«0*S«lQtSU«tC*fO> 

5 RCTU*N 

6 CM) 
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«KOUU)«CAM«; $Y* CRCATCO ON 11 UP 7? 

simncuTXHc cap TNS(N t nHcwi 



$Y* CRCATCO ON 11 SCP 79 AT 22t39r«9 
imUCUTXHC CAPTNSCN.hHCWI 
RITC TXOAL STAGE AT THAT TlHe««««« 



:8co;«xocwini m m m 

mraiitlS:t:::J:8:i8:a?M;a: .. 
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•fLOMULCAMN* S*' C&CATCC ON 11 %tP Tt AT 22s*ltSl 

J ^••••5S!T$ U JtO? flTLC^oS'llCTTO* OF OUftftAH fOK CALCUL ATCO ««•♦•• 

CALL StW*CLU.G#-t.*t.l«»l*MMO UW VtLOCITt, ♦0 1 i2 1 .^ „ mm% 

CALL $T*5oLflIOt-*.^t*l»t2»HA*CtOTt AKCHOGAGC if -liOCL I*G ♦ .0,2* I 

CALL *L0Tl-l«C,-2.} f JI 

CALL *L0TI-l.C f »?.},5t 

CALL >L0Tl^.S.O.$ t |l 

CALL »L0Tl*^.S 9 -2.)*2) 

CALL »10T<-I«0fl*2*3«2> 

ftCTUIEN 

CNO 
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c**«*«w*Irc *lot titlc on toiTo* or otao«am row calculatco *«•••• 

CALL ft«OOLU.Cfl.f».M*lf*f:« ¥. VCL0CITT f ♦O.lfl 

CALL $r«fi0L<l*Qt<-i«? t .lt*2tMANCL0U AftCMOftAO* IT «OOCLXftO**0«2t ' 

CALL »LOTf-t. 0,-2,5.11 

CALL ^L0TI*t«2t*T«i til 

CALL *LOTt •t*S»*T.9,2) 

CALL »L0n**.S»-2,St2! 

call »uort<»i«oM*s*2> 

ft£TUKN 



CNO 
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*;iOU<ll*CAPT»| 



CALL 
CALL 
CALL 
CALL 
CALL 



tri *LOT TITLK ON ftOVTO* Or OlAG*A* rot calcucatco clcvattom*««« 

'"*eLtlj^*-?.? t .l< v H»ANClQT£ AfcCHQAAtf ft «OOCL 1*4, *Ot*? I 

5«i!:i:?;;?:l:l! 



CALL ;iOT|.fc.j,.j.j,j| 



e*o 
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r 



-. ,t- rv**b?r*X<^* 



j-SWCMI***! 1 ^""'^T'Bfc-'VEH^t,: , 



* c*««««k<irc noT Tirtr on iottoi* of oiasmn re;* c*icui*fto Tr»»rMTu«« 
I SiSI- J T *5otn.c,-i.#,.i«,j7t«fi« iu»r*cc Tt»*c»»Tu»e..o,2Ti 

5 £?H I. T 5?? l ! l a" 5, : , : , i: l, «* ,MM,eww **cmo#»« by «©oeui»6,:c,*fi 

* iitt JbST T !:l:8:^:l:j 



It flO 
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